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ABSTRACT 

A phenomenological formalism is presented in which the apparent acceleration 
of the universe is generated by large-scale structure formation, thus eliminating 
the coincidence and magnitude fine-tuning problems of the Cosmological Con- 
stant in the Concordance Model, as well as potential instability issues with dy- 
namical Dark Energy. The observed acceleration results from the combined effect 
of innumerable local perturbations, due to individually virialized systems, over- 
lapping together in a smoothly- inhomogeneous adjustment of the FRW metric, in 
a process governed by the causal flow of inhomogeneity information outward from 
each clumped system. We discuss several arguments from the literature claim- 
ing to place sharp limits upon the strength of backreaction-related effects, and 
show why such arguments are not applicable in a physically realistic cosmologi- 
cal analysis. A selection of simply-parameterized models are presented, including 
several which are capable of fitting the luminosity distance data from Type la 
supernovae essentially as well as the best-fit flat ACDM model, without resort to 
Dark Energy, any modification to gravity, or a local void. Simultaneously, these 
models can reproduce measured cosmological parameters such as the age of the 
universe, the matter density required for spatial flatness, the present-day decel- 
eration parameter, and the angular scale of the Cosmic Microwave Background 
to within a reasonable proximity of their Concordance values. We conclude by 
considering potential observational signatures for distinguishing this cosmologi- 
cal formalism from ACDM or Dark Energy, as well as the possible long-term fate 
of such a universe with ever- spreading spheres of influence for its increasingly 
superposed perturbations. 

Subject headings: cosmological parameters — cosmology: theory — dark energy 
— large-scale structure of universe 
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INTRODUCTION AND MOTIVATIONS 



During the past decade or so, a number of complementary observational techniques have 
come together to point towards the existence of an exotic cosmic ingredie nt known as Dark 



Ener g y (DE). Luminosit y distance measurements from Type la supernovae (IPerlmutter et al. 



19991 ; iRiess et al.l Il998[ ) bear evidence of an apparent cosmic acceleration, which seems 



to indicate the existen ce of a substance capable of violating the strong e nergy condition 



(Hawki ng & El 



Caldwell et al. 



is 



200 



19731 ) . and perhaps even the dominant energy condition ( iCaldwelll 120021 , 
"^). O bservations of the Co smic Microwave Background (CMB) indicate 



a spatially flat universe (IKomatsu et al.ll2009l ); yet vario us lines of evid ence regarding large- 
scale structure formation and other measurements (e.g.. lTurnerll2002al ) indicate a total dy- 
namical matter content of only Qm ~ 0.3, far less than the value of unity needed for flatness. 
The standard conclusion from this, barring modifications to gravity, is that there must exist 
a smooth (non-clumping) component to the inventory of the universe supplying the remain- 
ing f2 DE ~ 0.7, which is capable of accelerating the universe without joining in the process 
of localized structure formation. 

Within the current resolution of the data, combinations of observations are still consis- 



tent (e.g., IKomatsu et al.ll2009l ) with the simplest Dark Energy model: flat ACDM with a 



Cosmological Constant (vacuum energy) as the Dark Energy 'substance'. A nonzero A den- 
sity would clearly avoid participating in spatial clustering, while also possessing an equation 
of state, w = wq = — 1, that would begin to drive the acceleration once Q\ had become the 
dominant cosmic component. 

There are well-known aesthetic problems with A as the missing piece to the puzzle, 
however. This includes a fine-tuning problem in which one must explain why p\ is some 
~ 120 orders of magnitude s maller than what would naively be expected from the Planck scale 



(e.g. 



Kolb fc Turnerlll990l ). And perhaps even more relevant to issu es of cosmic evolution, 
there is the "Coincidence Problem" (e.g.. lArkani-Hamed et al.ll2000l ). which questions why, 



given that p\/ Pm oc a 3 , that observers today should just happen to live right around the 
time (within a factor of ~2 in scale factor) when A began dominating the cosmic evolution. 
After all, if p\ had been just a little smaller, then cosmologists in this epoch would have had 
no hope of detecting it; and if p\ had been just a little larger, then the acceleration would 
have prevented structure formation from ever having occurred, and no observers would exist 
to measure it at all. 

Without delving deeply into anthropic issues, one possible way of avoiding such prob- 
lems is to give the Dark Energy an evolution in w(z), where such coi ncidences are elim 



inated or moderated; an example being tracker quintessence solutions ( jZlatev et al.l Il999 



Albrecht fc Skordid 120001 ). where there are natural reasons for Dark Energy to 'pop up' at 



-3 - 



a propitious, recent (relative to the Planck scale) time - such as being triggered by a phase 
transition like the onset of matter domination, for example. Moving away from a Cosmolog- 
ical Constant by giving the Dark Energy an evolving equation of state, however, also opens 
the door to other possible dynamics - su ch as the ability for the Dark Energy to be mo- 
bile and become spatially inhomogeneous ( iCaldwell et al.lll998l ). A dynamical Dark Energy 
(DDE) of this type is liable to join in structure formation to an observationally unacceptable 
degree, thus ruining it as a candidate to serve as the smoothly-distributed missing cosmic 
ingredient. 

An under-appreciated aspect of the notion of Dark Energy is what it really means when 
something is a 'negative pressure' substance. Thinking of the acceleration of the cosmic 
expansion, it is common t o refer to the DE with words like "repulsive", "antigravity" , and 
so on (e.g., lHeavensll2010f ). as if negative pressure substances were materials that naturally 



did not possess attractive forces or tend to clump on local scales, thus making them natural 
candidates for being a smooth cosmic component. But this is the exact opposite of the truth: 
a simple look at the 1 st Law of Thermodynamics, P = — d [Energy ]/d[ Volume], shows that 
for a negative pressure (P < 0) substance, an increase in system volume requires an increase 
in its energy content (dE/dV > 0). Consequently, since the environment must do work 
upon a local patch of DE to make it expand (and gets work back when letting it contract), 
a material with negative pressure should be self- attractive, not self- repulsive. Thus a true 
negative pressure species which is mobile (i.e., not A), able to move around and clump due to 
its own self-attractive forces, would naturally be expected to not only fail to remain smooth 
on galaxy cluster scales, but due to its extreme (w = P/p ~ — 1) negative pressure, it would 
clump relativistically, making it the most strongly clumping material in the universe! 

Being self-attractive would seem to be a counterintuitive property to have for a substance 
that is invoked in order to make the universe increase its expansion rate and thereby spread 
out faster, but in fac t it is due to a sign peculiarity in the FLRW acceleration equation (e.g., 
Kolb fc Turner1ll99oh : 

a AttG , „ „. 

- = =-(p + 3P , (1) 

a S 

where a(t) is the cosmic scale factor. From the negative sign on the right-hand side of 
this equation, we see that negative pressure does in fact cause a cosmic acceleration, not 
a deceleration; and thus we have the seemingly paradoxical (but true) result that a force 
which attracts locally , rep els globally. This qualitative view is supported by the work of 
Maor fc Lahavl (120051 ) and iMaorl (120071 ). where they found that for a two-component system 
made up of matter and dark energy, the virialization radius will increase when the Dark 
Energy is allowed to dynamically participate in the virialization, a sign of an additional 
attractive force (due to the increased virial kinetic energies occurring in a deeper poten- 
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tial well). One implication of this is that signs o f a cosmologically recent slow-down in the 



growth of clustering (e.g., IVikhlinin et al.l 120091 ). findings which are viewed as strong evi- 



dence in favor of Dark Energy, really only support a Cosmological Constant form of DE, not 
any more general DDE material with locally- attractive negative pressure. It therefore seems 
that a negative pressure substance like dynamical Dark Energy would be exactly the wrong 
prescription when looking for a cosmic component that can remain smoothly distributed in 
order to avoid enhancing the small-scale clustering of material - thus indicating an essen- 
tial conflict between the dual requirements of smoothness and acceleration for just about 
anything other than the immobile, coincidence-prone A. 

This potential instability of DDE to spatial p erturbations is well known by experts in the 



field (e.g., [Turner fc Whitelll9971 : lHulll998 



20051 ); and various forms of support pressure may 



be proposed in order to subvert the natural behavior of negative pressure substances towards 
small-scale clumping, thus allowing a DDE to remain smooth on scales large enough to evade 
limits on flu from observations of large-scale structure formation. But beating the laws of 
thermodynamics is a tricky business to be in, and often one errs by implicitly introducing a 
source of positive pressure somewhere in the problem - perhaps within the Dark Energy itself, 
such as by considerin g relativistic particles (with their ordinary, thermodynamic positive 
pressure) as the DDE ([Turner &: Whitdll997l ); or by using the (positive) de generacy pressur e 



of long- wavelength particles like neutrinos to keep them unclumped (e.g., I Wigmansl 120041 ) . 



Or alternatively, positive pressure can be generated by an auxiliary material invoked i n 



combination with the DDE to keep the latter from clustering (e.g., iBjaelde et al.l 120081 ). 
Either way, a positive pressure source sufficient for supporting the DDE against collapse must 
invariably contribute a decelerating contribution to the cosmic expansion (cf. Equation [I]) 
that counteracts the accelerative properties for which the Dark Energy was recruited in the 
first place. 

Even without any implicit, possibly inadvertent introduction of a positive (adiabatic) 
pressure source, it could still be argued that some form of nonadiabatic pressure within 
the DDE might succeed at preventing its small-scale clust ering without ru i ning the actual 



cosm ic acceleration. For a scalar field, it is well known (iGrishchukl Il994t ICaldwell et al. 



19981 ) that the sound speed is scale-dependent; and despite being imaginary (thus unstable 
against collapse) with c 2 s < for long-wavelength perturbations of a w < scalar field, for 
small wavelengths the sound speed instead approaches c 2 s ~ c 2 , thus seemingly solving the 
instability problem automatically when considering small-scale clustering. But the situation 
is not quite so simple, however, when considering the detailed processes of inhomogeneous 
structure formation. For example, in a study modeling the spatial e yolution of scal a r field 
DDE in the presence of a spherically-collapsing matter overdensity, iDutta fc Maori (120071 ) 
find that the initial response of the DDE is to join the matter in entering into a collapsing 
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phase, thus adding to the growing central overdensity. This collapsing phase does not last 
long, since the slower Hubble expansion in the central matter overdensity allows the scalar 
field to locally roll down its potential and lose energy, actually creating a weak (Sp^/p^ ~ 
—0.02) void of DDE there, thus appearing once again to solve the small-scale DDE clustering 
problem. But recalling the 1 st Law of Thermodynamics (P = —dE/dV), a material that 
loses energy when the universe expands is m ore of a positive pres sure substance than a 
negative pressure substance; and sure enough, iDutta fc Maori (120071 ) find that the effective 
equation of state w of their DDE does become substantially less negative as the system 
evolves, particularly in the center where the DDE void is developing. Thus we see that 
the less a DDE engages in small-scale clustering, the less it behaves like a negative-pressure 
substance capable of accelerating the universe; and it is not unreasonable to suppose that 
these two properties are fundamentally connected in a thermodynamically essential way. 
Now, we should not overstate the case that can be made from those results, since the scalar 



fields examined in IDutta fc Maori (120071 ) never actually lost their accelerative nature (i.e., 
w > — 1/3 is never reached): the equation of state only increases from w = —1 (everywhere) 
at z = 35 to w ~ —0.8 (at the center) at z = for their simplest scalar field potential (and 
increases only to w ~ —0.945 for their more theoretically-motivated double exponential 
potential). But then again, their numerical analysis only studied the system within the 
regime of linear perturbations; and it would be interesting to see what happens to the 
averaged equation of state of the total cosmic contents (matter plus DDE) - in order to 
evaluate the actual ability of such a cosmic mixture to effectively accelerate the universe 
- when a significant fraction of the DD E is located withi n regions of highly nonlinearly 
collapsing matter overdensities. (While iMota et al.l (120081 ) do perform a analysis that is 
fully nonlinear for the matter perturbations, they do not evaluate any effective changes 
to w which would occur as a result. Furthermore, the assumption made in both of these 
papers of a very small mass scale for the DDE scalar field, ~ O(H ), may inadvertently be 
recruiting degeneracy or relativistic pressure within the DDE its elf to help ke e p it u nclumped 
on small scales; this might explain the seemingly odd result in IMota et al.l (120081 ) that the 
DDE cannot clump onto pre-existing overdensities by itself, but instead stops increasing its 
density contrast once the nonlinear matter overdensity pulling the DDE in stabilizes its own 
collapse through virialization of the matter alone.) 

And even beyond the question of changes to the effective equation of state, w, there is 
a more subtle issue involved here as well, which we pose as the following question: if the 
clustering properties of a substance with nonadiabatic pressure do not solely depend upon w, 
then why would its accelera tive properties depend solely upon vol In the textbook derivation 



Kolb fc Turnerl Il990l ) of the acceleration equation, Equation [TJ one proceeds by first 



(e.g. 

assuming a perfect-fluid form for all cosmic constituents; but a DDE with its inhomogeneities 
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controlled by nonadiabatic pressure would clearly fail to behave as a perfect fluid, placing the 
very simple form of the FLRW acceleration equation itself into question (as well as muddying 
the physical meaning of the equation of state parameter, w = P/p). If one cannot trust the 
DDE to behave as a perfect fluid in one sense, then there is less reason to trust it to behave 
as a perfect fluid in any other sense; and thus it seems quite reasonable to believe that a 
dynamical substance sufficiently exotic as to avoid clumping despite P < 0, might also fail 
to generate a > despite P < —p/3. Rather than simply quoting a sufficiently negative 
value of w as demonstrating a claim that some non-clustering DDE is capable of accelerating 
the universe, a fully general-relativistic treatment of all of the inhomogeneous-DDE-induced 
accelerative and/or decelerative effects within a fair cosmological sample (clusters and voids) 
would need to be done to properly substantiate that claim. 

While the foregoing discussion does not constitute actual proof against the possibility 
of achieving acceleration with a negative pressure, mobile substance that somehow fails to 
clump, it should at least introduce a degree of skepticism into an idea that has generally 
been accepted at face value. And if these concerns about coincidences and negative-pressure 
instabilities do in fact cause one to abandon both A and DDE, then explaining the cosmic 
acceleration would seem to present a conundrum. Fortunately, however, the self-attractive 
nature of negative pressure presents us with the seed of an alternative idea: in its own way, 
normal gravitational attraction represents a form of negative pressure - specifically, the very 
same gravitational attraction which is involved in the growth of galaxies and galaxy clusters 
in large-scale structure formation. One might therefore make a virtue of necessity, and try 
to recruit structure formation itself as the driver (not just the trigger, as in some tracker 
quintessence models) of the cosmic acceleration. Such a step would completely solve the 
Coincidence Problem, since the remarkably contemporary onsets of both the acceleration 
and the existence of galaxies (and hence planets, life, and astronomical observers) could be 
viewed as nothing more coincidental than finding two neighboring apples that have fallen 
from the same tree. 

It was the realization of the paradoxically accelerative yet locally attractive nature of 
P < for Dark Energy that initially led this autho r to search for a structure- formation- 



induced solution to the cosmic acceleration problem (jBochnerl l2007al ) . During the past few 



years, a number of other researchers have also propos ed a variety of methods (e.g., ISchwarz 



20021 ; lRasanenll2004j ; iKolb et al.ll2005al ; IWiltshirdl2007l . as a few examples) of trying to obtain 



the observed cosmic acceleration from the breaking of FLRW homogeneity and the formation 
of structure, without requiring the intervention of non-standard gravitational or particle 
physics, or the non-Copernican approach of placing us near the center of a large cosmic 
void. This type of effort has generally become known as "backreaction" . The backreaction 
paradigm, in the consensus view, has so far been unable to completely replace Dark Energy 
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as the source of the apparent acceleration (we will discuss some of the reasons for this 
belief in the ne xt section); and the arguments which this author originally presented in 



Bochnerl (l2007al ) were advanced without the benefit of a plausible mechanism for realistically 



generating an acceleration. 

In this paper, however, we will introduce a physically plausible mechanism for produc- 
ing an observed acceleration; and though our method here is strictly phenomenological, it 
serves as the basis for a formalism that we can use to not only reproduce the cosmic accel- 
eration (as indicated by observations of standard candles such as Type la supernovae), but 
also, simultaneously, to reproduce several of the other key features of the apparent ACDM 
concordance. 

This paper will be organized as follows: in Section [21 we discuss several of the oft-quoted 
limitations of the backreaction mechanism, countering them, and then use such arguments 
as guideposts towards developing our formalism in Section [3j in Section HJ we describe our 
specific models and present our simulated Hubble curves; in Section [5j we show how to 
extend these results to the formulation of a new "Cosmic Concordance" ; in Section El we 
evaluate the success of this alternative concordance, and then suggest potential methods for 
observationally testing our paradigm in order to distinguish it from ACDM (and perhaps 
also from DDE); in Section [71 we discuss the possible long-term futures ('fates') of a universe 
dominated by our "causal backreaction" mechanism; and we conclude with a summary of 
these ideas and results in Section [HI 



VARIETIES OF "BACKREACTION" AND THEIR COSMOLOGICAL 

IMPACT 



There are several distinct ways of achieving an observation that looks like a cosmic ac- 
celeration, and one ma y categorize the dif f erent approaches or sce narios in various different 
ways (see, for example, iBiswas et al.ll2007l ; iBiswas fc Notaril 120081 ) . Here we will find it use- 
ful to characteriz e the various forms of o bserved acceleration according to the (qi, q 2 , qz, q^) 
terminology from lHirata fc Seljakl ( 120051 ) - where q\ refers to an 'actual' acceleration as de- 
scribed by an increase in the volume expansion 9; q 2 refers to an 'averaged' acceleration for 
inhomogeneous spacetimes, an effect strictly due to nonlinearities in general relativity (GR) 
when one has departures from N ewtonian metric perturbations (i.e., the "fitting problem", 
Ellislll984j ; lEUis fe Stoegerlll9871 ); and q%, q^ refer to an apparent acceleration, due to the lim- 
itations that one always experiences in the act of observation, attempting to infer the general 
cosmic evolutionary behavior from a circumscribed set of luminosity distance measurements. 
(The difference between g 3 and q^ has to do with how they perform angle- averaging for an 



-8 - 



anisotropic universe - for example, regarding the possibility of finding 'acceleration' in some 
directions, and 'deceleration' in others.) 

In general, q 3 and may be strongly dependent upon the observer's position in space, in 
relation to inhomogeneities, voids, etc., as well as being affected by anisotropic observations. 
For our formalism, however, which we will be describing as "smoothly inhomogeneous" , 
the observed acceleration will depend neither upon the position of the observer nor the 
direction in which one looks, but will instead depend solely upon the time of the observation. 
Technically speaking, only q 2 - which relies directly upon the nonlinear character of GR - 
was originally termed "backreaction" , although common usage has generally broadened the 
term to mean any form of observed or apparent acceleration from cosmic inhomogeneity. In 
this paper, we will use the term in its broadest sense, to mean any or all of (qi, q 2 ,q3, q±)- In 
the full complexity of GR for an inhomogeneous universe, it is not always a simple matter 
to disentangle these different terms in a clear-cut fashion; but our formalism presented 
here is likely best described as being the ultimate result of an averaged gi-type volume 
creation, perhaps in combination with (fe/ffe-type effects due to the fact that gravitational 
perturbation information coming in from far away - and thus from earlier look-back times, 
when the universe was less clustered - actually allows the late-time/small-distance cosmic 
expansionary behavior to be appear quite different from what we see from far away, but only 
in a purely observational sense (i.e., not due to any real local differences), and in a way that 
any observer would see for their own local cosmic neighborhood. 

At this present moment in the general theoretical consensus, backreacti on is not widel y 



considered to be a viable method of achieving the observed acceleration (e.g.. lSchwaral2010l ). 
There are a number of important reasons for this, which seem on the surface to be strong 
arguments; but we will explain here why despite being technically correct, several of those 
arguments do not properly apply to a physically realistic cosmological model. 



2.1. Smoothly-Inhomogeneous Backreaction: No Voids, Holes, or Special 

Boundary Conditions 



One simple way of achieving an apparent acceleration via FR W- violating inhomogeneities , 



is to imagine our selves located within an underdense void (e.g., iTomital l200ll ; iMoffatl 12005 
Alnes et al.l 120061 ) . Such a void, expanding at a faster rate than the average cosmic ex- 
pansion overall, would imprint Hubble curves with what looks like the signs of a recent 
cosmic acceleration, as incoming light signals from standard candles cross over from slower- 
expanding regions farther away from us, to the faster-expanding regions inside our local 
void. Such a void solution has difficulties, however - including a coincidence problem of 
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its own, regarding our fairly suspicious (if random) position very close to the center of the 
void, as would be needed to maximize 'acceleration' effects while minimizing the resulting 
anisotropies imposed upon the observed CMB radiation; as well as the necessary existence 
of a substant ially large and deep vo id that would be capable of achieving sufficient apparent 
acceleration ( lAlexander et al.l 120091 ). The existence of such a void woul d eventually becom e 
detectable using luminosity distance data over intermediate redshifts (IClifton et al.l 120081 ); 
and though it has been argued that a variety of cosmo logical observations remained com - 
patible with the possibility of such a large local void (jGarcia-Bellido fc Haugbolld 120081 ). 
recent Type la supernova data is co nsistent with a determination that no so-called "Hubble 
bubble" exists ( IHicken et al.ll2009bl ). a rejection of the void hypothesis furthe r substantiated 
by new Cepheid data used for calibrating such supernovae (IRiess et al.l 120 111 ). 



Perhaps fortunately, therefore, our formalism does not require the existence of a special 
local void centered upon us, but instead uses the overall, averaged effects of many inhomo- 
geneities, all summed together, in order to achieve the desired observational result. But the 
question necessarily arises as to how to represent a situation with multiple inhomogeneities 
via an appropriate cosmological model. One answer to this difficulty that has been studied 
in detail is a class of models known as "Swiss-Cheese" solutions, in which one starts with a 
homogeneous FRW background (the "cheese"), and carves out spherical regions (the "holes") 
in which matter can be distributed inhomogeneously - though still spherically symmetrically 
- in such a way as to represent many independent clumped structures of matter. One can 
use any spherically symmetric (dust-only) metric that one wishes for modeling the holes, 
such as the stati c (fully collapsed) Schw arzschild metric, as used in the original Einstein 
Straus solution (lEinstein fc Straus! 119451 ); or the radially-evolving LTB metrics (ILemaitre 



1933 



Tolmanl I1934J : iBondil 119471 ). as used in various recent studies (e.g., iBiswas fc Notari 
20081 ). The ke y here is that b y retaining spherical symmetry - thus recruiting Birkhoff's 
Theorem (e.g., IWeinbergl 119721 ) - and by carefully matching the mass-energy content of the 



hole to the exterior cheese at the hole boundary, one aims to completely zero out the effects 
of the existence of the hole on its surroundings. Thus the exterior regions continue to expand 
exactly as they would in an entirely homogeneous FRW cosmology, as if no holes existed, 
while giving one the ability to produce exact solutions for detailed study which do include 
deviations from homogeneity. Since each hole is presumed to have no effect on anything 
exterior to itself, one can include any number or variety of holes (as long as they do not 
intersect one another), even to the extent of filling the universe entirely with holes (and 
thus wit h essentially no che ese anywhere) in a fractal-like pattern, such as the "Apollonian 
gasket" (IMarra et al.l 120071 ) . Such a model would have the same average expansion as a 
FRW-everywhere model, despite actually being FRW virtually nowhere. 



Despite being mathematically convenient for study, however, such models have serious 
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disadvantages. One key shortcoming is that since they obey FRW dynamics on average, 
and since the holes have no effect on the exterior cheese, there is by definition zero actual 
"backreaction" : all effects causing apparent deviations from matter-dominated FRW expan- 
sion are purely observational, due to effects upon the photons passing through the holes. 
Ob viously, this will inhibit t he ability of the inhomogeneities to mimic a cosmic acceleration; 
and lBiswas fc Notaril (120081 ) find that such effects are indeed strongly suppressed, by a factor 
(L/R H ) 3 for passage through a hole, where L is the size of the hole and Rh is the Hubble 
radius. Even if one integrates the paths of light rays through many holes, this is clearly too 
small of an effect to replace the Dark Energy. 

Such suppression of backreaction effects is entirely artificial, however, since the Swiss- 
Cheese models themselves are artificial - they are, in fact, entirely inadequate for plausible 
studies of the real universe, since they lack so many key physical features. First, the spherical 
symmetry of the holes makes them unabl e to incorporate the viri alization and stabilization 
of forming structures (a point noted by iBiswas fc Notaril (120081 ). and which will take on 
more importance below). Second, the prescription of taking all effects of the hole back 
upon the exterior universe and setting them to zero is entirely unphysical - when condensed 
structures form in the real universe, rather than leaving exterior regions alone, they instead 
tend to become overdense attractors which keep growing in mass and perturbative cosmo- 
logical influence; co nsider for example our acceleration towards the Great Attract or/Shapley 
Concentration (e.g.. lLucey et al.ll20 05; Bo leiko fe Hellabvll2008t IColm et al.ll2010l ). and other 
evidence of large-scale bulk flows ( Feldman et al.l201ol ). Third, the notion of disjoint, non- 
intersecting holes - a necessary assumption for producing exact metric solutions, requiring 
special and implausible boundary conditions at the hole/cheese interfaces - is a fatal failure 
of the model, since in the real universe, the spheres of influence of individual inhomogeneities 
always continue to expand as they pull in more material (and have more time to spread their 
causal influences), until the separate perturbations eventually overlap and merge. The super- 
position of gravitational perturbations from many different inhomogeneities is a phenomenon 
that Swiss-Cheese models cannot represent at all. Simply speaking, if we recall the earlier 
suppressive factor of (L/Rh) 3 - and now consider that each hole (overlapping with all of the 
other holes) is about the size of the observational horizon out to which its gravitational pull 
can be causally felt - then this factor is much closer to unity, and not much of a 'suppression' 
after all. 



The challenge, of course, is how to represent such a model mathematically; since once 
the holes all merge together, it seems like we are back to the homogeneous FRW model 
again, in some averaged sense. The metric that one uses, in any case, cannot depend upon 
spatial position, since we are only considering averaged effects. But what we will justify 
below, in this section - and present a phenomenological model for, in Section [3]- is that the 
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perturbations can indeed be perpetually nontrivial functions of time (and of time alone); 
and the fact that we live merely in an averaged FRW universe (instead of in a truly homo- 
geneous FRW universe) can be modeled as alterations to the behavior of the average scale 
factor, a(t), in such a way that it looks like we live in a FLRW universe with matter and 
Dark Energy, rather than dust alone. This model, in which we re-establish the mathematical 
homogeneity and isotropy of the Cosmological Principle, while retaining the average, inte- 
grated effects of many overlapping inhomogeneities (incorporated as a perturbation to the 
Friedmann expansion), is what we refer to as a "smoothly-inhomogeneous" universe. As will 
be shown in Section HI these effects are capable (under realistic conditions) of reproducing 
the entire observed acceleration; and an important benefit of the smoothly-inhomogeneous 
model is that there is nothing special about our own position, and so all observers (barring 
strong local influences) would see the same apparent acceleration at the same cosmic time. 



2.2. Vorticity as a Large-Scale Player 

During the evolution from 'truly homogeneous' to 'smoothly inhomogeneous', a strong 
phase transition takes place throughout the entire universe. The question is whether or not 
the effects of this phase transition would have large 'accelerative' effects - or as some have 
argued, whether or not it can have any accelerative effects. Sharp limits have been placed 
upon such mechanisms by arguments in the literature, and so we must explain here why 
those limits are not applicable to realistic situations. 



It has been argued, for example ( lAlnes et al .1120071 1. that cosmologies dominated by 



pressureless dust do not permit acceleration (i.e., true volumetric acceleration, in the sense 
of qx). This conclusion was based upon a study of LTB inhomogeneity, which as we have 
noted is limited to spherical symmetry, and lacks the ability to model virialization; but, 
w hat is the true importa nce of virialization? Consider the q\ acceleration "no-go" theorem 



of Hi rata fc Seljakl (120051 ) . that was based upon the Raychaudhuri equation ( [Hawking Ellis 



19731 ). which assuming perfect fluids can be written as: 

df) f) 2 

j t = -- - - 4vrG(p + 3P) + , (2) 

where 6 is the volume expansion, a^ u is the shear tensor, u^ u is the vorticity tensor, and 
p and P are, respectively, the density and the isotropic pressure of all material contained 
within the stress-energy tensor T^ v . From this expression, it is obvious that if the vorticity 
is zero (or more precisely, if u 2 = u^u^ jl = 0), and if the strong energy condition (SEC) 
is obeyed (p + 3P > 0) for all species, then clearly d6/dt < 0, and no acceleration is possible. 
On the other hand, any violation of the SEC due to negative pressure such that p + 3P < 
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(i.e., w < —1/3), strong enough to lead to d6/dt > 0, would necessarily be classed as a form 
of Quintessence, Dark (or Phantom) Energy. Therefore, in a matter-dominated universe 
with negligible pressure (P ~ 0), with this dust being largely irrotational (u 2 ~ 0), all of the 
remaining nonzero terms would provide only negative (decelerative) contributions to ensure 
that dO/dt < 0. 

Here, however, is one of the key questions of this paper: Why is the vorticity set to zero ? 
This seems to be an especially unusual prescription in a universe where nearly everything is 
rotating and/or revolving about something. After all, all virialized structures in the universe 
larger than individual solid objects^ are stabilized against collapse by some version of vortic- 
ity or velocity dispersion. Henceforth referring to all varieties of such mechanisms as simply 
"vorticity" , this includes both the organized rotational motions of solar system planets and 
spiral galaxies, as well as the less organized orbital motions of stars in elliptical galaxies, 
and of individu al galaxies in galaxy clusters. The virialization mechanism itself (e.g., vio- 



lent relaxation; iLynden-Belll Il967t IShul Il978[ ) achieves such stabilizations by generating or 



concentrating all of this vorticity. The value of u 2 within galaxies and galaxy clusters must 
obviously be large - and in fact, physically dominant - since it is clearly sufficient to counter 
the gravitational attraction of the matter (including Dark Matter) to produce static regions 
of space with d6/dt = 0. 

It is certainly a bad physical approximation to take the dominant physical force (in op- 
position to gravity) that regulates cosmic structure formation^, and set it equal to zero. We 
must therefore consider the reasons for which this is typically done. Mathem atically, zero vor- 



ticity is necessary for the existence of the synchronous comoving gauge (e.g.. lHirata fc Seljak 



20051 ). such as is used for the homogeneous and isotropic FLRW cosmologies; and vortic- 



ity causes pro blems with the meshing together of the hypersurfaces of proper time (e.g., 



Rasanenl 120061 ). Thus one cannot define an idealized "Hubble flow" unless one sets the ro- 
tation to vanish; and so to use this common and useful simplification, it must be assumed 
that on large enough scales, one is justified in averaging any cosmological fluid to an irrota- 
tional, comoving state. In essence, therefore, incorporating vorticity explicitly would cause 
extreme complications for the entire mathematical machinery of how cosmological evolution 



1 Solid objects are supported by their internal pressure gradients and body forces, which also contribute 
(though likely on a cosmologically small level) to dO/dt > 0, but which have been dropped even before getting 
to Equation [2] as a result of the perfect fluid approximation. 

2 All of this vorticity requires great amounts of motion - i.e., kinetic energy. This energy comes from the 
gravitational potential energy of inhomogeneities originally pulled apart by the expanding universe. Thus the 
matter is 'pumped' by the expansion; and this accumulation of energy taken from the expanding universe is 
how the gravitationally-self-attractive cosmic matter effectively behaves like a 'negative pressure substance'. 
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and structure formation are usually studied in practical Such complications, barring the 
unlikely possibility of rotation on supercluster or even larger scales, would appear to be 
physically unnecessary to keep track of on a truly cosmological level. Thus vorticity has i n 
many cases been relegated to lesser status as a "small-scale player" (e.g., iBuchertl 120081 ). 



relevant only for cosmic averages performed over domains that are on or below the scales of 
galaxy clusters. 

The problem with this approach is that it is not actually the vorticity tensor, u^ u , which 
appears in the Raychaudhuri equation as a positive contributor to dO/dt, but the square of 
the vorticity, u 2 . And no matter what scale one goes up to in their averaging, the square of a 
quantity will never average away; in fact, the total amount of integrated u 2 in a volume will 
be strictly proportional to volu me, if the clustering and s tabilization of mass is essentially the 
same everywhere in space. As iBuchert &: Ehlerd (119971 ) put it, cosmic averages of positive 
semi-definite quantities like (u 2 ) (or in opposition to it, (cr 2 )) get "frozen" at the size of 
"typical subdomains" ; and thus they cannot be made to go to zero by averaging over larger 
domains, even if the spatially averaged value of the parent tensor, (u^ u ), itself becomes 
negligible on large scales due to collective cancellation from neighboring local subregions in 
the domain. Furthermore, the evolving cosmic value of (a; 2 ) is not subject to any particular 
conservation law (since uj^ throughout the early universe was not exactly zero), nor is uj 2 
limited from growing as large as it needs to be in any lo cation in order to viria lize a self- 
stabilizing region of clustered mass - formally speaking, lAsada fc Kasail (119991 ) find that 
vorticity is coupled to density contrast, becoming strongly amplified in locally collapsing 
regions, contrary to the expectations from standard linear perturbation theory. Systems 
demonstrating significant vorticity-dependent behavior should be quite natural in the real 
universe, as depicted generically in Figure [1] (with angular momentum and its squared value, 
L and L 2 respectively, serving as a proxy for u^, u 2 here). We can even observe the results 
of this kind of behavior in our nearby cosmic neighborhood by looking at the Local Group, in 
which the Milky Way and Andromeda galaxies each support themselves against collapse via 
their own individual spins; and yet Andromeda and the Milky W ay, like many pairs of large 
spiral galaxies, are counter- rotating with respect to one other (ISchaail Il975l . p. 174), thus 
partially canceling (u^ u ) for the Local Group as a whole (though not necessarily implying 



3 Thcrc are of cour se a few analyses which d o not immediatel y drop c osmoloKical vorticity at the outset 
of calculations : e.g., Ehlers fc Buchert ( 1997 ). Asada fc Kasai ( 19991 ). Christopherson et al. ( 20091) . and 



Rasanenl (|2010l ). as a few examples. But this does not guara ntee that vorticit y is being included to the degree 



necessary for generating causal backreaction. For example, Rasanen ( 2010l ) finds the relevant contribution 
due to vorticity to be simply a total divergence, leading only to a negligible surface term; but they restrict 
their analysis to first order in velocity v - and as will be made clear below in Section [2.31 an analysis at least 
up to 0(v 2 ) is necessary in order to incorporate the causal flow of gravitational information. 
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that it is small or conserved per se - e.g., iDunn &: Laflammd ( 119931 )). without reducing its 
total integrated value of (uj 2 ). 

Given the potential significance of vorticity in regards to both astrophysics and vol- 
ume exp a nsion, general warnings against neglecting it are given (multiple times in fact) by 



Rasanenl (120061 ) (and by various other authors), in the context of acceleration via struc- 
ture formation; but aside from such caveats, for typical studies in the literature the issue 
of vorticity is nearly always dropped before actual calculations are done, and models which 
completely lack the ability to represent vorticity or virialization (such as LTB) are usually 
used for the estimates of backreaction-related effects of inhomogeneities on the observed 
cosmic expansion. 



As iHirata fc Seljakl (120051 ) note, the absence of vorticity plays a key role in all of the 
"no-go" theorems used to prove the impossibility of accelerated expansion in the absence of 
some kind of SEC-violating Dark Energy. And therefore, due to the crucial importance of 
vorticity (and thus a substantially nonzero u 2 ) in structure formation for all of the virialized 
inhomogeneities pervading the universe, we can conclude that all acceleration "no-go" the- 
orems are cosmologically inapplicable. A smoothly-inhomogeneous cosmology with no Dark 
Energy, but with the overall effects of pervasive virialization, can indeed (as we will see) be 
able to produce enough of a backreaction to create (at the very least) an apparent acceler- 
ation, if not a full-fledged volumetric ('real') acceleration. Vorticity and velocity dispersion 
effects are not in fact a small-scale player in the physical universe, but represent the ma- 
jor player (along with gravity), and the only question that remains is to figure out how to 
properly incorporate them. 

The difficulties involved in explicitly including vorticity within cosmological models are 
very real, however, and it would be a highly nontrivial problem to try to precisely compute the 
detailed time evolution of (and volume expansion from) a realistically virializing structure. 
Fortunately though, as we contend here, there is no need to solve this exact problem in 
order to get a good approximation of the averaged effects of backreaction upon the observed 
expansion rate. Instead, we can exploit the fact that both the beginning state (nearly 
perfectly smooth FRW matter) and the ending state (a reasonably random distribution of 
discrete Newtonian mass concentrations) are extremely simple, making the net results of 
the phase transition to inhomogeneity very straightforward to estimate. Only the general 
evolution of how much of the cosmic matter has completed the phase transition to dumpiness 
as a function of cosmic time is important, not the detailed time dependence of the virialization 
of any individual object; and for this paper we will use a set of reasonable heuristic models 
to empirically determine which ones are best suited for fitting the observational data. This 
phenomenological formalism will be motivated and presented below, in Section [3j 
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2.3. Newtonian-Level Backreaction: Not Suppressed, Not Small, and Not 

Slow 



There remain debatable issues regarding the possible magnitude of such backreaction 
effects, however, given the fact that we continue to rely here upon (individually) Newtonian 
perturbations, and do not explicitly recruit the effects of nonlinear-strength gravitational 
fields. Powerful theoretical arguments have been advanced by previous researchers, claiming 
that cosmological backreaction from Newtonian-level perturbations must necessarily be very 
small, if not actually suppressed to zero, in a practical sense. Given that the development 
of Newtonian perturbation terms will provide positive contributions to the spatial metric 
components (to g rr , in particular, for the spherically-symmetric case), it is unclear why new 
volume expansion is not considered an inevitable result of those newly-developing Newtonian 
perturbations; but it is this assumed restriction for backreaction, apparently requiring one 
to appeal to higher-order gravitational terms, which is what has forced researchers to resort 
to desperate measures in order to achiey e a significant backreaction. This author's original 
attempt at backreaction (lBochnerll2007al ) . for example, postulated a toy model utilizing black 
holes, in order to generate a strong enough effect to make a difference in the cosmic expansion 
(a clearly unworkable model for several reasons, but no alternative was obvious at the time). 
Using another exotic approach, other researchers have attempted to recruit super-Hubble- 
sized density perturbati ons - first on their own (e.g. , iKolb et alKiiiiM , and then in tandem 
with sub-Hubble terms ( IKolb et al.ll2006l ) - in order to create a usefull y large inhomogeneity- 
induced backreaction, in spite of the causality-related difficulties ( jrlirata fc Seljak 2005 ) 
of trying to produce p hysically measurab le effects from perturbation modes larger than 
the observable horizon. IKolb et al.l (120081 ) find it necessary to go to the extent of arguing 
that a perturbed conformal Newtonian metric is not even an appropriate description for a 
universe observationally well-described by a ACDM model, when considering deviations from 
an Einstein-de Sitter dust model as the unperturbed background. It is clear, therefore, that 
the permissibility of a Newtonian description of inhomogeneous perturbations is commonly 
viewed as an insurmountable impediment for backreaction. We must therefore examine the 
reasons for this belief, and show how it comes about from an unwarranted oversimplification 
of cosmological physics. 



From the exten sive work of Buchert and collaborators (e.g., iBuchert &: Ehlersl 11997 



Buchert et al.ll2000l ). it is well known that the entire Newtonian backreaction, Qn, can be 
expressed mathematically as a total diverg ence: Qn = Qd%v = V-qrjy. We know from 
the divergence theorem (e.g., iJacksonl 119751 ) that for any such function integrated within a 
volume V (with boundary 5V), we can relate the volume- integral of Qn to the boundary- 
normal surface integral of qw- J v Qn = §gvQ N ' n - This surface integral will vanish if 
the universe is topologically closed (i.e., if there is no surface), or if periodic boundary 
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conditi ons are assum e d (as is typical for large-scale Newtonian cosmological simulations 



e.g., 



Springel et al 



( 120051 )). meaning that the total integrated Newtonian backreaction 
(J v Qn) would c onsequently have to vanish. But even if it does not completely vanish, 
Kolb et al.l ( 120061 ) argue that the overall effect will in any case be tiny, because the surface 
term involves the peculiar velocity u, which should be small in Newtonian situations; and 
perhaps another limiting factor here is that an integral over a surface of radius R will only 
increase as ~i? 2 , so that the volume average will go like ~i? 2 /i? 3 ocl/R, becoming negligibly 
small for arbitrarily large R. Thus the ability to represent Newtonian (i.e., gravitationally 
linear) backreaction as a total divergence, if valid, would act as a strong suppression of the 
magnitude of the effect. 

It is interesting, at this point, to consider the basic equations at the heart of the for- 
malism used for pl acing these limits upon the Newtonian-level backreaction - specifically, 
Equations la-Id in iBuchert fc Ehlerj ( 119971 ) - and to compare them with Maxwell's Equa- 
tions (along w ith the Lorentz force law and continuity equation) for electromagnetism (e.g., 
Jackson! Il975l ). A careful examination of these expressions shows that for the former, some- 
thing is missing: the magnetic fields. The Buchert formalism, and the various theorems 
derived fro m it, have all been d one having set what one may call the "gravitomagnetic' : 



fields (e.g., iMashoon et al.ll200fl ) to zero. Now, the neglect of the 'magnetic' fields in elec- 



trodynamics (or in Newtonian gravity) happens to be a very significant simplification, as is 
obvious from a quick examination of Ampere's Law, Vx£?= (4Tr/c)J+(l/c)dE/dt. If B is 
set to zero, then the time rate of change of E in a region is solely dependent upon the local 
current /momentum flow, J, inside that region. In other words (say, for gravitation), the 
contributions to the local Newtonian gravitational potential and fields by distant masses out- 
side the local volume are completely frozen-in, with any possible causality-driven dynamics 
being entirely suspended. This semi-static formulation (i. e., considering only the motions of 
local masses) is entirely in accord with the common view (IBuchert et al.ll2000l ) that: "In the 
Newtonian approximation the expansion of a domain is influenced by the inhomogeneities 
inside the domain." The main thrust of our arguments in this paper is to point out that this 
view is unacceptably incomplete, and in fact we rejected the Swiss-Cheese approach earlier 
on this very same basis, that it is physically unrealistic to bar different regions of spacetime 
from communicating with one another. While such an approximation may be reasonable ac- 
cording to what researchers usually think of as a 'Newtonian' approximation, it is certainly 
not compatible with the dictates of special relativity and cosmological causality. 

For an al ternate way of exp r essing Buchert-like restrictions on backreaction, consider the 
arguments of llshibashi fc Waldl (120061 ). in regards to the following, "Newtonianly- Perturbed 
FLRW" metric: 

ds 2 = -(1 + 2§)dt 2 + a 2 (t)(l - 2<S>)-f ij dx i dx j . (3) 
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In analyzing this metric, they adopt the various 'Newtonian' conditions for the potential 
function, $: 

|$| < 1 , (4a) 

\d$/dt\ 2 < — D%D^ , (4b) 



a 



(D^Di®) 2 < {D i D j ^)D i D j ^ . (4c) 

They point out that even in the regions of very large local density variations, 5p/p, the 
metric perturbation $ due to these spatially-varying matter concentrations should always 
remain small everywhere (with the exception of extreme regions, such as near black holes); 
and that given these conditions, the nonlinear co rrections to the expansio n will remain small. 



The linear effects are small by hypothesis, and llshibashi fc Waldl (120061 ) show that the use 
of Equations Hh-c implies that the remaining dominant backreaction term for $ satisfies an 
expanding-universe version of the Poisson equation, basically reducing it to a total divergence 
and thus subject to Buchert suppression. 

Now , one may concede that spatial variations in $ are indeed small as per the argu- 



ments of llshibashi Sz Waldl (120061 ); and if the temporal growth in $ were also small due to 
Equation I4b| then the perturbation potential $ (being negligible in the early universe) would 
thus have remained very small virtually everywhere and at all times in cosmic history to this 
point, so that little backreaction could ever have been generated from it. But even dropping 
the idea of large spatial variations - or of any spatial variations, which we do completely 
neglect in the case of a smoothly-inhomogeneous universe - why should one assume that 
Equation [4b] is a valid approximation? 



In computing the evolution of structure formation, one typically utilizes the Poisson 
equation, V 2 $ = AirGp (for now ignoring modifications to it needed for representing the 
cosmic expansion). It would of course be okay to follow this approach (as in the Buchert for- 
malism) if it were safe to ignore the gravitomagnetic fields, since the "gravitoelectric" fields 
thus become curl-free. But this is not the rigorously correct formula to use, because the 
Poisson equation is not causal - its (unphysical) solution is the integral of the instantaneous 
Coulomb potential (at that point) due to the charge/mass density distributed throughout 
all of space, integrated out to infinity. While such an approach may be appropriate for con- 
sidering a single, localized physical system embedded in an asymptotically empty universe, 
this picture makes no sense for a system embedded in an essentially infinite universe with 
important nonequilibrium processes going on everywhere, for which the expanding observa- 
tional horizon is always bringing new and important information (and perturbative forces) 
in towards the observation point of that system. 

For a causality-respecting approach, one must instead use the full wave equation (for 
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special-relativistically-consistent potential function <&sr, in analogy with the electric poten- 
tial in the Lorentz gauge, rather than the Coulomb gauge): 

1 <9 2 $« n 

V^ SR --^ = 4nG P . (5) 

Now, the usual impulse is to immediately drop the extra term in Equation [51 involving 
<9 2 $sr/<9^ 2 (equivalent to dropping the gravitomagnetic terms), because of its resulting pref- 
actor of v 2 / c 2 ; this factor would seem to make it very small given the (reasonable) assumption 
of nonrelativistic speeds for most matter flows, and thus (assumedly) ensuring it to negligible 
compared to the spatial variations term in any backreaction calculation. But this thinking is 
based only upon considerations of individual Fourier perturbation modes, not on the overall 
causal behavior of information flow in the struct ure-forming un iverse. If we instead keep all 



terms, and solve Equation [5] as-is, then one gets (I Jackson! Il975l . eq. 6.69): 



d> SR (x,t) = -G r^%^V, (6) 

J l x x I 

where the bracketed numerator is always evaluated at the retarded time, t' = t — |x — x'|/c. 
It is this retarded-time condition which restores causality, allowing different regions of the 
universe to communicate with (and gravitationally perturb) one another; and which provides 
the escape route from the Buchert suppression of Newtonian-level backreaction, because the 
backreaction resulting from this integrated perturbation potential is not expressible as a 
total divergence. We will refer to this propagation of gravitational perturbation information 
between distant (though communicating) regions as "causal updating". 

Now certainly, the factor of v 2 /c 2 multiplying such causal updating contributions is 
typically quite small; but there is a crucial difference between "small" and "suppressed". 
Many individually small contributions can be added together to produce a large overall 
result, if one has enough of them. And the perturbative contributions to the metric at a 
specific location are supplied by every virializing structure in the universe that is within 
the causal horizon of that location. Consider that the gravitational potential (and thus 
the Newtonian-level perturbation effects) of some particular mass at distance r decreases 
only like ~l/r, while the number of such masses within a spherical shell at that distance 
goes like ~r 2 . Summing up the small but nonzero perturbing contributions of all 'remote' 
virializing clumps, in their effects upon any particular observer's location, and integrating 
out to infinity, one gets a total effect oc J°°(r 2 )(l/r) dr = |°°r dr = oo ! Noting once again 
that each individual contribution is from a strictly Newtonian gravitational potential, we 
see that the net effect is not only non-small, but in fact it is infinite (in this very simplified 
picture), rendering any factors of v 2 / c 2 moot. The only effects which actually rein in the total 
integrated perturbation to make it finite in the real universe, are the finite causal horizon out 
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to which an observer can 'see' perturbations (with more distant regions of space obviously 
looking younger and less clumped, due to larger look-back times); and the cosmic (Hubble) 
expansion, which dilutes it by continually carrying the virialized, clumped structures farther 
and farther away from the given observation point. (This sit uation is reminiscent of the 



problem, and the solution, of Olbers' Paradox ( [Weinberg! Il972l ) ; except that the 'infinity' is 
even stronger here, since the Newtonian gravitational potential only fades with distance as 
ocl/r, rather than oc 1/r 2 , as with light intensity.) 

Though this infinity is not physically realized, what the infinite behavior of this simple 
integral actually does alert us to, is that the total perturbation may eventually get very large - 
ultimately defying even the assumed condition |$| <C 1 from Equation Ha]- meaning that the 
entire method of summing Newtonian potentials and using a Newtonianly-perturbed metric 
will break down completely as the backreaction grows strong enough to become general- 
relativistically nonlinear. Thus the summed effect of innumerable, individually-Newtonian 
perturbations is no longer 'Newtonian' in total. (We must therefore watch the calculated 
values of $sr produced in our numerically simulated models very carefully, recognizing the 
loss of accuracy of our formalism when this quantity approaches a value of order unity.) 
But even being finite and (in many cases) Newtonian in total, it is clear that the impact 
of distant shells of clumpy material upon the observer's metric can strongly outweigh the 
effects of inhomogeneities located in the actual vicinity of the observer. 

To the extent that a linearized-gravity approximation is valid, our smoothly-inhomogeneous 
approach therefore works by computing a time dependent, spatially averaged perturbation 
potential $sr(£), which we do not assume is slowly-evolving, as is (inappropriately) assumed 
for perturbation theory treatments via Equation libl Realizing that the summed effect of 
a spherical shell of perturbing clumps at distance r actually increases as ~r, this means 
that the dominant backreaction contributions will be delivered to an observer by masses 
at the farthest possible distances out to which that observer can still see significant matter 
inhomogeneities; while past that point, the retarded-time observations are from a look-back 
time too early in the cosmic history for substantial structure formation to have yet occurred, 
thus making the resulting integral finite. 

Considering the fact that the key metric perturbation function, <3>sR,(t), is predominantly 
affected by perturbation information coming in from distant locations, the retarded-time 
condition of an integrated formula like Equation [6] (suitably modified f or cosmo l oKical cal- 



culations) implicitly gives it the ability (contrary to the argument in iGromesI (120111 )) to 
exhibit relativistic behavior in what would otherwise appear to be nonrelativistic situations. 
To see this, consider the following argument: suppose that there is some approximate tran- 
sition time, tciump, when the universe evolves (not necessarily instantaneously) from being 
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essentially smooth to essentially clumpy. Now at any later time, t > t c iump, the physical 
radius out to where an observer can 'see' this transitional epoch is given (again ignoring 
cosmic expansion here, for simplicity) by r c i ump ~ c(t — t c i um p), a radius which increases over 
time, expanding outward from the observation point (say, Earth) at the speed of light. The 
backreaction effects get bigger and bigger (due to increasing shell size and thus a greater 
number of contributing perturbations) as this outgoing "wave of observed dumpiness" ex- 
pands outward from us at c, thus being fundamentally relativistic in nature and defying the 
assumption of small d$/dt\j. 

Thus we see that the potential function $sr can be fully (special-)relativi stic in charac- 
ter, vi olating; the assumed conditions placed upon 'Newtonian' perturbations by llshibashi fc Wald 
(120061 ). even if none of the matter in the universe is actually moving at relativistic speeds. 
This unexpected effect is due to the fact that cosmological systems (unlike anything else) 
are infinite in size, meaning that the inherently relativistic act of causal observation will end 
up encompassing volumes which are incomparably vast, causing even small effects - which 
never cancel out in gravitation, unlike in electrodynamics - to sum together to produce a 
dominant overall influence. Perhaps surprisingly, therefore, the apparent acceleration of the 
universe that we observe is not caused by any individually powerful or gravitationally excep- 
tional specific objects or regions, but rather from the summed influences of the weak-gravity 
Newtonian tails of innumerable mass concentrations, imposing their combined effects upon 
us from extraordinarily large distances. 

In this approach, it is clear that we are advocating a diametrical ly-opposed view from 
that of other authors who have sought to define a "finite infinity" (e.g.. lCoxll2007l ). represent- 
ing a very limited boundary from within which significant effects upon some specified local 
volume can have ever arrived. Such a limitation woul d restrict the history of 'important' 
interactions to be within an effective "matter horizon" ( lEllis fc Stoegerll2009l ) of only a few 
Mpc in size, delineated by the timelike world lines traveled by pressure-free matter due to 
scalar perturbations. But though they attempt to make the case that effects from outside 
of this matter horizon (yet within the fundamentally c ausa l "part icle horizon" ) are gener- 
ally very small, the chief effects actually estimated by ICoxl (120071 ) include the Weyl terms 
regarding geodesic deviation (oc 1/r 3 ), and the power radiated in terms of the Bondi news 
function (oc 1/r 2 ), both of which are vastly smaller at cosmological distances than the am- 
plitude terms (oc 1/r) perturbing t he m etric itself, which are what we study here. Notably, 
it has been shown by 



Ludvigsenl (119891 ) that even an arbitrarily small energy flux due to 



4 As an analogy, consider the contact point between the two blades of a very long pair of scissors. The rate 
at which this contact point moves outward from the central pivot does not represent the physical motion of 
any real object, and hence is not limited by the speed of the material in the blades as they come together. 
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gravitational waves can result in a finite amount of "geodesic deviation" (actually, long-term 
positional displ acements o f observers due to permanent metric changes) in the very distant 
radiation zone. ICoxl (120071 ) in fact concedes that for galaxies (or clusters, etc.), "the prospect 
of adequately treating such a diffuse body as isolated is doomed", due to such very-long- 
wavelength gravitational radiation (incidentally demonstrating once again the inadequacy of 
the Swiss-Cheese approach). 

The point is that even such small metric alterations can become significant (even dom- 
inant) when summed over the very many contributing sources within one's causal horizon. 
Perh aps t he mo st c oncrete way of demons trating the validity of this idea, is to point out 
that I Cox! (120071 ) and lEUis fc Stoegerl (120091 ) both acknowledge the Great Attractor/Shapley 
Concentration as a likely source of our bulk flow - i.e., as the completely dominant influ- 
ence upon the motion of our entire neighborhood of galaxy clusters - despite the clearly 
contradictory fact (for their claims) that those structures which control our bulk motion are 
exceedingly far beyond any reason able estimate of our so-called matter horizon. Thus the 
assertion in lEilis fc Stoegerl (120091 ) that, "the strengths of any other possible long distance 



influences... gravitational waves, or electric Weyl tensor components from sources outside 
that region - are insignificant compared to local effects," is simply incorrect. 

Despite the incontestable importance of those very distant mass concentrations in deter- 
mining such 'local' motions (and thus necessarily, our local metric), there are very straight- 
forward calculations using gravitational perturbation theory which are commonly believed 
to contradict all assertions of the significance of backreaction in general (and of nonlocally- 
acting causal backreaction in particular). The problem is that despite the standard (almost 
reflexive) tendency for experts to resort to perturbation theory when attempting to evaluate 
the effects of cosmic structure formation, the fact is that perturbation theory is singularly 
ill-equipped for dealing with the most important perturbative effects in an unbounded sys- 
tem like the (effectively) infinite universe. The central difficulty is that perturbation theory 
cannot make any predictions unless one pre-specifies which approximations to make, and 
thus which terms to drop. For example, iKolb et al.l ( I2005a[ ) neglect information-carrying 
tensor modes generated by nonlinear scalar perturbations (i.e., virializing masses) in the 
final expressions for their analysis (as well as assuming irrotational dust, thus neglecting 
vorticity), which leaves one with scalar perturbations only (vector modes are also dropped), 
thus once again limiting the calculation to effects co ming from withi n the inappropriately 
small (as we have seen) matter horizon. Alternatively, iRasanenl (120121 ) does comment upon 
the importance of the propagating degrees of freedom due to the "magnetic" gravitational 
component H a/3 (the same thing as our "gravitomagnetic" fields discussed above), and notes 
the difficulty of studying backreaction using a post-Newtonian scheme without them, given 
that the usual assumption of a finite and isolated system in such calculations is an inap- 
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propriate condition to adopt in a cosmological setting. Nevertheless, iRasanenl ( 120121 ) then 
effectively eliminates such terms by re-affirming the approximation that the time derivatives 
of metric perturbations are small - similar to the approximation above in Equation 14b \ which 
we have labored here to refute on the basis of supporting the causal flow of gravitational 
information. Now obviously, if one removes all "causal" aspects from "causal backreaction" , 
then nothing will be left, and it is unsurprising for one to then find that backreaction fails 
to work, given that the principal physics responsible for it has once again been set to zero. 

The essential problem with perturbation theory in regards to these issues relates to 
its basic program of singling out the 'important' physics by labeling the amplitude of each 
term as "large" or "small", and then dropping the small terms in order to focus upon the 
large ones. But this basic procedure is conceptually flawed in the case of causal backre- 
action, which deals with propagating modes, since such modes can bring in gravitational 
perturbation information to a local observer from vast spatial volumes, so that a term with 
a 'small amplit ude' can actually pro duce an enorm ous overall effect. For example, both 
Rasanenl (120121 ) (as noted above) and iGromesi (120111 ) assume that time derivatives of metric 
perturbations are small, leading to the conclusion that velocity-dependent terms can safely 
be neglected; and IGromesi (120111 ) specifically uses this point as a primary argument for dis- 
missing the effects of causal backreaction entirel}@. But conflating "small" with "negligible" 
is a continuing error in perturbation theory analyses of cosmological evolution, since the 
size of these effects, in terms of amplitude, effectively does not matter: it is practically 
irrelevant how small the time- derivative terms may be, when one realizes that such effects 
are being cumulatively contributed by all of the matter within a causal horizon that may be 
billions of light-years in extent, thus multiplying that amplitude by an amount of mass easily 
large enough to overcome its inherent smallness. Furthermore, even if the amplitudes of the 
relevant perturbative terms for causal backreaction were somehow 'magically' made even 
smaller than they naturally are, this still would not shut off the causal backreaction effect, 
but merely postpone it to begin a little bit later, since a somewhat larger causal horizon of 
self-stabilizing inhomogeneities would then be needed to supply a large enough integrated 
effect to cause an apparent acceleration. In that sense, one could never make the ampli- 
tudes of those perturbation terms small enough to avoid the eventual dominance of causal 
backreaction, since a big enough causal horizon (containing a sufficiently large amount of 



iGromesi (|201 lh also uses an argument referring to their "optimal gauge" , in justifying the neglect of 
causal backreaction; but as they themselves note, it can be a tricky thing to choose the right gauge when 
trying to connect theoretical expressions to real cosmological observables, and a different gauge would lead 
to a different amount of the kind of backreaction effect which we consider here. Thus the permissibility of 
choosing gauge conditions that set components of the backreaction to zero, a procedure which they engage 
in, is a highly nontrivial matter. 
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virializing mass) can always be attained after a sufficiently long period of time, which when 
multiplied by even the smallest amplitude, would eventually produce a term of order unity 
that then proceeds to dominate the cosmic evolution. Thus for causal backreaction in an 
effectively infinite universe, there is no such thing as "negligible" . 



3. DEVELOPING A PHENOMENOLOGICAL FORMALISM 

3.1. Observational Horizons and the Evolving Perturbation Potential 

As long as the universe remains homogeneous and isotropic enough to be considered at 
least smoothly-inhomogeneous, then as pointed out at the end of Section I2.2[ the net result 
of virializing structure formation is quite simple: it consists of the replacement of continuous 
FRW-distributed matter with a discrete collection of stabilized clumps, distributed (one 
assumes) fairly randomly. Considering the metric perturbation effects upon a specific volume 
V, a given stabilizing structure S will only produce nonzero backreaction within V for as 
long as new gravitational information is propagating from the structure to and through that 
volume. Once the contributing gravitational potential of the final, stabilized state of S has 
had the opportunity to completely propagate causally throughout the entirety of V, then as 
per the Buchert limitations on Newtonian backreaction discussed above, the backreaction in 
V due to S will be over. This variety of backreaction, due to a set of individually-Newtonian 
virializing structures, is therefore just a one-shot deal, unlike the kind of self-sustaining 
acceleration possible for truly general relativistic (i.e., nonlinear) backreaction effects. On 
the other hand, a never-ending collection of such one-shot contributions can by itself produce 
a persisting - and even growing - backreaction effect within V, as structures at ever-greater 
distances come into the causal horizon of V over time. 

We wish to come up with a simple representation of how the metric within V is altered 
by such one-shot perturbing contributions. An attempt by other researchers to estimate the 
total backreaction effects for a somewhat similar situation - for that of an infinite lattice 



of compact, static masses - was done by iGruzinov et al.l (120061 ). which once again found 



the main backreaction effect to be small, of order 0(H 2 l 2 /c 2 ), where I is the lattice size. 
But this result is only obtained as the leading-order backreaction term after the subtraction 
of an apparently divergent term due to the classical Newtonian gravitational energy of the 
masses, a 'divergence' which they assumed to be unphysical and thus removed via a "bare 
mass" renormalization step. But we have already seen this apparently infinite term above in 
Section 12.31 where we argued that the term was quite real, and rendered finite by the finite 
amount of time that it takes for inhomogeneities to initially form, and (more importantly) 
by the finite cosmic horizon out to which one can see those inhomogeneities. In a calcula- 
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tion using an eternal lattice of static structures, however, such a term would appear to be 
genuinely infinite and unphysical, and would thus unfortunately be dropped. 

In order to figure out how to properly implement Newtonian-level backreaction with 
causality in the absence of an exact cosmological solution, we do a little thought experiment. 
Consider our local volume V to be a homogeneous spherical region cut out of a matter 



domi nated, nearly perfectly homogeneous early universe. Now, it is well known (jWeinberg 



19721 . pp. 474-475) that the expansion evolution (i.e., the Friedmann equation) for V can be 
derived - using nonrelativistic Newtonian equations, in fact - without reference to anything 
outside of that sphere. Barring perturbations, the Friedmann evolution of V is determined 
completely internally, and all matter outside of it is gravitationally unmeasurable, as if the 
universe outside of V did not even exist. So let us remove it, leaving an infinite vacuum 
surrounding our spherical region V. 

This perfect isolation cannot last forever, of course, since later on, distant stabilized 
structures will form which certainly do exert real effects upon our local volume by becoming 
attractors which pull mass in towards themselves from all directions, including from within 
V. These gravitational effects, unknown to V before that time, will seem to appear as 
new from the 'empty' region outside of it, as if such clumped structures were brought in 
from infinity at some recent time during the structure formation epoch, to only then begin 
gravitationally affecting V. 

In most cases, one concerns oneself solely with the spatial gradients of the gravitational 
potential, <3>, caused by the distant perturbations; and so one is able to study overdensities 
and voids, bulk flows, and the like, even perhaps to the extent of trying to achieve cosmo- 
logical backreaction with them. But we have completely dropped those spatial gradients 
in our adoption of a smoothly-inhomogeneous version of the Cosmological Principle, since 
a randomly-distributed collection of inhomogeneities would exert forces in all directions, 
largely (if not perfectly) canceling each other out. 

But what does not cancel out, regardless of how the perturbations outside of V are 
distributed - because it is a non-directional scalar, rather than a vector - is their combined 
contribution to the overall level of that is, its actual magnitude, which does not matter in 
true Newtonian physics (only differences in potential do), but which does matter in general 
relativity. The contribution by all clumped masses to $ within V (i.e., $v) will always 
be negative (as gravitational energy is negative), thus adding together constructively and 
reinforcing their total effect upon our local volume, even for perturbing masses located on 
opposite sides of it. Furthermore, negative contributions to $v due to overdensities will 
not be canceled out by positive contributions due to underdensities and voids, because the 
situation is not symmetrical (and in fact is biased towards volume creation): underdensities 
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always keep expanding, and at a faster rate than the cosmic average; but overdensities do 
not keep shrinking, because they eventually stabilize themselves through vorticity and virial- 
ization. And these effects go far beyond the specific regions containing the clumped masses 
themselves, since by adding a nonzero net contribution to $v ; these exterior perturbations 
affect both the flow of time and the actual physical volume (consider the spatial metric terms 
of Equation [3]) within our local region of space. The fact that the effects of such distant per- 
turbations are individually small (oc 1/r) does not matter when one can sum the combined 
effects of all of the masses in spherical shells out to very large r. 

Our phenomenological approach is therefore one in which we model the inhomogeneity- 
perturbed evolution of V with a metric that contains the individual Newtonian perturbations 
to the potential <&v(£) from all clumped, virialized structures outside of V that have been 
causally 'seen' within V by time t, superposed on top of 'the background Friedmann expansion 
of V. Now, one may object that material outside of V actually falls within the observational 
horizon of V long before it became inhomogeneous, yet had no effect upon V whatsoever 
at that time; and thus such mass should not provide any new contributions to $v at a time 
later on, simply because it has spatially redistributed itself from being smooth into being 
clumpy But our reply would be that each individual structure - only after it has clustered 
and virialized itself - exerts its own individual pull upon the mass inside V, causing a 
peculiar-motion acceleration of the mass inside V towards itself, and that this is a new force 
(representing new "gravitational knowledge") not seen before within V. But one cannot feel 
a gravitational force, unless one simultaneously feels a gravitational potential; and so this 
perturbing potential must be new within V as well, approaching and entering V in causal 
fashion from this clumped object (and from all others) as they develop over time, everywhere 
in the surrounding universe. 

Now from the perspective of a general-relativist, the natural impulse would be to declare 
this entire procedure non-gauge- invariant, and then simply transform $v(0 away via a new 
time coordinate. But we assert here that this would be a physically improper procedure, since 
the act of making a cosmological observation is itself not a gauge invariant process. General 
relativity only guarantees a physical invariance under local coordinate transformations; but 
cosmological measurements are manifestly nonlocal - most importantly here in regards to the 
time coordinate - given that the Hubble curves used to demonstrate the existence of a cosmic 
acceleration are not measurements of a metric quantity at a single instant of time (or a single 
location in space), but are produced by integrating the motion of a light ray over its past light 
cone, an integration over timescales that are by design much longer than what could possibly 
be considered a 'local' range of t, in order to detect evolution of the cosmic expansion rate 
a(t). In cosmology, therefore, the time coordinate cannot be transformed to any other t'(t) 
time function with impunity, since the original t coordinate has a unique physical meaning: 
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it governs the rate of cosmic expansion (e.g., a(t) oc t 2 ^ 3 during matter domination), which is 
measurable in many ways. This cosmic time t thus serves as an absolute clock which we can 
use as a reference for comparison against the rate of local physical processes - for example, 
measuring how fast light rays travel through the universe, versus how long it takes the cosmic 
scale factor to double in size. Therefore, we claim that this $v(£) function - which affects 
local physics such as light-ray propagation, but (ironically, for backreaction) does not alter 
the cosmic expansion rate function a(t) - encapsulates real physics that produces observable 
consequences, and as such cannot be legitimately transformed away. 

The overall picture presented here is clearly a very heuristic one, and the growing per- 
turbation potential 3>v(£) would of course not represent the literal metric of the universe; 
rather, it is simply one that makes sense as a qualitative shorthand for the essential back- 
reaction effects generated by structure formation throughout the cosmos, as they act upon 
any particular patch of space. What is really happening, in a more complete physical sense, 
is that collapsing overdensities stabilize themselves and halt their collapse by concentrat- 
ing their local vorticity; this concentrated vorticity leads to real, extra volume expansion, 
representable (in the final state) at great distances by the tail of a Newtonian perturbation 
potential to the background FRW metric; and this Newtonian tail propagates continually 
outward into space by inducing inward mass flows towards the virialized object from farther 
and farther distances, as time passes. The total perturbation at any location in space (which 
will be independent of position, assuming similar structure formation rates everywhere) will 
then be the combined effects of innumerable Newtonian tails of this type, coming in all the 
way from cosmological distances, and from every direction. All that our thought experiment 
has allowed us to do is to show that it is reasonable to represent all of these complex physi- 
cal processes in a very simple way, as a cosmo logically uniform perturbation potential $v(£) 
superposed on top of the original FRW expansion, where the potential grows in time as 
gravitational information about virialized structures flows in from an ever-expanding "inho- 
mogeneity horizon" . It is the changing absolute level of this causally-developing gravitational 
potential function (referred to earlier as $sr(0> to emphasize special-relativistic causality) 
that we argue is the real main effect of backreaction; and we will show it to be sufficient 
for producing an observed acceleration completely on its own, in a smoothly-inhomogeneous 
context that does not require any other inhomogeneity-related effects. 



3.2. Inhomogeneity Evolution and Causal Integration for Metric Updating 



In order to derive an expression for the evolving perturbation function $sr(0> we start 
by considering the phase transition itself, during which an extremely symmetrical and smooth 
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( "unclumped" ) state during the early universe has gradually evolved to an almost entirely 
"clumped" state today During the transition, most of the matter in the universe is converted 
from something close to a completely smooth perfect fluid (say pressureless dust), to a 
discrete but infinite collection of randomly-distributed clumps of various sizes. We must 
define a physical quantity as a measure of 'cosmic dumpiness', to represent the extent to 
which this phase transition has proceeded to completion. One could perhaps use the ratio of 
the typical size of a stabilized mass clump divided by the typical distance between clumps 
as a symmetry-breaking parameter; or alternatively, one could define an order parameter, 
given as the fraction of cosmic mass in the clumped (i.e., virially-stabilized) state, versus 
that remaining in the smooth (i.e., freely-expanding FRW) state. We choose the latter 
approach; and as the order parameter grows large, this indicates the impending breakdown 
of the matter-dominated FRW evolution of a pure-dust model, and the growing relevance of 
inhomogeneities. 

To phenomenologically model the evolution of clustering, we therefore define a "clump- 
ing function", ^(t), defined (as a fraction of the total density) over the range from = 
(perfectly smooth matter), to ^(t) = 1 (everything clumped). Note here that we are assum- 
ing a spatially flat background universe everywhere in our calculations; and given that our 
formalism is specifically designed to eliminate the need for Dark Energy, this implies that all 
of the cosmic contents are treated as pressureless dust (for convenience ignoring relativistic 
matter and energy, which do exist but in lesser amounts), and that such dust adds up to 
Q M = 1. (This latter condition is only in reference to the unperturbed background, however, 
given that the apparent observational value of Qm will evolve to become quite different from 
unity, as we will see.) We will not be rigorously deriving \I/(t) from first principles, but will 
rather (as described in later sections) be utilizing a variety of functional forms and parame- 
ters motivated by basic physical principles and measurements, and then empirically evaluate 
their performance in light of cosmological observations. 

Obtaining a metric which can explicitly represent the full dynamics of a self-stabilizing 
structure, complete with vorticity and virialization, is an enormously complex task; it is so 
challenging, in fact, that we have already seen how exact solutions (like Swiss-Cheese models) 
and more general methods (like the Buchert formalism) have both failed in this task, due 
to the overly-restrictive assumptions that they were forced to adopt in order to simplify the 
situation enough to provide a tractable analysis. We therefore seek a practical way in which 
the key backreaction effects of clustering (described by clumping evolution function ^(t)) 
can be simply but effectively modeled, given our smoothly-inhomogeneous formalism with 
causal updating. 

We must begin by considering a single clumped mass (representing a virialized, self- 
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stabilized inhomogeneity) that is embedded within the expanding universe. It is actually 
quite easy to embed a single clum p within an empty, coasting universe by considering the 
"Expanding Minkowski Universe" ( [Robertson fc Noonanlll968l ) with a(t) oc t; this metric is 
mathematically equivalent to the (static) Minkowski spacetime of sp ecial relativity^ and c an 
be transformed to it via the transformation: t' = tyl + 



tr f Misner. et all 11973 



P- 

743, using (t',r') instead of their (t,x) coordinates). The inverse transformation is readily 
obtainable as: t = a/ (V) 2 — (r') 2 , r = r'/ a/ (t 1 ) 2 — (r') 2 (noting that the entire coasting 
universe is covered by just the t' > r' portion of the Minkowski metric). By using the 
simple trick of applying this inverse coordin ate transformation to the static Schwarzschild 
metric^ of a Black Hole (e.g., Weinberg 1972 ). rather than to pure Minkowski space, instead 
of producing a completely empty coasting universe we now obtain an (essentially) coasting 
universe which possesses a clumped mass at the origin. There is the drawback that such 
a cosmology is still empty everywhere outside the central mass, and thus has a negative 
global curvature and cannot represent a matter-filled cosmology; but a few appropriate 
approximations and ansatzes can readily convert this metric into a spatially flat, matter- 
dominated metric with a clumped mass at the origin. 

If one objects to such extrapolations, the n the exact sam e resulting metric can be 
obtained by linearizing the McVit tie solution (IMcVittid Il933l ). as can be seen from the 
perturbed FRW expression given in lKaloper et al.l ( 120101 ) . One thus gets, for the Newtonian 
approximation of a single clumped object of mass M, embedded at the origin (r = 0) in an 
expanding, spatially flat, matter-dominated (MD) universe: 

ds 2 « -c 2 [l + (2/c 2 )$(t)]dt 2 + [a M D(t)H^ , (7) 

where $(£) = {—GM/[aMD(t)r]}, and omd(^) t 2 ^ 3 is the unperturbed MD scale factor 
evolution function. Note that this McVittie solution actually includes no mass accretion - 
i.e., M in the above metric is a constant - even though such accretion would certainly be 
expected to occur in a realistic cosmology due to inflows onto any overdensity embedded 
within the matter-filled universe (and this was indeed one of our main arguments earlier for 
disregarding the Swiss-Cheese limits on backreaction). In any case, this "M" does not simply 
represent the total quantity of mass (i.e., conserved 'dust') within some specified volume, 
but rather embodies the amount of clumped and virialized mass in that volume (which grows 



6 Given the key role of spinning and vorticity in our discussion of backreaction, one might expect the use 
of the Kerr metric here; but the vast majority of mass in the universe does not possess relativistic amounts 
of angular momentum - merely a sufficient degree of it, acting continuously for a cosmologically-long time, 
to ensure that the object remains stabilized. Thus we treat the basic existence of persistently-stabilized 
Schwarzschild-like point masses as the biggest perturbation here, without needing to include the further 
metric perturbation which results from their actual rotation. 
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in time as a reaction to inflows onto overdensities); and so our phenomenological model will 
therefore effectively assume a time-dependent expression, M(t), in the above metric, where 
this time dependence is derived from the clumping evolution function ^f(t), as we now work 
out in detail. 

Assuming that various mass concentrations Mi will be randomly and (sufficiently) evenly 
distributed throughout the universe, one must integrate over distance (from a particular 
observer's location) for the various clumps, and one must also angle- average over them to 
get the metric for the combined effects that would represent an averaged ds 2 typical for that 
observer experiencing displacements in any given direction. (Angle-averaging the effects 
of an ensemble of randomly-distributed clumps will also cancel out any nonspherical Kerr- 
type behavior.) Since a displacement of magnitude \dr\ will have a radial component (with 
respect to any particular mass concentration) of \dr\ cos#, with a randomly-distributed value 
for the angle, the angle-averaging of the perturbation term therefore requires one to average 
over dr 2 = \df\ 2 cos 2 9 = (dx 2 + dy 2 + dz 2 ) cos 2 9. Since (cos 2 9) = (1/3) in three spatial 
dimensions, this means that the summed, averaged value of the spatial part of the overall 
metric perturbation due to a large number of randomly-located clumps will be multiplied 
by the factor (1/3); though the temporal part (i.e., the perturbation to g tt ) will remain 
unaffected by such angle-averaging. This is due to the fact that only the radial projection of a 
given translation (with respect to a particular mass concentration) will 'feel' the perturbation 
potential from that clump in its contribution to the interval ds 2 (cf. Equation [7]), while the 
full strength of the gravitational potential always contributes to the temporal part of the 
metric, regardless of directional configuration. 

Consider now the observer to be located at the origin, surrounded by a set of (roughly 
identical) discrete mass concentrations with total mass M, all located within a particular 
spherical shell at coordinate distance r' from the origin, but randomly distributed in (9, (ft). 
From the above arguments, we can write the total, linearly summed and angle-averaged 
metric for this observer in 'isotropized' fashion, as^: 

ds 2 = -c 2 {l - [R Scb (t)/r'}} dt 2 + [a MD (t)] 2 {l + (l/3)[R Sch (t)/r']} \df\ 2 , (8) 

where R Sch (t) = {{2GM / c 2 ) f[a MD {t)]}, and \dff = {dr 2 + r 2 d9 2 + r 2 sin 2 9d<ft 2 ) = \dx\ 2 = 
(dx 2 + dy 2 + dz 2 ). 

7 Note that this result for the overall summed perturbation is mathematically no different from what one 
would get by integrating a smooth spherical shell of continuous matter; but by thinking of mass M as being 
an assortment of many individual, vorticity-stabilized discrete sources, it becomes clearer to understand how 
the same actual matter - considering its mass plus its dynamical effects - can be responsible both for the 
background matter-dominated FRW metric, and for the perturbative metric terms superposed on top of it. 
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Looking at this expression, it seems physically reasonable to regard the term multi- 
plying \dr\ 2 as a 'true' increase in spatial volume, when comparing the volume of a spatial 
hypersurface at two different cosmic times; and to regard the g u term as an 'observational' 
term, slowing down the perceptions of observers - relative to the expansion of the universe, 
governed by cosmic time t - at all times after inhomogeneities have begun to develop (includ- 
ing at the current time, t = to)- Significantly, even if the spatial term by itself may not be 
enough to generate a real volumetric acceleration, once it is coupled with the temporal term 
the entire perturbation may indeed be enough to create a so-called "apparent acceleration" 
that is sufficient to explain all of the relevant cosmological observations. 

The factor of (1/3) in the perturbative term for the spatial metric components (relative 
to that for q tt ) may seem intriguing because it resembles the initial claim (later retracted) by 



Beanl (120091 ) of (I/77) = (ip/<p) = (^Time/^Space) — 3 from t heir analysis of the w eak lensing 



shear field of the Hubble Space Telescope COSMOS data (IMassey et al.l 120071 1 . appearing 
(briefly) to indicate a deviation from general relativity, which predicts ip = <fi in the absence of 
anisotropic stress. But this factor of (1/3) in our Equation[8]is not 'fundamental', but merely 
the result of approximating the linearized sum of many individual 'Newtonian' solutions; this 
effectively spreads the total spatial perturbation among all three spatial metric terms, rather 
than confining it solely to usual when considering a single inhomogeneity. Thus 

the general relativistic expectation of equal temporal and spatial potentials - i.e., ip — (j) - 
is not really violated here, and no actual new physics or modified gravity is implied by it. 

Implications aside, Equation |H] just represents the perturbations to the metric due to 
masses at some specific coordinate distance r' - and thus from a specific look-back time t' - 
as seen from some particular observational point. The total metric at that spacetime point 
must be computed via an integration over all possible distances, out to the distance (and 
thus look-back time) at which the universe had been essentially unclustered. Finally, a light 
ray reaching us from its source (e.g., a Type la supernova being used as a standard candle) 
travels to us in a path composed of a collection of such points, where the metric at each point 
must be calculated via its own integration out to its individual "inhomogeneity horizon"; 
and only by calculating the metric at every point in the pathway from the supernova to our 
final location here at r = z = can we figure out the total distance that the light ray has 
been able to travel through the increasingly perturbed metric, given its emission at some 
specific redshift z. 

Consider a light ray emitted by a supernova at cosmic coordinate time t = isN; which 
then propagates from the supernova at r = r SN , to us at r = 0, t = t . We refer here to the 
geometry depicted in Figure [2j 



For each point P = (r, t) of the trajectory, the metric at that point will be perturbed 
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away from the background FRW form by all of the virialized clumps that have entered its 
causal horizon by that time. Consider a sphere of (coordinate) radius a, centered around 
point P, with coordinates (a, t Tet ) (where t ret < t is the retarded time), defined such that the 
information about the state of the clumping of matter on that sphere at time t ret will arrive 

- via causal updating, traveling at the speed of null rays - to point P at the precise time 
t. To compute the fully-perturbed metric at P, we must integrate over the clumping effects 
of all such radii a, from a = out to a max , the farthest distance from P out from which 
clumping information can have causally arrived since the clustering of matter had originally 
begun in cosmic history. 

To really determine the relationships between (a, i re t), « max , and (r, t) with precision, 
we would need to compute the propagation of null rays from such (a,t ret ) to (a = 0,t) 
recursively or iteratively, since the propagation time of a null ray carrying perturbation in- 
formation would itself be affected by all of the other perturbation information coming in to 
cross its path from everywhere else, during all times prior to arrival. In other words, causal 
updating is itself slowed by the metric perturbation information carried by causal updating, 
creating an operationally nonlinear problem. But out of necessity for this initial, proof-of- 
principle paper, we avoid this "recursive" (as opposed to gravitational) nonlinearity through 
the convenience of assuming an unperturbed, flat, matter-dominated FRW cosmology for 
computing all causal updating effects for the metric at (r, t). This crucial simplification - in 
addition to the simplification of ignoring gravitational nonlinearities as well, since we com- 
pute the total metric perturbation terms via a trivial summation of individual contributions 

- is necessary for the purposes of this paper, though a more complete analysis (considering 
both nonlinearities) must eventually be employed in order to meet the standards of Precision 
Cosmology. For now, however, what we can do is to check at the end of calculations that 
these nonlinearities have not grown too severe as z — > (they should not be for most of our 
simulation runs, below), and to note that the output cosmological parameters and fits from 
our calculations here will unavoidably possess some systematic theoretical uncertainty as a 
result of these approximations. 

Now, for a FRW metric with a(t) = ao(t/to) 2 ^ 3 , the coordinate distance traveled by 
a null ray in the cosmic time span from t\ to t 2 will be a = (c/a ) (t/t )~ 2 ^ 3 dt = 

[(3c/a )(to) 2 ^ 3 (4^ 3 ~~ ^ 3 )]- E)efining ao = 3ct = 2c/ H , and with t 2 = t, ti = i rct , we 
thus have: a = [(t/to) 1 ^ 3 — (t rc t Ao) 1 ^ 3 ]- We then turn this into a prescription for computing 
t rct as a function of t and a (and implicitly of the present time, t ), as follows: 



t rct (*,<*) = t [(t/t ) 1/3 -a] 3 



(9) 
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Similarly, we can determine a max , given some initial time t^t & t which structure forma- 
tion can be reasonably said to have started (i.e., \l/(t < t init ) = 0): 



Now, how the metric is affected at P by a spherical shell of material at coordinate radius 
a depends upon the state of clumping there at the appropriate retarded time: ^[t Tet (t, a)]. 
The total effect is computed by integrating all shells from a = out to a = ct max (t, tinit); 
but in order to compute the metric perturbation from each shell quantitatively, it is first 
necessary to relate this clumping function to an actual physical density of material. 

As discussed above, we define the ^(t) function as representing the dimensionless ratio 
of matter which can appropriately be defined as 'clumped' at a given time, expressed as 
a fraction of the total physical density. With the convenience of assuming a flat FRW 
cosmology as the initially unperturbed state, the total physical density at all times will 
merely be an evolved version of the unperturbed FRW critical closure density from early 
(pre-perturbation) times. 

Recalling Equation[HJ we have the perturbation term [Rsch(t)/r'} = {{2G M / c 2 ) / [r' gimd^)]}, 
with omd (t) = ao(t/to) 2 / 3 = c[18t 2 / Ho} 1 / 3 . The value of M to use here is given by the 
clumped matter density at coordinate distance a, times the infinitesimal volume element of 
the shell. The clumped matter density at time t, as implied above, will equal (t) p CTit (t)]; 
and the volume element in the integrand for that shell is given by: 4tt R 2 hys dR p h ys = 
4vr[aMD(^) a] 2 [aMD(t) da]. (Note that any density-dilution effects in [\l/(t)p cr i t (t)] due to 
volume creation by virializing inhomogeneities will be precisely canceled by the correspond- 
ing volume increase of that spherical shell, leaving its (backreaction-effective) differential 
mass element, Ll dM" , unchanged. On the other hand, there should be an extra distance 
factor multiplying the denominator of [Rs<±.(t)/r'], due to $(£), -Rsch(^) 7^ 0; but we must 
drop it here in this simplified treatment in which we neglect "recursive nonlinearities".) 

Collecting these terms (and letting t' = t ret (t, a)), the integrand will thus be equal to: 



"max (t, tinit) = [(t/t ) 



/3 -(W* ) 1/3 ] • 



(10) 



r'=Q— ¥{a+da) 



{(2G/c 2 ) dM I [a MD (t) «]} 



(11a) 
(lib) 
(11c) 
(lid) 



{(2G/c 2 ) [auo(t) a] 1 Pcrit 

{(8nG/c 2 ) [a MD (t) a]' 1 [p crit (t) [a MD (t)] 3 ] [a 2 da]} 

{(8ttG/c 2 ) tf(t') aMD^ 1 [pcriM a 3 } [ada]} 



{(SttG/c 2 ) *(0 [{to/tf 3 (Scto)- 1 ] {[3H 2 /(8nG)] (3ct ) 
{12 *(f ) [(toA) 2/3 ] [ada]} , 



3 



} [c(dd$ 
(iif) 
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where for simplification we have used Ho = (2/3)fg 1 an d the fact that [p(t)a(t) s ] is constant 
(or effectively so, as explained above), both true for a matter-dominated universe. Note 
also that only \1/ is evaluated at the retarded time, t ret (t, a). The strength of the metric 
perturbation (at time t for point P) for a point-like Newtonian perturbation embedded 
in the expansion actually depends upon its instantaneous physical distance from P at t, 
as is obvious from $(t) = { —GMj [omd (f ) r]} ~ [-GM/i? phys (t)] in Equation [3 The 
only "relativistic" piece of propagating information which is causally delayed is the state 
of clumping, \P[f ret (f, a)], that has just then arrived from coordinate distance a to observer 
P at (r,f). 

From this result, we can now determine the total integrated metric perturbation function 
due to clumping, J(t), as experienced by a null ray (or any observer) passing through point 
P at (r,t): 

I(t) = / {12 ^[t iet (t,a)} [(t A) 2/3 ]} « da , (12) 

Jo 

with 1(f) implicitly being a function of t init (with I(t) = for f < ti n ; t ), as well as of f . 

Finally, we can insert this result back into the formalism of Equation [HJ to obtain 
the final clumping-perturbed metric that we will use for all of our subsequent cosmological 
calculations: 

ds 2 = -c 2 [l - 1(f)] dt 2 + {[a MD (f)] 2 [1 + (l/3)/(f)]} \df\ 2 , (13) 

Representing a smoothly-inhomogeneous universe, as it does, this metric depends only upon 
time, and is thus equally good anywhere in the modeled universe - in particular, at every 
point in the trajectory of a light ray from a distant supernova to us. 

A few important comments must be made about this result. First, note that the inte- 
grand for 7(f) in Equation [12] would actually become infinite as a — > oo, were it not limited 
by the finite causal horizon for seeing perturbations (i.e., finite a m ax), and by the lessened 
degree of clumping as one looks back to earlier retarded times, deeper in the past (i.e., 
^(f — > fi n i t ) ~ 0). This implies that the integrated result for I(t) can indeed become quite 
large at late times for Newtonian-level perturbations alone, given a large degree of clumping, 
in accord with our discussion from Section 12. 3j and also that the dominant contribution 
to I(t) will typically come from the largest coordinate distance out to which one can still 
see a substantial degree of clumping at its associated f ret , indicating (as we will see) that 
clustering models \l/(f) with stronger clumping earlier on will have a much more powerful 
overall perturbative effect. 

Second, we should reiterate the various approximations that have been made in or- 
der to obtain this smoothly-inhomogeneous cosmological metric - beyond, of course, the 
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core assumption of spatially-random clustering. These include the gravitationally-linearized 
treatment in the summing together of independent metric perturbations (only valid for I(t) 
sufficiently less than unity, as will be examined later in Section [6.ip . and the dropping of 
the recursive nonlinearities inherent to the physical process of causal updating. Such effects 
would eventually need to be included to produce a high-precision model of this type of back- 
reaction, from individually-Newtonian, virialized structures. But too-sophisticated a model 
along these lines does not necessarily make sense, since at some point the basic assump- 
tion of spatial randomness itself breaks down, requiring the abandonment of this smoothly- 
inhomogeneous formalism e ntirely in favor of a fu lly 3D cosmic structure simulation program 



- perhaps along the lines of ISpringel et al.l ( 120051 ). for example, with Newtonian- level backre- 



action effects and causal updating added in. Since such a 3D simulation model is far beyond 
the scope of this paper, we will consider it sufficient here to stick with Equation [T3] as a 
reasonable first-order approach to the problem, with all approximations and caveats kept in 
mind. 

Lastly for this subsection, we note that the metric in Equation [13] has been made to 
look very much like an ordinary FLRW metric; and it can be made to look exactly like 
one with the transformation dt' = \/l — I{i) dt, and with an appropriate redefinition of 
the scale factor, a' {if) = {cimd yH- {l/3)I[t{t')}}. One might then be tempted to 
conclude that observing an acceleration in a' (if), just like in the usual FLRW case, still 
requires some form of Dark Energy violating the strong energy condition. But this would be 
incorrect, because the cosmological model that we have produced here is not simply a true 
FLRW metric shown in different guises via coordinate transformations, but in fact is a very 
physically different (dynamically inhomogeneous) model that is merely being approximated 
with FLRW-like averaged parameters. The detailed astrophysical effects being averaged into 
it, in contributing to the 'real' volumetric effects in g rr and/or to the 'observational' signal- 
delaying effects in g tt , are themselves capable of combining together sufficiently to create an 
apparent acceleration in a'(t') - even without any SEC- violating component, or any actual 
accelerating patch of spacetime as might be measured by observers in a local reference frame. 



3.3. Light Propagation, Redshifts, and Luminosity Distances with Causal 

Updating 

Due to the altered mathematical form of Equation [13] with respect to the pure FLRW 
case, the observed cosmological evolution and parameters will differ from what would nor- 
mally be expected given some usual a(t) function. One must distinguish between parame- 
ters which merely refer to the underlying theoretical FRW (unperturbed, no Dark Energy) 
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model that holds for the pre-clumping universe - that is, the "bare" parameters - versus 
those ("dressed") parameters which are obtained from current-day observations. We will 
use the superscript "FRW" for the former, and the superscript "Obs" for the latter, and 
will compute the relationships between these different parameter sets for observables like 
redshift, the Hubble Constant, the Age of the Universe, and other important variables. 

We begin by focusing upon parameters used for computing the supernova-based lumi- 
nosity distance function that most directly traces out the cosmic expansion history. First, 
one must calculate how observed redshifts are altered in this model. Looking at Equation [TBI 
which is in the form ds 2 = —gtt{t) dt 2 + g rr (t) \dr\ 2 , we note again that one could transform 
away the gttif) function with a redefinition of the time coordinate, and thus redshifts can be 
calculated here simply by taking y/ g Tr (t) as the new scale factor. In other words, given the 
standard relationship for the "FRW" variables: 

^FRW (t) = _ 1 = (Wt)2/3 _ 1 ^ (14) 

we can similarly write: 



^Obs 



(t) = _ ! = [ /l + (l/3)/(to) 2/3] _ x (15) 



where we define to = ^o RW > ^mit = ^£it W j an d we wm mean t x = t FRW , z x = ^ RW , for time 
coordinate and redshift values in general, except when expressly referring to them in the 
form i° bs or z® hs . 

Now, in order to calculate observed luminosity distances, we must compute the coordi- 
nate distance r of a supernova going off at coordinate time t, for which a light ray would 
be arriving here (at r = 0) precisely at to- The modification of r(t) imposed here as a 
perturbation to the FRW result, essentially due to the summed Shapiro time delays (e.g., 



Weinberg! Il972l ) contributed by all causally-seen virialized clumps, is the major effect of in- 



homogeneities th at we consider in this paper, since we neglect other effects like lensing along 



beam paths (e.g., iKantowskil 120031 ) . and so on, which are likely too small to generate an ob- 
served cosmic acceleration. For a null ray, with ds 2 = 0, and considering pure inward radial 
motion (\df\ 2 — > dr 2 , dr/dt < 0), we have dr/dt = — -\/ gtt(t) / g r r(t) ■ We can thus compute 
the coordinate distance of the supernova from US, cLS 8b function of t, as follows (with I(t) 
still as computed via Equation [T2~]) : 




,KRW (t) = | r FRW (to) _ rF RW (t)|= I { X l^±}dt> (16a) 
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*°{ — C -7-^ \ 1 , f wl a ) (16b) 



aoit LV 7 ' V 1 + (V3)/(f) 



This coordinate distance function can then be converted into an expression for the 
observed luminosity distance. For the perfectly homogeneous FRW case, the luminosity dis- 
tance of a standard candle is given by the current physical distance to it, times a redshift 
factor: g?l,frw = [ a o r SN (1 + z)\. For our inhomogeneity-perturbed model, all of the ap- 
propriate time dilation/redshift factors are taken care of by using z° hs (t) from Equation [T5| 
and the physical distance can be given by ^^(t) and the modified scale factor, as follows: 

4,Pert(*) = Ml + (l/3)/(t )] r^ w (t) [1 + z 0bB (t)] (17a) 
1 + (V3) 




'i + mm 

1 + ( V3) 



>/l+[/(*r)/3] 3 k " T> V 1 + Wr)l*\ 



(17b) 
(17c) 



where Iq = I (to), and t r , t' r are dimensionless time ratios (e.g., t r = t/to), with no change 
to the essential form of I(t) (i.e., I(t) = I(t r -t Q ) =>• I(t r )). Note that this o?L,Pert(^) is simply 
proportional to (ct ) = [(2/3) c/Hq RW ], as would be expected. 

With this expression for c?L,Pert(^), and Equation[T5]for 2;° bs (t), one could in theory com- 
bine them analytically to produce c?L,Pert(^ 0bs ), the function actually observed in supernova 
luminosity distance (i.e., Hubble) curves. But this is neither analytically nor computation- 
ally practical, and so we have instead performed our numerical calculations using arrays for 
many discrete points in t, utilizing the above formulae to evaluate arrays for G?L,p e rt(X) an d 
z (t), which we then simply combine together into one array as c?L,Pert(-2 ;0bs )- It should 
therefore be noted that all of our subsequent plots of luminosity distance curves in this pa- 
per for different clustering evolution models, though presented as smooth curves, are in fact 
connected dots of points evaluated for discrete values of t FRW . 

Our pixelization, which we have rigorously tested to ensure accurate results, typically 
employs ~2150 discrete time values, most of which (~1500) are concentrated in the more 
recent times (z FRW < 1) where most of the dynamical evolution is happening; fewer pixels 
are reserved for earlier times, going back to before clumping had started when the simple 
FRW model was still correct. Increasing the number of discrete time values used, even by 
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a factor of 5, results in only a small change (e.g., < 0.15%) for all cosmological fits and 
parameters, even those requiring up to three derivatives of the luminosity distance function. 

Given that this discrete version of c^PertO 2 ) must in fact be differentiated to obtain 
cosmological parameters from it (see Section [5] below), we do so by using the definition of 
the derivative for each pixel. That is, to get the i th pixel entry (i.e., evaluated at t^y) for 
the N th derivative of expert with respect to z° hs , we compute it as follows: 

rWV/ 1 - _r_I r l( N - 1 )'] - Q L,Pert - »L,Pert 

L^PerdW ~ d zQhs L fl L,Pert J« " ^Obs (t{i+i}) _ 2 Ob s(t{ . }) " W 

Note that each subsequent derivative array has one less pixel than the prior derivative array, 
since we lack an extra pixel at the (high-z) end to compute the last pixel of the differentiated 
array (e.g., if d' LPert has 2149 pixels, then d'l Pert has 2148 pixels, and so on). 

The calculation of currently observable cosmological parameters requires derivatives 
which are evaluated at z = 0, t = t . We simply use the first (lowest-z, highest-t) pixel for 
that derivative value: 

[^PeJo^O, t-+t ) = [d^Pertlfl} ' ( 19 ) 

This appears to be a robust procedure, since the arrays of derivative values are well-behaved 
everywhere (except for transient jumps at transitions of our piecewise-continuous input \I/(t) 
functions at high-2;), and are smooth heading towards t — > to- Moreover, many simulation 
runs with a variety of different clumping function time-dependencies and amplitudes have 
given results for these cosmological parameters which seem to be reasonable and mutually 
consistent. 

Ultimately, luminosity distance functions for standard candles are plotted as residual 
Hubble diagrams, where one logarithmically plots the ratio of g?l to the luminosity distance 
function for a coasting universe of the same current expansion rate. More precisely, we 
are interested in the function A(m — M) = 5{Log 10 [<iL,Data(-2 ; )] — Log 10 [<iL,coast( 2 )]}, where 
^L,Coast(-2 ; ) = \{. c /Hq){\ + z)Ln(l + z)]. But, which H does one use here for subtracting off 
the coasting universe from our numerically simulated perturbed-universe models? Although 
we see from Equation [T71 that expert is indeed proportional to to oc (l/iJQ RW ), that is not 
sufficient to tell us the real observed expansion rate that is asymptotically approached as 
z — > (i.e., H° hs 7^ if FRW ), since the integral expression could evaluate to practically 
anything, given the appropriate clumping evolution function. 

T o determine the prope r value to use for the z — > expansion rate, consider Equation 
2.51 of Kolb fc Turner! ( 119901 ). from which one obtains the general approximation d^ ~ cz/H 
for very small z. This gives an operational definition for the observed expansion rate, H{? hs , 
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via: 

< Pert ,o = {d[d L ^ t }/d[z Ohs }} z ^ = c/H^ . (20) 

(We will later use this general method to explicitly relate H® hs to i/ FRW , and to find similar 
relationships for other cosmo logical observables, in Section [5j) Using this prescription, we 
can finally define the residual Hubble diagram function for our clumping-perturbed model, 
as follows: 

A(m - M) Pert (z° bs ) = 5 { Log 10 [d L , Pert (z Obs )] - Log 10 [« )Pert)0 )(l + ^° bs )Ln(l + ^° bs )] } . 

(21) 

A straightforward series expansion demonstrates that we have correctly specified the formula 
for the observed expansion rate, i?° bs , in order to match the expansion velocities of the 
coasting and clumping-perturbed models at to- 

Using Equation [5T] and the rest of the expressions in this section, we are now able to 
convert any clumping evolution function \l/(t FRW ) into a residual Hubble diagram that can 
be compared with any theoretical FLRW model that one chooses (e.g., ACDM models), as 
well as with real standard candle data from Type la supernova observations. 



4. MODELS, PARAMETERS, AND RESULTS 

4.1. Selection of Clumping Evolution Functions 

In order to produce cosmological predictions with our formalism, one must first choose 
a set of \l/(t FRW ) functions to evaluate which are reasonable models of how the fraction of 
cosmic matter in the clumped state has evolved over time. Choosing likely functions is not 
as trivial as it may seem, however. For example, while the contrast of a density variation 
will evolve as Sp/p oc a(t) oc t 2 ^ 3 in the linear regime for a matter-dominated universe (e.g., 



Kolb &: Turner! 119901 ). this only represents the linear evolution of a single clump; it says 
nothing about the initial development of new clumps (often due to collisions), or about the 
nonlinear regime and virialization for very dense clumps. Our clumping evolution functions, 
however, must serve as proxies for emulating all of these effects combined together. 

As one possibility, one could perhaps look at clumping evolut ion functions that ar e out- 



put from large cosmological structure formation simulations (e.g.. ISpringel et al.ll2005l ). But 
since such models basically use the instantaneous Poisson Equation formalism (recall Equa- 
tions Hl|6]), without causal updating, they lack the realism necessary to properly simulate the 
evolution of clustering without the use of a Dark Energy "fudge factor" , which complicates 
the interpretation of such results. 
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This author's first approach in defining \I/(t FRW ) was to consider observational data 
directly, for guidance, treating the star formation rate (SFR) as a tracer of the rate of 
the increase in clumping - i.e., d[*f>(z)]/dz oc SFR(z). This method of choosing likely 
functions was indeed capable of finding some which reproduced the supernova Hubble curves 
with reasonable success (results not shown here); but giv en the large uncertaint ies in the 



observationally measured SFR power-law parameters (e.g.. iGlazebrook et al.ll2003l ). and the 
large number of arbitrarily-tunable parameters in the ^{t) functions adapted from such data, 
the actual statistical significance of a 'good fit' was too difficult to meaningfully determine. 

Given these difficulties, we decide here to opt for simplicity in choosing which clumping 
evolution models to use for the main runs of our numerical simulation program in this analy- 
sis. Obviously, there is nothing simpler to use than $?(t) oc t, a fairly sensible choice to begin 
with, since one would assume that the amount of clumping that can occur should to some 
degree depend directly upon how much time is available for that clumping to develop. For 
some alternatives, we have also chosen to use models with \l/(t) oc t 2 / 3 , as this is proportional 
to the linear density contrast evolution in a matter-dominated universe; as well as models 
with an 'accelerating' clumping rate, \l/(t) oc t 2 , to test whether that would possibly help in 
creating an observed acceleration. This latter time-depen dency also potentially corresponds 



to the final nonlinear evolution of a density perturbation ( IKolb fc Turner! Il990l p. 322) 



Quantitatively, we define our three different classes of clumping evolution models as 
follows: 

= * (^FM (22a) 

to ~~ Hnit 

*md(*) = ^o(^^) 2/3 (22b) 

to ~ Hnit 

W*) - *o(^y^) 2 , (22c) 

where t init represents the beginning of clumping, such that \l/(t < i init ) = for all models; 
and \I/o = ^(t Q ) represents the current state of clumping today. Note that all of these 
functions are defined in terms of t FRW (and hence z FRW ), rather than in terms of t° bs or 
z 0hs ; the latter would of course be preferable, though it is not possible here because t° hs and 
z obs depend recursively upon *&(t). 

For these three different classes of models, we have two physically meaningful parameters 
to vary: \& an d tinit (though we will usually express the beginning of clumping in terms of 
Zinit, which can be obtained from t init via Equation Hlj) . Now, while this gives us a fairly 
wide region of parameter space to explore through in trying to match the observed supernova 
data, results that succeed in matching the data are still meaningful, since these model input 
parameters are constrained by astrophysical considerations; and also because (as will be seen 
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below) these classes of models have characteristic behaviors that quite naturally look very 
much like ACDM cosmologies over a wide range of redshifts and model parameter choices. 



4.2. Numerical Results: Residual Hubble Diagrams and Supernova Data Fits 

In designing a suite of simulation runs to test the effectiveness of our model at repro- 
ducing the observed cosmic acceleration, we use established observational data as our guide 
for specifying interesting choices for parameters z init and \l/o- 

With Zinft representing the beginning of the cosmic phase transition from smooth to 
clumped, it seems reasonable to associate z init with another symptom of the beginning of 
structure formation: the epoch during (or slightly before) the onset of cosmological reion- 
ization. 



The five-year WMAP Data analysis (IDunkley et al.ll2009l ) supported a value of z 



11 in the case of an instantaneous reionization; but the data also suggests the possibility of an 
extended period of partial reionization, perhaps beginning as early as z ~ 20, and extending 
no later than z ~ 6. To broadly cover this range (and to bracket it, to be conservative), we 
have chosen values of z^ = (5, 10, 15, 25) for our simulation runs, with the larger values of 
Ziait generally being more astrophysically interesting as starting times for the transition to 
clumping. 

Specifying appropriate values of is a more subtle task, though, since one has to 
quantitatively characterize the much more complex dynamical situation of late-time clus- 
tering with a single number. Furthermore, we must rely a great deal upon our assumption 
of a smoothly-inhomogeneous universe, since the existence of a local void or bubble would 
seriously alter the limiting behavior of ^f(z — > 0). Nevertheless, we are able produce a range 
of parameter values which make general astrophysical sense, and also turn out to produce 
good cosmological results. 

First, consider that while observations may tell us that, say, fiM bs = 1 — Q® bs ~ 0.27 
with n£ hs ~ 0.04 (and thus figjg ~ 0.23), our model here is one of a flat, apparently 
accelerating universe, with no dark energy. Thus we exploit the fact (details to be given in 
Section I5T21 below) that if° bs ^ Hq RW , in order to achieve f2^ RW = 1 without ever changing 
the actual physical matter density (cum oc pm) that can be more directly determined through 
other observations (growth of structure, cluster mass measurements, etc.). Furthermore, 
since the physical baryon den sity uib, and its rela tionship to wdm, are fairly well tested by 



the CMB peak height ratios ( IBennett et al.l 120031 ) regardless of whatever model is used to 



produce acceleration at later times, it is therefore wise not to change this ratio, o^/cudm- A 
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simple scaling up to flatness (fi^ w = 1) thus gives us 1^ RW ~ [0.04(1.0/0.27)] ~ 0.15, and 
fi™ w ~ 1 - n F b RW ~ 0.85. 

We must then determine reasonable estimates for how much of each species, Dark and 
baryonic matter, might be clumped at the present time, noting (as per Section 1372]) that 
is given simply as a dimensionless fraction of the total matter density. 

At one extreme, we may treat the universe as 'completely clumped' at the present time 
-i.e., = 1. 

At the other extreme, we may consider almost all baryonic matter to still be unclumped, 
due to 'gastrophysics' like shock heating and its resultant thermal pressure; exceptions being 
only that amount of baryonic matter clumped into stars, and perhaps (depending upon how 
to appropriately define 'clumped') the amount of gas bound into virialized galaxies. Dark 
matter, on the other hand, not being subject to ordinary thermal pressure (and being able 
to viriali ze gravitationa l ly), w ould be almost entirely clumped (for Dark Matter clustering, 
see, e.g., iGilmore et al.l (120071 )). except for some small portion of it that would actually be 



Hot Dark Matter (neutrinos). 

Gett ing prec i se num bers for these quantities is not trivial, but estimates are available: for 



example, iTurnerl (l2002bl ) estimates ~l/8 of the baryonic matter to be contained in stars, so 



that (for our model) f2 sta rs ~ 0.0185. Also. iHinshaw et al.l (120091 ) limits the neutrino physical 



density to being less than either ~12 — 13% of the total dark matter density (WMAP data 
only), or ~5— 6% of it (WMAP+BAO+SN combined data), which corresponds (respectively) 
in our model to < 0.1 or Q u < 0.05. 



Thus, summarizing the total density budget, one ends up with ~2% locked up into 
stars, ~13% remained as either clumped or unclumped baryonic gas, < 10% (or < 5%) as 
unclumped neutrinos, and > 75% (or > 80%) as mostly clumped Cold Dark Matter. 

We can therefore classify ~77 — 82% of the mass as 'probably clumped', < 5 — 10% as 
'probably unclumped', with ~13% or so as some mix of both. With this information, and to 
space out our parameters fairly evenly, we have chosen values of ty = (0.78, 0.85, 0.92, 0.96, 1.0) 
for our simulation runs, with the mid-range values of \l/o likely being the most astrophysically 
sensible. 

With four different values of z init , five different values of ^q, and three different clumping 
evolution models, this gives us 4 x 5 x 3 = 60 simulation runs in total. Residual Hubble dia- 
grams have been computed for all of these runs, and are plotted below, with Figure [3] showing 
the results for the runs, Figure @] depicting the \l/sqr runs, and Figure[5]depicting the \1/md 
runs. Each figure includes four panels, with the panels representing z init = (5, 10, 15, 25) in 
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order from top to bottom. Within each panel itself, the FLRW cases of the flat SCDM model 
(f^M = 1) and a Concordance ACDM model (Q\ = 0.73 = 1 — Qm) are shown for comparison 
against five of our models with different degrees of clumping: ty = (0.78, 0.85, 0.92, 0.96, 1.0) 
in order from the lowest curve to highest. 

A detailed quantitative analysis of the best of these runs in terms of quality of fit to 
real supernova data, and with several important cosmological parameters being computed 
for each run, will be presented below in Section [6J For now, though, we make the important 
qualitative observation that this formalism works. Not only does it yield a selection of several 
models possessing a perturbative effect strong enough to reproduce the observed (apparent) 
acceleration, but it produces curves that clearly do behave very much like ACDM - in 
particular, the ^sqr runs look like flat ACDM with f2 A ~ 0.3 — 0.4, the runs look like 
flat ACDM with fi A ~ 0.5-0.8, and the # M d runs look like flat ACDM with fi A ~ 0.65-0.97. 

(Ironically, the 'accelerated' clumping models, \I/sqr, produce the weakest observed accel- 
eration effect, because they have less clumping at early times. Clumping is more important 
the earlier it occurs, because of a geometric effect: a longer look-back time for the beginning 
of clumping results in a much larger horizon out to which an observer can see substantial 
inhomogeneities; and it is the huge volume of this outer shell of perturbations that produces 
the strongest effect on the observer, as pointed out towards the end of Section [3721 ) 

Most importantly, even without doing a search over our parameter space for 'best-fit' 
models - just by choosing a set of simple clumping evolution models and astrophysically- 
reasonable parameters for input into our model - we find that a significant number (~10 or 
so) of these 60 simulation runs manage to fairly precisely reproduce the Concordance ACDM 
Hubble curve. Furthermore, as will be shown shortly, these runs are able to fit the supernova 
data essentially as well as (and in certain cases, with some model parameter optimization, 
even better than) such best-fit flat ACDM Hubble curves that have heretofore been used to 
argue for the existence of Dark Energy. 



5. FORGING A NEW CONCORDANCE FOR A 
SMOOTHLY-INHOMOGENEOUS UNIVERSE 

The qualitative ability of our formalism to reproduce ACDM-like Hubble curves was 
shown in the previous section, and a more quantitative analysis demonstrating how these 
models are able to fit the observed supernova data will be given later on. 

But beyond just succeeding at explaining the apparent acceleration seen in data like that 
from Type la supernovae (SNe), it is widely recognized that a true model of the universe must 
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satisfy the constraints imposed by several different, complementary cosmo logical data sets, 
while simultaneously producing a set of cosmological parameters that are consistent with 
all other relevant astronomical observations. Only such a fully consistent, astronomically- 
correct cosmological solution would attain the status of a 'concordance', sufficient to replace 
the well-known "Cosmic Concordance" representing the range of ACDM models that appear 
already to be cosmo logically consistent, given the state of the data at this time. 

In order to extract the appropriate cosmological parameters from our simulated Hubble 
curves, we must compute each of the relevant cosmological observables from G?L,Pert(<2 0bs ) as 
it was defined previously in Equation [17] (and subsequent discussion), with all derivatives 
(with respect to z° hs ) and limits of this (discrete) simulation output array performed as 
indicated via Equations [TBlfT^l Brief derivations will be presented below of the expressions 
needed for converting model simulation results into observable parameters. 

With such expressions in hand, we will be able to determine how these ("dressed") 
observables relate to unperturbed ( "bare" ) model parameters such as to RW , -^o RW > an d ^m RW > 
for different choices of z init , and clumping evolution model ^un, ^md, or ^sqr- These 
relationships are interesting not only because they give us hard numbers to use for comparison 
with real observations, but also because they reveal how much causal backreaction due to 
structure formation has altered the evolution of our universe from that predicted by purely 
homogeneous FRW models. 

While it is beyond the scope of this paper to construct a completely new concordance 
in all of its aspects, in the subsections below we will show how to calculate several impor- 
tant cosmological parameters from our numerical simulations. Then, in Section [6j we will 
demonstrate the observed consistency of our best-fitting models with several of the key ob- 
servational parameters of the ACDM Cosmic Concordance, without the use of any negative 
pressure species like Dark Energy. 



5.1. Hq , Cosmic Proper Time, and the Age Problem 

As discussed earlier in regards to Equations |2"0H2~T| an operational definition for the ob- 
served Hubble constant in terms of our simulation results can be given as if ° bs = c/d' L Pert . 
For a normal, matter-dominated flat SCDM cosmology with no Dark Energy (and neglecting 
radiation), one has if FRW = [a(t) / a(t)] t ^ to = (2/3)tQ 1 . We can put this together to get the 
straightforward result: 

9 1 

rrFRW _ rrObs r A' \ (OQ\ 

11 ~ n Q I g C £FRW u L,Pert,0J • l zc V 
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Now, we recall from Equation [TTb that G?L,Pert(£) (and hence o^Pert^ 0138 ) , and all of its 
derivatives with respect to z° hs ) are simply proportional to ct FRW ; thus any dependence upon 
the parameter t FRW cancels out (as it must) from Equation [221 and this formula merely gives 
us a dimensionless ratio between the observed and unperturbed Hubble Constants, H® hs and 
if FRW , the value of which depends upon the result of that numerical integration. 

Note that it is if FRW which we place on the left hand side of Equation [231 as the 
'unknown' parameter; the value of H® hs is set from real Hubble recession measurements 
(getting, for example, 72 km s _1 Mpc~ 1 ), which one can then translate into if FRW for any 
particular clumping evolution model, in order to obtain its 'true' cosmic expansion rate, as 
would have been observed if perturbations had never altered the observed expansion rate 
from its FRW value. 



The importance of this relationship is that the expression in braces in Equation [23] will 
be less than unity (since I(t) > 0) when perturbations exist, thus resulting in if FRW < H® hs . 
This permits one to have a low value of if FRW (say, in the 40's), while still retaining 
H® hs 72. This possibility explains why several cosmological measurements may be 
quite concordant with a low Hubble constant, while removing the contradiction that such 
a result would seem to create for late-time measurements with standard candles (e.g., 
Freedman et al.l 120011) th a t clea r ly indicate a high H ® hs . (For low-ffo discussio ns, see for 



example Blanchard et al. ^2003 



. ISpergel et ajj (120031 ) , and lHunt fc Sarkarl (120071 ); and also 



see Figure 14 of lLarson et al.l (120111 ). showing the consistency of WMAP-only data with 
Qyi = 1 models for ff in the ~30's-40's.) 

Having if FRW ^ H° hs also makes ot her apparent conflicts go away, such as the classic 



Age Problem/Crisis in cosmology (e.g., iKolb fc Turnerl Il990l ; iTurnerl l2002bl ). in which a 
matter-dominated SCDM universe appears to be younger than some of its oldest constituents 
(e.g., globular clusters). The age of such a universe is t = {2/?>)Hq 1 , which for ff ~ 
72 km s~ 1 Mpc _1 gives only t — 9 GYr, requiring one to assume an accelerating universe 
that had slower expansion in the past, in order to lengthen to to ~13 — 14 GYr. 

For our model, on the other hand, we can use the metric given above in Equation [TBI to 
relate the 'observed' age of the universe, t® hs , to t FRW , as follows: 



t 



Obs 



.FRW 

JO 



(24) 



Now, even though this results in t$ hs < t FRW , we also have if FRW < if Obs , so that 
the value of t FRW = (2/3)(l/if FRW ) will be much larger than that expected from FRW 
considerations (i.e., significantly larger than ~13 — 14 GYr); and thus when these two factors 
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are combined together, the result for t® is able to fall precisely within the range necessary to 
solve the Age Problem, as we will show below in Section IfTTl for several of our best simulation 
runs. 



5.2. Spatial Flatness and the Observed Matter Density: f2^ RW versus f2^ bs 

Another major argument used in favor of Dark Energy as part of the Cosmic Concor- 
dance, is the apparent contradiction between CMB peak data showing the universe to be 
spatially flat (f2 Tot ~ 1), while actual searches for the required amount of matter persistently 
turn up short on the overall density of clustering mass, looking instead like pm/ Pait — 0.3. 
The usual conclusion is that Qm — 0.3, and that the gap of (0x o t ~ ^m) — 0.7 is filled by 
the existence of Dark Energy. 

The unstated assumption in this reasoning, however, is that it is always accurate to 
use Qm = ^M/h 2 = PM/Pcrit to relate the closure density value of the matter to its actual 
physical density. But because if FRW ^ Hq, and thus p FRW PcriT> this ceases to be true. 
Specifically, oc um/ (Hq rw ) 2 may very well be equal to unity despite the low value of the 

physical density, u M = [p M (8^/3) (100 km s^Mpc" 1 )- 2 ], when # FRW < H° hs is properly 
taken into consideration^. This is important because the spatial flatness of the universe - 
particularly as determined using data observed from the ancient and very homogeneous CMB 
epoch - is dependent upon the value of the unperturbed parameter, ^m RW ; relevant to the 
very early universe; not upon the observationally-defined parameter, fly^ s oc WM/(-f^o >bs ) 2 , 
reflective of the more recent, post-structure- forming epoch. 

To obtain an expression for the relationship between fi^ bs and fi^ RW m our smoothly- 
inhomogeneous cosmological formalism, we begin by assuming that the universe actually 
is spatially flat in terms of its FRW-defined parameters, and would have appeared to be 
so to observers in the past, before the onset of the apparent acceleration due to clump- 
ing. Now, at such an early time, t E - say, £cmb <C t E "C t FRW , f° r matter-domination, 
but yet small inhomogeneous clumping - one had (as functions of t E ): p M ,E = Pcrit^E] = 
3(iJ FRW [tE]) 2 /87rG = (l/67rG)(tE)" 2 - As the universe evolved to the current epoch, the 
volumetric dilution of matter went like g r -r^ = {[GM D (t)r 3 [1 + (l/3)/(t)]" 3 / 2 } (cf. Equa- 
tion [13]) . yielding a matter density today of: pm,o — Pm,e {(^eAo^) 2 H + (W3)]~ 3 / 2 } = 



8 Note that we assume the physical matter density from observations to be accurate despite causal back- 
reaction - i.e., Pm RW = /°M bs = Pm ~ since it can be measured at lower redshifts and in less cosmologically- 
dependent ways than Clw Thus we treat = = wm as unchanged in our formalism from the usual 

FLRW value, defining all of the discrepancy to be within HQ hs and ^M bs - 
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{(1/6ttG)(C W )- 2 [1 + (/o/3)]- 3 / 2 } = {[3(if FRW ) 2 /8vrG] [1 + (V3)]" 3 / 2 }. 



To turn this formula into a value for fi^ bs , one must consider it in terms of what we 
think the critical density is today, observationally. Now clearly, p° b t s = 3(if^ )bs ) 2 /87rG ! ; and 
simply by taking the ratio ^m* 3 = /°M,o//°SttO! we § e ^ the following result - along with an 
alternative way of expressing it, including a convenient definition for f2™ w : 



Now, the proper way to view these relationships is as more of a consistency check, 
than as an independent prediction. If one assumes an initially-flat FRW universe, and 
adopts some favored inhomogeneity clumping evolution function *&(t), then one can use 
it to compute J (via Equation [T2l and (if° bs /if F rw ) (via Equation 123]) ; in comparison 
with this, one takes the observationally measured values of the physical mass density wm 
and expansion rate HQ bs = 100 h° bs km s _1 Mpc , and puts them together to form fi^ bs = 
WM/(^° bs ) 2 - One then checks to make sure that Equations I2"5"a.b are satisfied - i.e., that fl^ s 
computed numerically from the model via Equation [23k matches ^M bs fr° m observations. 
Or equivalently (and more conveniently), as we do below, one may adopt some reliable value 
of ^M bs fro m observations, and put it together with Iq and (if ° bs /if FRW ) from the model, 
in order to use Equation [25b to check that fl^™ = 1. If this test is not satisfied to some 
acceptable level of error, then either: the early universe was not spatially flat; ^M bs nas been 
poorly estimated (via measurement errors in lom and/or HQ hs ); the chosen clumping model 
^(t) is not optimal; or there is a problem with the formalism itself (either fundamentally 
or with its simplifying approximations). The goal is therefore to find a \I/(t) that fits the 
supernova data well, while simultaneously achieving fi^j RW — 1 f° r an appropriately specified 
value of tt™ s . 

One last comment about ^m RW 5 i s th & t while having it be equal to unity does represent 
a (primordially) 'flat' universe, this does not necessarily represent a 'critical' universe (i.e., 
steady expansion at a rate smoothly asymptoting to zero), since the future evolution of 
such a universe will strongly depend upon the detailed effects of causal backreaction. These 
issues will be discussed further in Section [71 though we note here that the ever-increasing 
strength of backreaction over time (i.e., I(t) monotonically increasing towards unity for 
large t) would likely require the development of a fully gravitationally-nonlinear treatment 
of causal backreaction in order to determine the true fate of the universe. 




(25a) 
(25b) 
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5.3. More Key Cosmological Observables: g° bs , w° bs , and j® hs 

In conjunction with the fitting of cosmological evolution models to Hubble plots of the 
SNe data, it is also useful to extract a few standard cosmological parameters that characterize 
the da t a in a generalized way (e.g., 'accelerating' versus 'decelerating', etc.). Following 



Visserl (120041 ). we give the definitions of (respectively) the Hubble, deceleration, and jerk 
(or jolt) functions as: H(t) = a/a, q(t) = — (a / a) H \t)~ 2 , and j(t) = (a/a)H(t)~ 3 , where 
a = a 0hs (t 0hs ) is the 'observed' volumetric scale factor (equal to gll 2 , as could be read off 
from the smoothly-inhomogeneous metric, Equation [T3|) as a function of observable time, 
and overdots represent derivatives with respect to t 



Obs 



The limiting values of these functions as z — > 0, t — > to are the well-known parameters 
HQ hs , g° bs , and j° bs , which do not depend explicitly upon the entire cosmic evolutionary 
history, but can be described mathematically in terms of a series of Taylor expansion coeffi- 



ohs \ defined for low-z . As given in iVisser 



(2004 


); 


Riess et al. 


(2004) 



(26) 

Extracting these cosmological parameters requires multiple derivatives of - in our case 
'differentiating' (as described above in Equations [TBlfTU]) the simulated d^ curve for each run. 

The observed and unperturbed Hubble Constants, if^ bs and Hq rw respectively, have 
already been defined in terms of one another and such derivatives via Equation [23j The 
other cosmological parameters can be computed independently of those specific values, by 
taking ratios of the derivatives, as follows: 



and: 



% ° bs = (27a) 

a L,Pert,0 

^o° bs = l(Vo hs -\), (27b) 



Ml 

Jo° bs = -1 + C S + 3(g ° bs ) 2 - • (28) 

a L,Pert,0 



With Equation [2Tb . we have also included a reference to the usual Equation of State 
(EoS) function, w(z), which is usually interpreted (though obviously not in our formalism) 
as a measurement of the pressure properties of the cosmic contents, particularly that of Dark 
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Energy. Note, though, that this w° bs above represents the observed EoS (or whatever effect 
mimics it) of the current epoch of the universe as a whole; it is not the same thing as a 
parameter characterizing the Dark Energy component alone, i.e., Wq. In particular, for a 



matter plus Cosmological Constant 



-1) cosmology with a given Q/ 



'1-Q 



M, 



at 



z = 0, one has w° bs 



In terms of evaluating the third derivative of the luminosity distance function, we have 
deliberately chosen to characterize its behavior in terms of this jerk/jolt parameter, j® hs , 
rather than using alternative formulations. This term in d'£ (i.e., the 0[(z 0bs ) 3 ] term in the 
expansion) is the lowest-order term containing information required to characterize changes 
in the EoS relations of the cosmic contents over time, such as might be the result of an 
evolving, quintessence-like Dark Energy. Parameterizations are therefore often chosen to 
highlight or simplify such an analysis, by defining the Dar k Energy EoS as a function of 
z via parameterization s such as w A (z) = w j 



+ w A 'z (e.g., iRiess et all 120041 ) . or w 



[w A z/ (1 + z)] (e.g.. lKowalski et al.ll2008l ). But we do not do this here, for two important 



reasons. 

First of all, jo is a purely empirical parameter, allowing our analysis to be completely 
agnostic with respect to the physical cause of the apparent acceleration; as our formalism does 
not include any form of Dark Energy, it would be less productive to use parameterizations 
(such as [wq, w a ]) which are optimized to determine the EoS of a Dark Energy which is 
nonexistent in our models. In the language of IVisserl (I2004T ). we are choosing a "retrodictive" 
approach, rather than a "predictive" one, by considering "cosmographic" fits without any 
prior ass umption of Friedmann dynamics. And as has been noted in previous cosmological 
analyses (jCattoen fc Visserll2008t IRiess et al.ll2007l ). the choice of parameterization can have 
a significant impact upon the best-fit results obtained, especially for data with large scatter 
and uncertainties. 

A second useful feature of the jerk parameter, is the property that both Qm = 1 SCDM 
and Cosmological Constant ACDM have j(t) = j(z) = jo = 1 for all time. That is, spatially 
flat ACDM models with w A (z) = —1, containing only pressureless matter ('dust') and vac- 
uum energy, will always have a jerk parameter of unity; a condition that will be true for any 
value of Qa = (1 — ^m), regardless of whether the cosmology is dust-only, dust-dominated, 
A-dominated, or A-only. (And presumably even a slowly-evolving Dark Energy fairly close 
to A, with \dw(z)/dz\ <C 1 and w(z) never too far from —1, would yield jo — 1-) This 
apparently coincidental result - it is not the case when significant radiation is present, for 
example - allows one to conduct a signature test of the entire SCDM/ACDM set of cosmolo- 
gies, for a Dark Energy that is anything close to a Cosmological Constant. Searching for 
deviations from jo = 1 therefore represents an (essentially lowest-order) test of the ACDM 
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(and FLRW) paradigm itself, rather than simply narrowing down the parameters within that 
paradigm. This also represents a falsifiable test of our formalism, since our best-fit simulated 
cosmologies (as will be shown below) generally produce strong deviations of jo from unity. 



5.4. Cosmic Microwave Background Observations and the CMB Acoustic 

Scale 

Testing a theoretical model in Precision Cosmology requires the comparison of model 
predictions against data from several different, independent observational methods, in or- 
der to reduce the effects of large measurement uncertainties, and to disentangle parameter 
degeneracies. One of the most powerful probes of the universe is the Cosmic Microwave 
Background, so we consider here the characteristic angular scale of the CMB acoustic peaks, 
I a, which roughly controls the positioning of the peaks in the CMB power spectrum. 

Our formalism of a smoothly-inhomogeneous universe may be described as minimally- 
disruptive for the CMB, in that few details of the very early universe are altered from the 
standard FRW case: the primordial spectrum of fluctuations is not changed, and neither is 
the spatial flatness (on average) of the universe, the present-day physical density of matter, 
the ratio of baryonic matter to Dark Matter, or just about any parameter affecting the 
physics of the CMB epoch. The only major physical change is that made to the angular 
diameter distance to the last scattering surface; plus some other 'apparent' modifications due 
to differences between the observed (i.e., dressed) cosmological parameters, and the 'true' 
(i.e., bare) parameters. 

The observed acoustic scale, l® hs , is determined by the ratio of two values: the 'standard 
ruler' provided by the CMB sound horizon (r s ), and the angul a r diam eter distance (c?a) to 



the last scattering surface. As discussed in lEfstathiou &: Bondl (119991 ). the exact projection 



of a three-dimensional temperature power spectrum to a two-dimensional angular power 
spectrum is complicated, and depends upon the Doppler peak number m and the shape of the 
primordial power spectrum; but to simplify matters, we may settle here for the approximate 
(flat-space) relationship given in their Equation 21a: l m ~ mirdx/rs, and thus: 

(29) 

Given this simplification, as well as some othe rs (e.g., n eglecting the contribution of the early 



ISW effect to the location of the first peak (jHulll995l ). ignoring radiation in the evolution 



of the scale factor, etc.), our results will therefore not be directly comparable to the precise 
peak positions actually observed in the CMB. However, we will always be consistent in 
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comparisons of this /a parameter between different models, including comparisons of our 
numerical results against the usual SCDM and ACDM models. 

Now, for computing dj±, we note that the angular diameter distance is actually just the 
physical distance to that surface as would be measured at the time of the last scattering. So 
for flat FLRW models, the result would simply be: Oa.flrw = [o(*cmb) ^cmb] = [ao r CMB/(l + 
zcmb)], where tcmb is the coordinate distance traveled by null rays from the CMB to us; 
i.e., the integral of [c/a(t)} from ^cmb to t , where a(t) represents whichever FLRW model 
one chooses. 

Our inhomogeneity-perturbed model, on the other hand, requires several alterations. 
First, we note that the (unperturbed) scale factor, redshift and (coordinate) time of the 
CMB recombination will be changed. Taking z^mb as a gi yen measured parameter, and 
noting that the (pre-clumping) CMB epoch has /(£cmb) — 0, we use Equations [TUfTol to get: 



1 I ^Obs 
1 t Z CMB 

\A + ( V3) 



^FRW _ x ' ^CMB -, / Qn \ 

2 CMB — r, . ,t 7^T ~ 1 ' \ 6K} > 



and hence, since o M d(*cmb) = «cmb = [qmd(^o) / '(I+^cmb)]) anda M D^) = [a (t FRW /t FRW ) 2/3 ]: 



a 



fl CMB - " obs 
1 ' ^r!Mi 



\A + (V3) (31a) 



CMB 



yi + (V3) 13/2 



L Z CMB 

Using this last result above, we can now compute the coordinate distance to the last 
scattering as ?" frw (£cmb)' where this value comes from an integration as defined earlier in 
Equation [16j given whatever clumping evolution function *$f(t) that one is doing a simulation 
of. 

Noting finally that the physical distance to the last scattering, at tcMB, is given by 
[ocmb rFRW (^cMB)]i we can now combine these previous results to get the final expression 
for the CMB angular diameter distance in our inhomogeneity-perturbed formalism: 



VI + ( V3) 

«A,Pert — «0 : obs r \lcMB) > \ 6Z ) 

1 ' Z CMB 



where the arbitrary (though dimensionful) factor ao ultimately cancels out due to the factor 
of 1/ao in r(t), as per Equation [TBb. 

Now by using the purely matter-dominated (MD) expression for the evolution of the 
scale factor in this above calculation - e.g., using a(t FRW ) oc t 2 ^ 3 to obtain Equation [3Tb 
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for t^MB - we have dropped the effects of radiation upon the cosmic expansion rate. We 
have also made this same approximation for the integration in Equation [16] for r FRW (t), 
and for all other integrations and calculations in our formalism. This (greatly simplifying) 
approximation makes little difference for our results, since the dumping-related perturbations 
that we model only become important deep into the MD epoch, by which time radiation 
has become a fairly negligible cosmic component. The exception to this, however, is for 
calculations in which one goes back to very high z - such as in computing expert f° r the 
CMB, or in quoting a total observed age of the universe. But even in those cases, the error 
is only around the ~1% level, certainly accurate enough for a proof-of-principle study in 
which one is comparing different paradigms against one another (as we do in this paper), as 
opposed to a high-precision analysis attempting to extract best-fit cosmological parameters 
from the data. 



Next, to calculate r s for spatially flat FLRW models, lEfstathiou fc Bondl (119991 ) give the 
expression: 

= _c_ 1 1 r cMB d(a/a ) 

r s ,FLRW ^ HoVn - 1 + {[{a/ao) + (1 + R)]l/2 , I ) 

where a eq represents the scale factor at equality between matter and radiation (the latter 
including nearly massless neutrinos), and R = (3pb/4p 7 ) oc (a/ao) is the r atio determining 



the s ound speed of the photon-baryon fluid via c s = [(c/\/3) (1 + R) 1//2 ] ( iHu fc Sugiyama 



19951 ). Using numerical estimates for the present-day densities and properties of radiation 
and neutrinos, they then give the values: (a eqi FLRw/ a o) — (24185 cum) -1 (for three light 
neutrino flavors), and -Rflrw = [30496 (a/a )]; and it becomes a simple matter to perform 
the integration and calculate r Sj FLRW- 

For our model, once again, several modifications must be made. Note fi rst that we have 



alread y made two alterations in Equation l33lfrom the precise formula given in lEfstathiou fc Bond 



( 119991 ). which will be needed for clarity in our following calculations. First, we do not im- 
plicitly normalize all FLRW scale factors (e.g., a eq , ocmb) to ao = 1, but rather normalize 
them explicitly (when necessary) as a/a , etc.; this is necessary because a(t) ^ a/ g rr (t) in 
our formalism when I(t) ^ 0, so that setting a = 1 no longer properly normalizes g rr {to) 
to unity. And second, while their r s is typically called "the sound horizon at decoupling", 
what it actually is, is the sound horizon at decoupling measured today. To convert it to the 
sound horizon size then (for proper comparison to g?a above), we have had to multiply it by 

(aCMB,FLRw/ao,FLRw) = (1 + ^CMB,FLRw) ■ 

In order to generalize these results for our formalism, all of these ratios of scale factors 
have to be modified to account for inhomogeneity-induced perturbations. Once again we 
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have ocmb = [ao\/l + (W3)/(l + ^cmb)] as P er Equation l3Ta. and we can similarly write: 

«eq = r^Oh; V 1 + ( V3) , (34) 



eq 1 + z° bs 

where z° bs is the usual redshift of equality computed from the observed densities according 
to the FLRW formalism. 

Relatedly, the baryon-to-photon ratio as a function of scale factor will need to be ad- 
justed. Noting that R will evolve here in precisely the same way as it does in the FLRW 
formalism when measured as a function of evolving z° hs (rather than a or t), we write: 



i? Pert (t FRW ) = Ro - / / T / \ , 05) 



a_ yj + [I(t FRW )/ 3] 
«o V 1 + (V3) 

to represent the proper volumetric dilution, where the present-day baryon-to-photon ratio 
(assumed fixed by observations) factors in here as Ro = (3pb,o/4p 7 ,o)- We must include the 
factor y/l + [J(t FRW )/3] in the formal definition above to get the correct value of -Rp er t at 
t FRW , though it reduces to unity in these CMB-related calculations because I(t < t init ) = 
in our models (justifiable since 7(£cmb) — I{teq) — in the real universe). 

With these relations and definitions, we can integrate from scratch to find the sound 
horizon using c Si p ert = [(c/v3) (1 + i?p ert ) -1 / 2 ], along with the normal evolution of a radia- 
tion/matter early (FRW) universe. Keeping careful track of all perturbation factors, we get 
the result: 



yi + (W3) 

?"s,Pert — ^ rrFBW . b S X 

'CMB 



V3 H™ 1 + jf» 



— « + V1 + (W3) + a Bo 1/2 a 

a l + ^ q bs a y/l + (I /3) a 

Note that expressing the prefactor of the integral in terms of i? FRW , as opposed to HQ hs , 
automatically eliminates the f^ 1//2 factor (i.e., fi^[ RW is implicitly set to unity, as required 
for our spatially flat dust-only model). 

Inserting the results of Equation [36] (which can be integrated numerically) for r s p crt 
and Equation [32] for expert Equation [291 we finally obtain an expression for Z^ bs = 
(tt rfA,Pert/ r s,Pert) to be calculated for each of our simulation run outputs. To get numer- 
ical results, however, one must assign values to the parameters z^ bs , ^cmbj an d -^o = 
[(3/4) (o;b,o/^7,o)] = [(3/4) (wb,o/wM,o) (1 + 2~ bs )]. (Note that the direct dependence upon 
i/ FRW oc l/t FRW cancels out of the ratio (dA t p ert /r s ); cf. Equations ITSb. [T7b. I3"2"| 131)1) 



Using the best-fit (WMAP-only data) cosmological parameters given in iHinshaw et al. 

,oi 

eq 



(|20Q9j), we have z° bs = 3176, z^mb = 1090.51, and (u hfi /u Mfi ) = K,o/(^cdm,o + Ub.o)] 
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(0.02273/0.13263). We take these as actual measurements of physical observables, which are 
not altered by our smoothly-perturbed universe model; rather, it is the theoretical ("FRW") 
model parameters which must be modified, in order to match these observed values. 

To compare the l® hs values from our models to those of unperturbed FLRW cosmologies, 
we use a flat ACDM cosmology that is best-fit to real SNe data. Further discussion of this 
best-fit optimization step is given below, in Section 15.51 

In order to construct a FRW SCDM model to compare these CMB (and other) calcu- 
lations to, we increase Qm from ~0.27 to 1.0, we multiply Wb,o, wm,o, (1 + ^ bs ) (and hence 
Ro) by (1.0/0.27), and we leave (1 + Zcmb) unchanged. 

With all of these parameters and formulae, we can now calculate the CMB acoustic 
scale l® hs (and all of the other cosmological parameters specified previously) for each of our 
numerically-simulated cosmologies, and compare them to the (homogeneous) ACDM and 
SCDM cases. These results will be collected and discussed soon below, in Section 16.11 

One last remark about the calculations in this subsection, though, is that one must not 
overestimate their precision in estimating CMB observables. There is a great distance in 
time and space between our era and the decoupling epoch, and there will be many effects 
in a universe described by an inhomogeneity-perturbed formalism like ours (e.g., possible 
modifications to lensing, to the ISW effect, etc.), that are not included in these above calcu- 
lations. Such complex effects are useful, since they should ultimately help in distinguishing 
our formalism from FLRW/Dark Energy models; but for now, it must be noted that the 
preceding formulae for observed cosmological parameters - and especially these calculations 
of the CMB acoustic scale - are being derived from a highly simplified and averaged model 
of the universe. 



5.5. Computing x 2 Values and Fit Probabilities for our 
Inhomogeneity-Perturbed Hubble Diagrams 

To quantitatively assess how well our simulated cosmological models (plotted earlier in 
Figures [3]l5]) manage to reproduce the apparent acceleration, we must analyze how well they 
fit a sample of reliab le supernova data. The publicly available "SCP Union" compilation 



( iKowalski et al.l 120081 ) of 307 SNe la (after selection cuts) will serve as the fiducial set of 



supernovae for all of the analyses in this paper, except where stated otherwise. 

The output of each of our numerically-simulated models can be converted into a distance 
modulus function via /z Pert = (m - M) Pert = {5 Log 10 [d LiPert (z Obs )] + 25}, with <i LiPert (z 0bs ) 
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computed according to Equation [171 It can then be turned into a residual distance modulus 
function, AyU Pcrt = A(m — M) Pert , by subtracting off a coasting universe model (as per 
Equation |2~T|) for analysis and for plotting in residual Hubble diagrams. 

Now since this formula for Afip ey - t is auto-normalized to give a zero y- intercept as z — > 0, 
there is no need to specify a particular Hubble Constant in the theoretical model. Real data, 
on the other hand, must be normalized by assigning some value of H to the coasting universe 
model that gets subtracted from the data in order to compute the set of A/z SN = A(m — M) SN 
values. (Generally speaking, there will always be an adjustable factor if° bs needed in order 
to create a comparison between a model, which does not have any spe cific Hubble Constant 



intrinsic to it, versus the observed data, which does.) iRiess et al.l (J2004J) state that the chosen 
value of H (equivalent to the absolute distance scale) is arbitrary, since it only shifts the plot 
of A(m — M) up or down by a constant amount, and that their analysis only depends upon 
differences in magnitude. In practice, however, choosing a poor value of Hq not only makes x 2 
worse for every model, but there is also no guarantee that the relative goodness-of-fit between 
different models will stay the same when this constant offset is changed, particularly in the 
case of data with large uncertainties and scatter. Furthermore, each different theoretical 
model requires a somewhat different vertical offset (i.e., its own individualized if° bs value) 
in order for that specific model, and the dynamical cosmological evolution that it represents, 
to fit a particular SN data set as well as it possibly can. In our analysis, therefore, we 
optimize H^ 8 separately for each simulation run (done simply by adding, and optimizing, 
a constant offset value to A/^sn; for its comparison to that particular AyUp crt curve). Such 
optimization of the observed Hubble Constant is naturally subject to external constraints due 
to other astronomical observations which limit its acceptable range of values; but this is in 
fact advantageous, since the best-fit value of if ° bs that is output for each specific simulation 
run gives us one more independent check upon whether or not that model with those input 
parameters is an acceptable approximation of reality, above and beyond its ability to properly 
fit the SNe data. 

One other challenge in comparing our numerical models to the SNe data is that our 
simulated results for d L are only lists of points for discrete values of z° hs , not continuous 
functions. To evaluate differences between // Pert and /isn for any given z SN , some form 
of interpolation is necessary to make /z Per t continuous. Since our simulations use enough 
points to sample z very finely - with, for example, inter-pixel gaps of Az < 10~ 3 — 10~ 2 for 
<2° bs ;$ 1 — 3, and Az < few x 10~ 4 for z° bs < 0.1 - any decent interpolation scheme should 
give good results. In all of our runs, we compare three different functions for interpolating 
between data points, using polynomials of order 1,2, and 3 in z° hs ; and we find that these 
three interpolation schemes always produce \ 2 values that are identical to one another for 
at least 5 or 6 significant figures. 
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Given such a continuous interpolation function fip ert (z° hs ) for a particular simulation 
run, we can compute the \ 2 value for that theoretical model fit as follows: 



XFit = > , To ! > ( 37 ) 



SN,i 



where //sn,i and <7sN,i for each SN are as given in iKowalski et al.l ( 120081 ) and associated 
SCP data files. (Note that we do not separately fold in additional SNe dispersions due to 
peculiar velocities, lensing, or other intrinsic or systematic effects, as they do in order to 
be conservative in estimating ranges of cosmological parameter values; such additions here 
would simply make it harder to distinguish between the quality of different theoretical fits, 
without providing any new useful information.) 

Once this Xmt value (optimized with respect to -P° bs ) has been calculated for a given 
inhomogeneity-perturbed model, one may compute the likelihood of this fip eit curve by in- 
tegrating the cumulative distribution function for the \ 2 distribution with A DoF degrees of 
freedom (i.e., x^ DoF [X]), as follows: 

P Fi t = 1 - Paw{0 < X < j&J = 1 - / X ^ DoF [X] dX . (38) 

Jo 

This Ppit represents the goodness-of-fit probability that the given theoretical curve, if it 
actually is a correct description of the universe, would give a value of \ 2 as high (or higher) 
than the Xfh that was found. (Although we will see that the Ppit values calculated for 
the models considered here - including best-fit ACDM - are relatively small (~0.3 — 0.4), 
they are not small enough to be a serious concern, since adding in the necessary systematic 
uncertainties not already folded into these cxsn,* values would make all of the fit probabilities 
larger.) 

The relevant number of degrees of freedom here is given by Ad f = (iVgN — iVpi t ), where 
iVsN = 307 for the SCP Union data set, and N Fit is the number of model fitting parameters, 
once a particular type of theoretical model has been chosen. Now first of all, H® hs being 
individually optimized for all theoretical and simulated Hubble curves gives us one fitting 
parameter for every model. Additionally: flat ACDM has Np^ = 2, since only the value of Qa 
(along with H® hs ) is optimizable; and flat SCDM, with no remaining adjustable parameters, 
has iVpit — 1- Alternatively, our ^un, ^Sqr, and $md models each have N Fit = 3, since both 
\l/o and Zi n n can be varied for fitting the data. Thus the number of degrees of freedom for 
evaluating the likelihoods of our inhomogeneity-perturbed models is Ad f = 304; whereas 
the "Concordance" ACDM fit (Q A ~ 0.73 = 1 - M ) has A DoF = 305, and flat SCDM has 
A DoF = 306. 
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One thing that we do not do in this analysis, however , is to compute t he "reduced 



itDoF — (XFit/^DoF)- While it is quite common (e.g., iRiess et al.l ( 12004 ). and many 
other sources) to use Xmt DoF ~ 1 as a proxy for indicating a good model fit to the data, 
it is a poor statistical practice for informally estimating likelihoods. The use of XFit.DoF 
implicitly assumes the approximation, Pi{0 < X < xlitDoF} — Pn OoF {0 < X < Xfh}- But 
XitxIit/^DoF] is a very poor representation of xfv D pbdkt] f° r large -/V DoF , because the former 
distribution is much broader and more gradual: it has a longer right tail, thus leadin g to a 
much larger chance of a Type II error (incorrect acceptance of a false hypothesis, e.g., [Ross 
(119871 )); and the probability distribution increases more slowly as one goes towards lower \ 2 
values, thus also leading to a larger chance of a Type I error (incorrect rejection of a true 
hypothesis). In particular, in the case of the fitting results to be given below for our best 
simulation runs and for the ACDM model, in which the results span xlitDoF ~ 1-02 — 1.05 
for iV DoF ~ 304 — 305, the use of the "reduced" Xmt DoF would cause us to unknowingly mis- 
estimate model-fitting probabilities of Ppit ~ 26% — 38% as probabilities of P F it ~ 30% — 31% 
- a significant loss of comparative information about how well the different models fit the 
data. In this paper, therefore, we will stick to the more accurate statistical practice of simply 
quoting x Fit , N DoF , and the resulting P Fit for each theoretical or simulated model. 



6. OBSERVATIONAL TESTS OF THE FORMALISM 

In order to firmly establish our smoothly-inhomogeneous backreaction formalism as an 
acceptable paradigm for understanding the cosmic evolution, we must ultimately accomplish 
three goals: (1) Explaining the 'already- known' - i.e., reproducing the most important obser- 
vational results that have formerly been interpreted as signs of Dark Energy; (2) Providing 
'falsifiability' for our formalism, by establishing new predictions that can clearly distinguish 
our model from the conventional ACDM Concordance paradigm; and lastly: (3) Convinc- 
ingly argue plausibility or 'naturalness', in the sense of showing that our model does not 
suffer from similar fine-tuning problems as the Dark Energy approach. For this third point, 
in particular, while the linking of the onset of the apparent acceleration with the beginning 
of widespread structure formation clearly removes the Coincidence Problem of Cosmological 
Constant Dark Energy (recall Section [TJ, it still remains to be shown why an alternative, 
non-A explanation of the cosmic acceleration, such as ours - even if true - should hap- 
pen 'fortuitously' to look so much like the action of a trivially simple (if aesthetically and 
coincident ally unpleasant) Cosmological Constant. 

The first two of these three tasks will be discussed below, in Subsections 16. II and I6.2[ re- 
spectively. The third task - addressing the basic necessity, itself, of employing an alternative 
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cosmology in the face of current cosmological data sets that appear to be fairly consistent 
with the ACDM Concordance Model - will be reserved for a future treatment, both because 
of space limitations here, and due to the fact that this is a generic question for all alter- 
native cosmologies, not a unique issue for the causal backreaction formalism introduced in 
this paper. We will therefore present our detailed discussion on that topic elsewhere; but in 
any case, it is likely that all of these above issues cannot be resolved to a satisfactory de- 
gree without much more comprehensive and precise cosmological data, along with continued 
input from the broader expertise of the entire cosmological community. 

6.1. Supernova Fits and Cosmological Parameters from our Numerical 
Simulations: Evaluating the New Concordance 

Using the fit probabilities, Pfh, obtainable from each of the simulation runs presented 
earlier in Section HJ in conjunction with the cosmological parameters calculated from each 
run according to the formulae in Section 0, we now have a quantitative context for judging 
our causal backreaction paradigm and the clumping evolution functions that have been 
modeled for this paper. Given those results, we have specified an informally-chosen set of 
'best' runs for further detailed discussion here. Of the sixty cosmological models plotted 
previously in Figures [3J1S1 twelve of them (six ^ Lin runs and six $ MD runs) appear to this au- 
thor as being 'very good' at replicating the apparently accelerating behavior of the universe, 
while also having fairly good cosmological parameters. 

Now, given that these so-called "best runs" are chosen simply from a discrete set of 60 
runs performed overall, it would seem likely that a truly optimized search over the (^>q, z^) 
parameter space might find model runs that do an even better job at fitting the SNe data. 
Such an effort is not really called for here, though, given the theoretical uncertainties of our 
formalism - particularly the neglect of both gravitational and recursive nonlinearities, as 
well as the ad-hoc nature of our chosen clumping functions, themselves - but such an 
optimization can indeed be done, producing many fits with values of Xfh that are just as low 
(and often slightly lower) than that of the best-fit flat ACDM model. This is most naturally 
accomplished by choosing parameters to weaken the $md runs, since strengthening the 
'accelerative' effects of the ^tia or ^Sqr runs requires one to implausibly resort to 2> 25 
or > 1- Exploring the range of models with f MD clumping evolution, there turns out 
to be an extensive 'trench' in parameter space extending at least from (\& , z mit) = (0.89, 5) 
to (0.765, 14) (and likely beyond) for which Xfh ~ 311.7 always remains true, allowing one 
to fit the SNe data with these models slightly better than one can with ACDM over a wide 
range of model input parameters, thus simultaneously providing good fits and opening up an 
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extended range of output cosmological parameters to more flexibly match the model results 
to other external observational constraints. 

For illustration purposes, we choose one of these input-parameter-optimized runs as 
an additional example model for comparing to the SNe data and observed cosmological 
parameters. All of the $md runs in this phase-space trench are virtually degenerate in their 
xlit values (there is an extremely weak gradient in Ppit among them), so as a somewhat 
loosely defined "semi-optimized" run we select the (ty , z init ) = (0.768, 14) case, which has 
the lowest xpit °f an y *md run that we found among those tested (in a non-exhaustive search 
through the parameter space), while also having sufficiently different parameters from our 
discrete sample of 20 $md runs for it to be of further cosmological interest. We thus add this 
case to our twelve other "best runs" discussed above, giving us a set of thirteen highlighted 
runs, in total, for use in the more in-depth discussion of the results of our simulations which 
follows now. 

Residual Hubble diagrams of these thirteen chosen runs are shown in Figure El plotted 
along with the (now SNe-best-fit) SCDM and Concordance ACDM models. Shown against 
these curves are the SCP Union SNe data, with the 307 SNe averaged here into 40 bins, 
equally spaced in Log[l + z ]; binning is necessary for clear visualization of this data set, 
given not only the large number of SNe data points, but also the very high scatter for those 
data, and the large error bars for each individual supernova magnitude. 

It is obvious by inspection that these thirteen runs represent cosmological models that 
produce good Hubble curves, being visually almost indistinguishable from one another (and 
from Concordance ACDM) in the SNe-data-rich region of z° hs ~ 0.1 — I. This is especially 
evident when plotted on this y-axis scale; a scale that clearly shows the large separation 
between all of these models from SCDM, as well as the large scatter and error bars of the 
SNe data themselves - even after binning and averaging - when contrasted with the tight 
overlap between all of these models and the Concordance ACDM cosmology itself. 

For a more quantitative analysis, the comprehensive output data for these thirteen runs, 
along with corresponding output parameters from the best-fit ACDM and SCDM runs, are 
given in Table [TJ 

First, considering the goodness-of-fit of the models to the SNe data, we see that the 
results for this discrete selection of inhomogeneity-perturbed models are nearly as good as 
that for ACDM in terms of Xmv anc ^ they are comparable in terms of Pfu- Our thirteen 
'best' models range from Xfh Pert ~ 311.7 — 319.0, with the "semi-optimized" run (and others 
near it in the 'preferred $md parameter space trench' discussed above) having xfit Pert < 
X Fit acdm = 311.9; and the Ppit values for most of these thirteen models are almost as good 
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as the Concordance ACDM result here, with ten of them having PFit,Pert > 0.3, and five 
of those going as high as -PFit.Pcrt ~ 0.36 — 0.37 (as compared to Pfh,ACdm = 0.38 for the 
fully-optimized case of Q\ = 0.713, H = 69.96 km s _1 Mpc _1 ). In short, it seems fair 
to conclude that the smoothly-inhomogeneous formalism presented in this paper - even at 
this 'toy model' stage, given all of the simplifications discussed above - is already able to 
successfully reproduce the apparent acceleration of the universe essentially as well as the 
more standard ACDM paradigm can do. 

Next, to check on the general validity of these models all the way to t = t , we must 
consider the diagnostic parameter Iq = I (to), which varies from Iq ~ 0.52 — 0.72 for these 
thirteen best runs, and stays bounded at Io < 0.66 for eleven of them. As discussed in 
Section 13.21 where I(t) was formally constructed, this function (representing the causally- 
integrated influence of inhomogeneities) is only valid in the linear context of Newtonian 
gravitational terms, for which the metric perturbations from different masses can be simply 
summed together. To the extent that I(t) approaches unity, the linearized-gravity approach 
of our formalism breaks down, and full general relativity becomes necessary for describing 
the combined effect of all perturbations, in total. 

Now, seeing Iq take on values roughly midway between and 1 for these runs is in fact 
a good sign for trusting the general predictions of these models. If Iq (and thus d[I(t)]/dt 
as well) were too small, then the total perturbative effect of these inhomogeneities on the 
average cosmic gravitational potential would be too weak to cause much of an observable 
effect at all (recall Equations @Ji,b and the related discussion) - certainly not a reaction strong 
enough to explain the apparent acceleration. For example, the very weak FRW-perturbing 
behaviors of the \l/s qr models (cf. Figures HJ), which have Iq < 0.2, is symptomatic of their 
limited accelerative effects (i.e., w° bs > —0.5). On the other hand, too large a value - such as 
I ~ 0.93, from the strongest of all of our models (the $md case with (\& , z init ) = (1.0, 25)) 
- is so large that the detailed quantitative results of such a model cannot really be trusted. 

The fact that the perturbative effects from our thirteen best-fitting model cosmologies 
do not completely overshoot the linearity approximation, and thus remain fairly trustworthy 
in their output results, is not completely a matter of luck: one key feature of the cosmic 
acceleration (and the source of the Coincidence Problem) is that it seems to have just re- 
cently begun, being due (in the backreaction paradigm) to the same structure formation 
that recently (in cosmological terms) has created us, as well. But this 'luck' will not hold 
out in the long term, since I(t) should continue to grow in the future for many billions of 
years, until I(t 3> t ) eventually becomes so large that the evolution of the real universe 
itself (and not just our phenomenological representation of it) should cease to be describable 
by any simply-perturbed FRW/FLRW model. The dynamical behavior of such a universe 
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would likely become quite extreme and difficult to predict, and we will speculate upon the 
possible cosmic futures in Section 

Moving on to the observable cosmological parameters for each run, we first note the 
changes to z° hs , so that it is now no longer precisely equal to z FRW . Table [1] gives results 
corresponding to z FRW = 1 (i.e., [ao/a(t)] = 2) as a sample epoch, demonstrating that while 
the differences are not huge (e.g., z° hs ~ 11 — 14% larger than ^ FRW for these thirteen runs), 
they are also non-negligible, implying that the backreaction effects of inhomogeneities have 
a real ability to alter the relationships between observed cosmological time, volume, and 
redshift, which must be taken into account for any precision understanding of the evolution 
of astrophysical objects. 

The effects upon the Hubble Constant, on the other hand, are much larger in magnitude, 
as we have both expected and required, in order to create an alternative concordance. First, 
considering the individually-optimized values of the observed Hubble Constant for each run, 
we see that our thirteen best runs (along with the flat ACDM model) have a range of 
H° hs ~ 68.8 - 71.8, which compares well w ith the result of H n = 70.1 ± 1.3 km s x Mpc 



from combined (WMAP+BAO+SN) data (IHinshaw et al.ll2009l ). This is an encouraging 



cross-check of these runs, in contrast to the utter failure of flu = 1 SCDM (with optimized 
H® hs = 61.35) to properly fit either the SNe data or the Hubble Constant; similar to the 
tendency of our model simulation runs with poor SNe Ppit values (not presented here) to 
also have observationally discrepant H® hs values. 

Next, considering the unperturbed/ "FRW" Hubble Constant values for these thirteen 
best runs, we get results of if FRW ~ 34.2 — 42.8, with six of the very best values lying within 
Hq RW ~ 36.1—41.6. These results are entirely in accord with our discussion from Section l5TTj 
in which we noted how the findings of other authors have indicated that a variety of other 
cosmological measurements are seemingly concordant with — 1 models that have a low 
Hubble Constant, only to be stymied by direct observations apparently indicating H ~ 70. 
But here we see that the ability of our inhomogeneity-perturbed formalism to achieve a value 
of HQ hs in the ~70's, despite having a 'true' Hubble Constant - i.e., the one that actually 
mattered during the pre-structure-formation cosmic epoch - of if FRW in the ~30's-40's, 
completely removes this contradiction. Furthermore, it also helps in the establishment of an 
alternative, matter-only concordance in several other ways. 

One important, related aspect of this low value of if FRW is that we have now solved the 
Age Problem in consequence, as described above in Section 15. 1[ without needing a period 
of recent 'real' acceleration that would nominally be provided by Dark Energy. Our results 
for these thirteen runs are t® bs — 13.2 — 14.5 GYr (and t® bs ~ 13.4 — 14.2 GYr for the nine 
best fits), which compares very favorably to the SN-best-fit flat ACDM value of 13.64 GYr 
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(noting again that the time estimates here are not exact, since none of these models include 
radiation; though they are all mutually consistent for = 0) . Now, a cosmic age difference 
of a ~billion years either way is certainly not negligible, and careful observations of objects 
like globular clusters, etc., could be expected to someday strongly discriminate between 
ACDM and our formalism (as well to as help optimize parameters within our formalism, for 
different clumping evolution models). But given the typical state of cosmic age lower-limits, 
taken from the me asured ages of a strop hysical objects - for example, like that from a review 
on Cosmic Age in ISpergel et all ( 120031 ) , where the strongest lower limit that they quoted 
was 12.7 ± 0.7 GYr, from observations of cooling White Dwarfs - it seems reasonable to 
conclude that the values given here for t^ bs from our formalism, even at this 'toy model' 
stage, are good enough to have effectively solved the Age Problem possessed by SCDM; 
especially considering that to,scDM ; i n contrast, is as low as ~10.6 GYr, as seen in Tabled) 

Given if FRW (along with iJ^ bs and Iq) for each run, we can also now address the 
question of the spatial flatness of the universe, by computing fly^ via Equation [25b . This 
first requires us to choose some value of fl^ s from cosmological observations. Now, this 
value need not be equal to (1 — 0.713) = 0.287, which one might infer from the SN-best-fit 
flat ACDM model above, since it would be better to use a value that comes from a more 
comprehensive combination of different data sets; and a value of ^M bs — 0.287 would in fact 
be rather high, given most recent estimates of ^M bs - 



The actual amount of matter as a fraction of the apparen t (i.e., observational ) closure 
density is still fairly difficult to pin down with great precision. ISpergel et all (120071 ) demon- 
strated in their Tables 5 and 6 that which best-fit value of ^M bs they obtained depended upon 
which external data set (if any) that they chose to c ombine the (3 r d -Year ) WMAP data set 
with, with variations of Qm ~ 0.226 — 0.299 possible. iDunkley et al.l ( 120091 ) noted the tension 
in the preferred (high versus low) values of the matter density observed by SDSS (flu = 
0.265±0.3), and 2dFGRS {Q M = 0.236±0.02). And while they described how the uncertainty 
in the matter density had dropped with each new analysis of the growing WMAP data set, it 
is also true that there has been some oscilla tion in those estimat ed values, with the 1 st - Year 
WMAP analysis giving fi M = 0.27 ± 0.04 (feennett et al.ll2003l ). the 3- Year WMAP mean 
dropping to Qu = 0.241 ±0.034, and the 5- Yea r WMAP mean rising back up (despite a max 
likelihood of only 0.249) to fi M = 0.258 ± 0.03 ( faunklev et~al"]l2009[ Table 2) Furthermore, 



that 5- Year WMAP-only result actually increases to VIm — 1 — 0.721 = 0.279 (IHinshaw et al. 



2009 



, Table 6) when complementary data sets are included (WMAP+BAO+SN). Similarly, 
Larson et al.l (120 111 ) give (Q c + Q b ) — 0.274 for their (updated) 7- Year WMAP values. 



For the calculations of flYF^ in this paper, we choose to adopt what seems to be a 



generally reasonable value for our search for a new concordance: = 0.27. What must 
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be realized in these calculations of the 'true' (pre-perturbation) spatial curvature fiM RW 5 
however, is that although we have adopted one specific value of ^M bs f° r quoting values of 
f2^ RW , we can only trust our estimates of spatial flatness (i.e., ^m RW ~* 1) as being accurate 
to within the observational uncertainties in ^ly^ s ] for example, to within a (very approximate) 
range of, say, ~(l/0.27) • [±0.1] ~ ±[0.3 - 0.4], or so. 

Looking at the fi^ RW results for our thirteen best runs in Table [TJ we can clearly say 
that to within these uncertainties, our smoothly-inhomogeneous cosmological formalism has 
indeed achieved spatial flatness; and that this has been done without requiring the incorpo- 
ration of any "missing mass" or "missing energy" (beyond the usual Cold/Hot Dark Matter), 
such as Dark Energy, to bring f2xot up to unity. 

Specifically, we see that these V&Lin runs yield fi^ RW ~ 0.90 - 1.04, and that the * M d 
runs yield ^m RW ~ 0-99 — 1.64 (narrowing down further to fi^P ~ 0.99 — 1.41, if we drop 
the worst of these seven $md cases). Thus the ^Lin runs are especially good at achieving 
flatness. And the $md runs, even with their stronger effects, are not too far off either: if 
one (reasonably, as we have seen) chooses a lower value of the observed matter density - 
say, = 0.24, instead of 0.27 - then this drops the range down to ^ RW ~ 0.88 - 1.25 
for the six best $md runs, making them even more consistent with f2™ W = ^Tot = 1- 
Even given the observational uncertainties, as well as the many theoretical simplifications 
of our formalism, it is a notable step towards the achievement of an alternative concordance 
that such a trial-and-error set of sixty simulated cosmological models (with astrophysically- 
motivated parameters) has been found to include about a dozen runs that succeed in bringing 
^frw f rom ^0.3 up to a reasonably close range around 1, using matter alone. 

Next, we consider the quantitative amount of apparent acceleration - i.e., the degree to 
which g° bs < 0, w° bs < (—1/3) - that our models have produced. Table CD shows that each 
of these thirteen best runs have produced a strong amount of 'acceleration', with 
(-0.71) - (-0.82) [i.e., q$ hs ~ (-0.56) - (-0.73)] for the ^ Li n runs; and w$ bs ~ (-0.75) - 
(-1.0) [i.e., q$ hs ~ (-0.62) - (-1.0)] for the ^ M d runs. 

Now, a flat ACDM universe with a given value of (Cosmological Constant) Q\ will have 
w obs = _^ a _ Taking the v alues of Q A = . 742 ± 0.03 (WMAP-only) and Q A = 0.721 ±0.015 



(WMAP+BAO+SN) from lHinshaw et al.l (120091 ) as a guide, it seems like we would look for 



our models to reproduce values in the range of, say, w° bs ~ (—0.70) — (—0.77), in order to 
look like Concordance ACDM. As we see, four of these ^Lin runs, and one of these $md runs, 
do in fact fall within this range, with the others producing somewhat stronger 'acceleration' 
effects. 



While this is already a fairly good success rate for our models in reproducing the appar- 
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ent acceleration, it is important to note that the precise value of WQ hs which is 'observed' is ac- 
tually just the result of some particular best-fit procedure, given some assumed cosmological 
model; and output parameters from best -fits depend non-trivially upon the fitting function 



used. For example, iKowalski et al.l (120081 ) quote a mean value (from SN+BAO+CMB data) 
of Wq = —0.969 for a const ant-w A analysis, which in conjunction with their fl\ = 0.713, 
turns into w® hs = (0.713) • (—0.969) = —0.691. Alternatively, in their varying Dark Energy 
EoS analysis using w A (z) = Wq + [w£z/(l + z)], they get (Rubin, D. 2008, private com- 
munication) a mean value of = —1.13 (with u> A = 0.73), which in conjunction with a 
re-fit value of = 0.718, yields WQ hs = —0.811; this is clearly a substantial change from 
WQ hs = —0.691, just given a change in the fitting assumptions. The conclusion here is that 
it is not necessary to precisely match the amount of 'acceleration' (i.e., WQ hs , w° hs (z)) found 
with the best-fit Dark Energy models, in order to say that one has reproduced the apparent 
cosmic acceleration. Rather, all that one has to do is to produce an apparent acceleration 
that is very roughly in the correct range of the old Concordance u>° bs , while simultaneously 
producing a fit to the SNe data that is essentially as good or better than ACDM. As seen 
above, we have achieved this by obtaining ~5 — 10 runs here with fits of a quality comparable 
to that of ACDM, with nearly all of them having an apparent acceleration within the fairly 
good range of w° bs ~ (—0.7) — (—0.9). 

One further way of characterizing the cosmic acceleration, in particular its 'sudden' 
onset, is to descr ibe the universe as having recently experienced a strong, positive "jerk" 



(IRiess et al.ll2004l ). Given our definition of the jerk/jolt parameter in Equation we find 
our six runs in Table [T] to have j° bs ~ 2.5 — 3.5, while the seven \1/md runs have 
jo hs ~ 2.6 — 5.5. Recalling from Section 15.31 that spatially flat cases of both SCDM and 
Cosmological Constant ACDM models (for any value of I^a = 1 — ^m) always have j° bs = 1, 
it thus appears that these universally high values of j° bs - universal, that is, for those of our 
models also capable of fitting the SNe data and w° bs properly - is the result representing the 
most discriminating test that we have found so far which could potentially distinguish our 
formalism from traditional Concordance Models with a Cosmological Constant (or anything 
close to it) that are capable of generating the apparent acceleration. Given that the quantity 
j° bs - or its corresponding alternative in 'dynamic' parameterizations, w A - comes from the 
third-derivative term in the expansion for in z° hs (cf. Equation l2~6"j) . its value is not yet 
observationally well constrained. In that sense, we will regard the output values of from 
our runs as a prediction of our formalism, rather than as an attempt to match known results; 
a matter which we will explore further in Section 16.2.21 below, where we consider (in brief) 
the observational situation for j° bs as it stands. 

Unfortunately, though measurements of j° bs > 1 would conclusively rule out ACDM, 
we cannot conclude with certainty, based only upon these above results, that high jerk 
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parameters are a robust feature of our formalism. Too many theoretical uncertainties exist 
within our models to make this an 'iron-clad' prediction of our calculations just yet. A 
number of different effects should all act to moderate the strength of the apparent acceleration 
as z — )■ 0, t — >■ t - thus weakening the late-time "cosmic jerk" - both within our model 
calculations and in the real universe. 

First, there is the issue of "recursive nonlinear ities" discussed earlier in Section 1372] the 
fact the pre-existing inhomogeneities serve to slow down the causal updating which brings 
in later information about new inhomogeneities, thus softening the apparent acceleration 
effect by some undetermined degree. As calculating these effects properly would take an 
algorithmically nonlinear simulation program, we cannot quantitatively estimate them here; 
all that we can say is that they would provide some significant amount of damping of the 
apparent acceleration effects due to causal backreaction, especially at later times. 

Second, even if we could safely rely upon this and other simplifications in our formal- 
ism, and even if the functional forms of the and/or $md clumping evolution models 
(cf. Equations l2"2"k.b) generally do make sense for the long-term growth of structure in 
the universe, there is still no guarantee that those cluster evolution functions would remain 
meaningful all the way to t — > to- As structure forms, the large-scale collapse of material 
and feedback from star formation act to shock heat cosmic baryons to millions of degrees, 
particularly at late times (z° bs < 3), thus inhibiting clumping^ by keeping and /or sendin g 



a significant portion of the baryons into the superheated IGM (e.g.. lCen Sz Ostrikerll2006l ). 
Now, our models can account for some of this effect, simply by using a lower value of 
but if this shock-heating is strong enough to really soften the evolution of growth of cosmic 
structure (thus significantly altering the functional form of ^(t), itself), then this would be 
a moderating effect that is not taken account of by any of the (highly simplified) clumping 
evolution functions that we use in this paper. 

Any of these effects which weaken the apparent acceleration as z — > are generally bad 
for our models, since this could make at least some of them too weak to produce a new 
cosmic concordance; alternatively, they are generally good for our f MD models, by reducing 
their computed values of (if^ bs /i^ RW ), f^ RW , and WQ hs . In any case, however, such effects 



9 On the other hand, given that simulations 
decreasing as z — > due to continued relaxation, this means that at least some virialization does continue 
all the way to the current time, even for structures which largely are 'already clumped' - thus implying 
a continuing "backreaction" of the kind which we phenomenologically model here. But whether or not 
this continuing level of virialization would be strong enough to substantially counter the slowing of such 
backreaction at late times due to shock heating, it would in any case represent yet more complexity in 
determining the most physically realistic functional form to use for ^(t). 



Ho et al.ll200ffl show that galaxy cluster ellipticities are still 
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would likely lead to a reduction in our calculated output values of j® hs . The magnitude of 
such a change is not certain, since the primary effect of such alterations might simply be 
to change which choices of model input parameters (i.e., which (\I/ , z init ) and/or which type 
of \I/(t) function) manages to fit the cosmological data best, without changing the general 
values of j° bs (or of other output parameters) that tend to emerge from those models that 
do achieve good SNe fits. But all we can say for certain here, examining Tabled] once more, 
is that for the seven runs 'most acceptable' on every front - i.e., those runs with Xfh < 315, 
t° bs ~ 13.4 - 14.0 GYr, fi™ w ~ 0.97 - 1.3, and w Obs ~ (-0.75) - (-0.87) - that we have 
the range of j^ bs ~ 2.6 — 3.8 for the jerk parameter. It therefore seems reasonable, at this 
stage, to propose that such a range of j® hs may be a generic prediction from our formalism 
for models which succeed in achieving a good alternative concordance - a conclusion subject, 
of course, to further theoretical development of this paradigm. 

Lastly in our evaluation of these thirteen best model runs, we consider the acoustic scale 
of the CMB peaks, /° bs , as derived above in Section I5\4l As noted there, this parameter is the 
one incorporating the most far-ranging assumptions and simplifications, and the one most 
prone to error due to the great look-back time and distance to the last scattering surface. 
The theoretical uncertainties in our results for Z° bs are therefore expected to be large, though 
we cannot precisely quantify them here. 

Table [1] shows that these six best $Lm runs have Z£ bs ~ 284.2-291.3, while the six best 
$md runs (dropping the worst of these $md SNe fits) have Z° bs ~ 270.5 — 288.7. Narrowing 
it down even further, the six runs most acceptable on every front (i.e., the seven 'most 
acceptable' runs referred to previously, now minus the "Semi-Optimized" ^md run) have 
/g bs ~ 277.5 - 288.7. Those values succeed very well in bracketing the SN-best-fit ACDM 
value (for no radiation) of /aacdm = 285.4, while simultaneously producing good values for 
all of the other cosmological parameters computed in this paper, as well. 

Such results, while serving as a fairly good match of our models to the CMB acoustic 
scale, may not quite match the precision of the actual ±~0.8 — 0.9 measureme nt uncertainties 



quote d in observational CMB parameter estimation results for Z^ bs (e.g., iHinshaw et al. 



20091 ). Nevertheless, it is significant to note how much better these estimates are than any 
fit achievable with a I ow-^m open CDM model without Dark Energy. For example, using the 
rough rule-of-thumb ( Turner 1999 ) of lf hs ~ 200 ^To/ 2 f° r the ^ TS ^ CMB peak, shifting from 



a flat SCDM universe (Om = 1) to an oCDM universe (Qm ~ 0.3, = 0) increases /!p bs by 
~200[(l/-\/oT3) — 1] ~ 165 = A/° bs . Keeping ^M bs ~ 0-3 without some mechanism to bring 
about spatial flatness - such as Dark Energy, or our smoothly-inhomogeneous perturbation 
formalism - therefore leads to errors in Z^ bs of over a hundred; whereas insisting on flat 
SCDM mostly saves Z° bs (see Table [[]), but at the cost of ruining the fit to the SNe data 
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(Ppit — 3.4 x 10 -22 ). Our formalism, on the other hand, preserves the good SNe fit (and 
good values for all of the other cosmological parameters calculated here), while producing 
deviations from the Concordance ACDM model of only Al® hs < 15 for twelve of these 
thirteen best runs (and to within Al® hs ~ 0.2 — 7.9 for the six very best runs), all while 
maintaining Q\ = 0. 

To sum up the results of this subsection: from the basic set of sixty simulation runs 
done with our inhomogeneity-perturbed formalism using astrophysically-motivated input 
parameters, we have obtained twelve runs that produce quantitatively good fits to residual 
Hubble diagrams of the SCP Union supernovae, many of them very close in Xfh to the 
best-fit done with flat ACDM. In addition to these, we include a thirteenth run, taken as 
one example of a long trench in (\& , z init )-space containing many f MD runs which produce 
nearly equal values of Xmti an lower than that achievable with any flat ACDM model. These 
thirteen-plus causal backreaction models have therefore successfully reproduced the signs of 
apparent cosmic acceleration that are seen in Type la SNe data sets. Furthermore, as it is not 
surprising that good fits to the SNe data also imply good values for other parameters, we see 
that six of these runs - specifically, ^Lin with (^o,z init ) = (1.0,25), (1.0, 15), and (0.96,25), 
and f M D with z ini t) = (0.78,10), (0.85,5), and (0.92,5)0- not only provide excellent 
SNe fits, but also reproduce several other cosmological parameters to within acceptable 
'Concordance-level' bounds, including tQ bs , fi^ RW , u^ bs , and /° bs . And beyond just achieving 
an agreement for already-known observable parameters with these models, we also find them 
to yield j® hs ~ 2.6 — 3.8, indicating that observations of a jerk parameter value significantly 
greater than unity would as of now appear (pending further theoretical refinements of these 
calculations) to serve as a potentially powerful way of distinguishing our formalism from 
Cosmological Constant ACDM, as we will examine further next. 

We note once again that all of this is achieved in our models without the incorpora- 
tion of any Dark- Energy- type (i.e., negative pressure) species, but is instead done simply 
with a spatially flat, matter-dominated universe that is perturbed by causally-propagating 
information about self-stabilizing inhomogeneities. 



10 Recalling once more that larger Zi n a and mid-range V&o are the most 'astrophysically sensible' values for 
these parameters, generally speaking. 
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6.2. Distinguishing the Smoothly-Inhomogeneous Formalism from 

Concordance ACDM 

A candidate paradigm for an alternative Cosmic Concordance, such as our formalism, 
must simultaneously satisfy a variety of fairly independent observational constraints, includ- 
ing the parameters considered above plus several other types of observations. 

Assuming that such a goal can be fully achieved, it becomes rather difficult then to find 
distinctions which would be capable of demonstrating a statistically strong preference for 
the alternative concordance over one achievable with Dark Energy (DE); especially so in the 
case of data with high scatter and uncertainties, and a paradigm as malleable as DE with an 
optimizable equation of state. Nevertheless, we will discuss a few potential methods here to 
distinguish our formalism from Dark Energy in general, where possible; and more feasibly, 
from Cosmological Constant DE in particular, which despite being aesthetically questionable 
and more constrained (and thus easier to falsify) than evolving DE, still appears consistent 
with a wide range of observations. 



6.2.1. Direct and Indirect Effects of Inhomogeneities: Observational Anisotropics and 

Other Cosmological Signatures 

One obvious departure of our formalism from any non-clustering version of DE is the 
importance that our model places upon the combined influence of many localized inhomo- 
geneities. Though our calculations are done using a "smoothly" inhomogeneous ansatz in 
which spatial variations are averaged away, in the real universe these effects will not perfectly 
average to smoothness, and thus it becomes more interesting in our paradigm to look for 
evidence of unexpectedly large anisotropies on a variety of scales. 

Though there seems to be no substantial rejection yet of the large-scale homogeneity 
assumed by the Cosmological Principle and typical FLRW cosmologies, some intriguing 
observation al results have t urned up using various methodologies. A prelimina ry study by 



this a uthor flBochnerll2007bl ) of anisotropies in the Riess gold04 SNe compilation ( IRiess et al. 
2004 ) showed some marginal, positive signs of the existence of real anisotropies; and studies 
using SNe, galaxy clusters, etc., to map anisotropies (or find e n hanced variances) in th e 



Hubble Flow by other resea rc hers - e.g. , Kolatt fc Lahav J2001 ); McC 



Schwarz & Weinhorsi 


t f 


2007) 


(2009 


); 


Colin et al. 


( 


2O10h - 



ure & Dverl (120071): 



varying statistical significance. 

Alternatively, considering inhomogeneities seen via the CMB, there have been stud- 
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ies co ncerning; the possibi lity of significant an i sotropies ( Hanse n et al.l l2004t iBernui et al. 
2007 ). non-G aussianities (|McEwen et al. 2005 ; Tojeiro et al. 2006 ). a nd other unexpected 



features fe.g.. [Larson fc Wandeltl 12004 ) present in the WM AP data ( Hinshaw et al. 



Dunkley et al. 2009h . such as the low CMB quadrupole (e.g., Jarosik et al. 201(1) 



Axis of Evil flLand & Magueiiol 120051). and the W MAP Cold Spot (e.g. 



2007 



Cruz et al 



the CMB 
2007h. 



Also, some CMB glitches (e.g JHunt fc Sarkaii2007h . though wiped out by binning (IHinshaw et al 
20071 ). may perhaps be signs of real inhomogeneity effec ts that should not b e averaged away. 
(Though one must not ignore counterarguments (e.g.. iBennett et al.l 120101 ) that anomalies 
such as these mentioned here might in large part be due to a posteriori selection effects and 
technological limitations in the observations.) 

Still more directly, one may investigate major cosmological structures themselves in 



detai l, such as the Shapley Concent ration and /or Gr eat Attractor regions (e.g.. lBardeUi et al. 



1994 ) and the Sloan Great Wall (IGott et al.l 120051 ). and use them as part of the effort to 



Colin et al. 


2010). Such studies 


"Dark Flow" ( 


Kashlinskv et al. 



20081 . |2010| ) which may be due to a tilt across our entire current observational horizon (if 
not due to some causal- inhomoge neity-driven mechani sm) , as well as other evidence of bulk 
flows out to 100/i _1 Mpc or more (IFeldman et al.ll2010l ). 



Searches like these for cosmic inhomogeneities and/or anisotropies are already fully 
underway through the work of many researchers, as a way of determining the fine details 
of the cosmic evolution. Our work here merely serves to add to the impetus behind such 
investigations, given our proposal that inhomogeneities may very well determine the average 
observed cosmological parameters themselves, not just the deviations from the averaged, 
so-called best-fits. 

Even without direct evidence of inhomogeneities, though, the method of apparent cosmic 
acceleration that we propose here would undoubtedly impose changes upon several other 
cosmological measurements, changes which we have not explicitly estimated in this paper. 
Particularly interesting are observations which are not explained by (or which directly conflict 
with) the ACDM paradigm. For example, the fact that causal backreaction does not occur in 
completely smoothly-distributed fashion, but is instead concentrated mostly near virializing 
masses, may lead to interesting clustering-induced feedback behaviors affecting the formation 
of stars, galaxies and galaxy cluster s. Conceivably, th is may have some useful application 



to issues such as galaxy downsizing ( jCowie et al. 



the possible dearth of dwarf satellite galaxies ( IPrimackl 120031 ); but such connections are 



1996k t he cus py CDM halo problem, and 



speculative, and require a detailed quantitative analysis to see if causal backreaction truly 
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succeeds (or significantly helps) in explaining these issues. 



Other examples of interesting processes or o bservables whic h might be affected by causal 
backreaction include the Late ISW effect (e.g., iHo et al.l 120081 ); the shape parameter, V = 



Om^ (jTurnerlll999l ); observations of peculiar velocities constraining cts^m 6 (lLahav fc Liddle 
20081 ); large-scale velocity flows constraining = Q®f/b, with b as the l inear bias facto r for 
galaxies versus the overall matter distribution ( jFukugita fc Hoganl|2000l ; lFukugitall200ll ); ob- 
servations of weak lensing due to large-scale structure, constraining ~<r 8 f2^ 53 ~ ' 64 , depending 
upon the angular scales used ( Fu et al. 20081 ); Baryon Acoustic Oscillations, which constrain 
the combination A = {v^m [Ho/ H^*)} 1 / 3 ^*/ z*] 2 / 3 }, wit h r* being the dime nsionless co- 
moving distance to sampled data having typica l redshift z* (IMarassi et al. I2010T) : indications 
of ac celeration from the Alcock-Paczyhski test f Alcock fc Paczyhski 19791 : Marinoni fc Buzzi 
20101 ); and observations of the growth and evolution of lar ge galaxy clusters, analyzed in con- 



junct ion with several other types of measurements (e.g., iMantz et al.ll2008l ; IVikhlinin et al. 
2009h . 



A serious study of how all of these (and other) observations would be modified by our 
apparent acceleration mechanism goes far beyond the scope of this paper, and would likely be 
premature with these calculations being essentially in the toy model stage. But any eventual 
finding that our formalism does as well (or better) than simple Dark Energy models at 
achieving a concordance, based upon a variety of such measurements, would be convincing 
evidence in favor of our paradigm over the more physically exotic FLRW approach. 



6.2.2. Testing Our Formalism with Estimates of the Jerk Parameter, j$ 



Obs 



Given the importance of the jerk parameter in potentially falsifying ACDM (via jo ^ 1) 
and supporting our formalism (if jo significantly exceeds unity), it is useful to provide an 
estimate of the actual value of j here. Accurately measuring j is quite difficult, however, 
even with the best of today's Precision Cosmology data sets, as one is required to go to third 
order in the expansion of the luminosity distance (cf. Equation I2"6"j) for cosmographic studies; 
or equivalently, one must go to fits with four terms, (H ,Q\,Wq, w£), in dynamical studies 
with an evolving Dark Energy EoS, in order to calculate jo as a function of (Q\,Wq, w£). 
The difficulty of produci ng stable, pr e cise e stimates of j and higher-order cosmological 
parameters is evident in Ku fc Wangl ( 120101 ). in which they use their own cosmographic 
methods to obtain values like jo = — 4.996^7;3^ 3 and jo = 15.665lg3; 81 2, for two cases using 
different com binations of external da ta sets in conjunction with the recent SCP Union2 SNe 
compilation flAmanullah et al.ll2010f ). 
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Several previous estimates of jo by this author and by other researchers (details not 
given here), obtained from different SNe data sets and produced with a variety of methods, 
show a general tendency for j > 1 (ranging broadly from j ~ 0.9 — 5.5 in most cases); and 
indeed most data sets and analysis methods give values for j on the hi gh side - though ther e 



are important exceptions, such as the Constitution SNe compilation ( iHicken et al.ll2009al ) 



which prefers the rather low value of jo ~ 0.5. This latter result r eflects the significant 



statistical differences between the first SCP Union SNe compilation ( jKowalski et al.l 120081 ) 
and the Constitution set; for the former, most of the allowed parameter space (see their 
Figure 16) prefers th e Phantom Energy regime (wq < — 1), whereas for the latter, Figure 4 of 



Hicken et al.l ( )2009bl ) shows the Constitution set to be more centered within the Quintessence 
region (wq > —1) of DE parameter space. (Conditions of Wq < —1 tend to prefer jo > 1) 

A 




barring a reversal due to large values of of the wrong sign; and vice- versa for > — 1 



favoring j < 1.) 

The Constitution papers themselves ( Hicken et all2009a] jb1) do not an include an analysis 



of DE with a time- varying EoS; and here we will focus upon SNe from the original SCP Union 
compilation (our benchmark data set for cosmological model fits), along with the more recent 
SCP Union2 compilation (which supersedes the original Union and Constitution data sets), 
using a couple of straightforward methods for calculating j . 

As our theoretical motivation is to estimate the purely observable parameter j° hs devoid 
of any assumptions about a DE equation of state, the most direct, cosmographic approach 
is to do polynomial best-fits to the SNe luminosity distance values. Using the Taylor series 
expansion for di,(z 0hs ) (cf. Equation 126]) . it is straightforward to invert the coefficients from 
such a polynomial best-fit in order to obtain values for HQ hs , qQ hs (or equivalently WQ hs , cf. 
Equation [27b ) and jQ hs from the first three terms in the expansion. 



The main difficulty, as explored in a key paper (jCattoen fc Visserl 120081 ) examining 
earlier SNe data sets, is that such estimates of cosmological parameters (even for the second- 
order term qQ hs ) are "distressingly" unstable, depending very sensitively upon which redshift- 
distance relation one selects (the luminosity distance, d^, is not a unique choice), which 
redshift variable one uses, and how many terms one employs in the best-fit polynomial. This 
parameter instability (even when the polynomial expansions used all provide quantitatively 
very good fits) renders standard uncertainty estimates useless, as we will have to dismiss 
cases with 'unrealistic' parameters out of hand, based upon subjective prior judgments of 
likely ranges for vjQ hs and j^ 138 - 

First, considering the most standard case - polynomial fits to d L (z) - it is obvious from 
Equation [261 that we must include at least three terms (i.e., up to 0[(z) 3 ]) in order to obtain 
a value for jo- But that does not mean that we must use only three terms; in fact, we could 
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include any number of terms, N > 3, in our expansion, and use only the first three of them to 
obtain a value of j calculated independently of those higher-order terms (thus disregarding 
their weak information regarding higher-order cosmological observables) . The advantage of 
using N > 3 expansion terms is that this gives the best-fit polynomial additional terms for 
modeling the high- z behavior of the data, without forcing that constraining responsibility 
onto the third term from which we must calculate jo- As put by D. Rubin (2008, private 
communication), the N = 3 case "truncates the series expansion of a(t) after j , so even a 
LCDM universe will not give a fit of jo = 1"; a problem alleviated via N > 3. The downside 
of using too many terms, however, is that they become f ar too weakly constraine d by the 



data, making the effects of statistical uncertainties worse (ICattoen fc Visserl 120081 ) . so that 
even empirically good fits will yield more and more unrealistic cosmological parameters for 
larger and larger N. In practice, we perform polynomial fits using N = (3,4,5,6) terms, 
with the iV = (3, 4) cases usually producing the best parameter values (and often the only 
sensible ones). 

Applying these polynomial fits to d^z) for the original Union SNe compilation, we 
obtain (w$ hs , j$ bB ) = (-0.746,1.32) for the N = 3 case, and (w$ bs , j Obs ) = (-0.822,2.49) 
for the N = 4 case. Things appear to start going bad for the N = (5,6) cases (with 
(w? bs Jo° bs ) = (-0.625,-1.76) and (-0.918,7.19), respectively); but the N = (3,4) cases 
make a range of jQ bs ~ 1.3 — 2.5 (or more broadly, ~1 — 3) for the best- fit values seem 



quite reasonable. This is in decent accord with the varying-?/; analysis in ( iKowalski et al. 



20081 ). where the allowed parameter space depicted in their Figure 16 shows the ACDM point 
(i.e., (w ,w^) = (—1,0), equivalent to j = 1) to lie just around the inside edge of the lex 
statistical ellipse. 

Applying these same fits to d L (z) for the Union2 SNe compilation, one gets a startling 
result: the N = 3 case produces (w Obs ,j Obs ) = (-0.729,0.991), and the N = 4 case yields 
(w° hs , jo° bs ) = (—0.733,1.050). (With the N = (5,6) cases again being less physically 
reasonable, yielding (w$ hs , j$ hs ) = (-0.697,0.187) and (-0.419, -8.619), respectively.) Such 
estimates of j® bs = 0.991 and 1.050 would appear to be amazingly good for Cosmological 
Constant DE, seemingly virtually to prove ACDM. But in fact, while these results are most 
definitely consistent with ACDM, they are misleadingly precise - there is no way that the 
Union2 SNe data alone could reliably yield estimates of jQ hs to within 5% or less of the 
ACDM value, even if that model is true. Such values are impossibly accurate, requiring us 
to dig a little deeper in order to get some sense of the real parameter estimate uncertainties. 



Following ICattoen fc Visserl (120081 ) . we note that d-^(z) is not really a good function 
to use for a Taylor series expansion. First of all, the series will not even converge for 
z > 1 (a relevant redshift range for at least ~20 of the Union2 SNe), rendering the high-;? 
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behaviors of the best-fit polynomials meaningless. Rather, they recommend the use of the 
"2/-redshift" variable, y = z/(l + z) — [1 — (a/a )], which is bounded above by unity for 
past epochs - i.e., y G [0,1) for z G [0,oo). Second, in order to remove the dominating 
influence of the "nuisance" parameter, H® hs , which does not encode any of the dynamical 
cosmological information, they recommend instead fitting to a function such as In [d^y) jy\. 
Series expansions of this type of function have a first term (the only one containing the 
Hubble Constant) that is decoupled from the higher terms and contains no powers of y, 
and can thus be subtracted off as a constant offset. Lastly, they note that there is nothing 
necessarily unique or special about the luminosity distance function, g?l; and that several 
other cosmological distance functions - e.g., the "photon count distance" dp, the "angular 
diameter distance" d\, etc. - can be obtained from it simply by dividing by different powers 
of (1 + z) = 1/(1 — y). This is troublesome because fitting data with large scatter and 
uncertainties (such as the SNe data) to different distance functions will produce extremely 
different values of {qo,jo) ~ especially for fits to polynomials with a small number of terms, 
iV - and because (as they claim), "There is no good physics reason for preferring any one of 
these distance variables over the others." 

For our analysis here, we will consider fits of the Union2 SNe to \n[dp(y)/y], where 
dp = d L / (l + z) = d L (l — y). We use dp for three reasons. First, we argue that dp = (a r SN ) 
actually happens to be the most physically appropriate distance function, being simply equal 
to the present physical distance from the SN to us; it is the only distance function without 
an (artificial) factor of (1 + z) in its definition, and thus it will be finite for any observed 
physical object (even back to the Big Bang, assuming finite cosmological horizons), and 
goes to zero as z — > simply as the real distance to such an object would become zero. 
Se cond, as dp is in the midd le of the five physically-motivated distance functions considered 



by ICattoen fc Visserl (120081 ) . it yields the median values of their parameter estimates, and 
quoting the median estimates that one may find would seem to be a conservative strategy. 
Third, as a practical matter, we find that none of our polynomial fits to \n[d^(y)/y] yield 
reasonable values for j° bs (not to mention u>° bs ); whereas our N = 3 fit (though only that 
one) to ln[dp(y)/y] is more well-behaved. 

The result of this N = 3, \n[d P (y)/y] fit to the Union2 data set is (w° bs ,j° bs ) = 
(-0.750, 1.338). This value of j Obs is also consistent with ACDM, while being less unbe- 
lievably accurate, and gives us some (very qualitative) sense of the uncertainties involved. 
Simultaneously, it remains slightly high of ACDM, as our formalism predicts; and if not 
nearly so high as those values of jQ hs from our simulated cosmologies given above in Table [TJ 
the incorporation of recursive nonlinearities into our models will probably bring their cal- 
culated values of j° bs somewhat lower, as discussed previously. Nevertheless, Cosmological 
Constant ACDM remains entirely consistent with these results. 
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One other issue for these cosmological parameter fits relates to the "outlier rejection" 
process performed on these SNe data compilations, done as a final supernova cut, after all 
other cuts based upon the internal quality of the data points themselves have been com- 
pleted. This outlier rejection cut is the removal of SNe data points based upon their poor 
(>3cr) fits to flat ACDM reference cosmologies that are individually optimized and applied to 
each of the component d ata sets that go i nto th e Un ion and Union2 compilat ions. Extensive 



arguments are given in iKowalski et al.l (120081 ) and lAmanullah et al.l (120101 ) justifying this 



outlier rejection process as an important technique for robust statistical analysis; and yet, 
since the reference cosmologies used for the cuts are invariably the same type of flat ACDM 
models that are at the heart of the Cosmic Concordance Model, there is the unavoidable pos- 
sibility that such outlier rejection cuts may remove some legitimate evidence against ACDM, 
perhaps thereby introducing some pro- Concordance bias into the results. It is useful here 
to provide a simple demonstration of the extent to which cosmological parameter estimates 
can be dependent upon the decisions made regarding such data cuts. 

This 3a outlier rejection cut removed 8 SNe from the SCP Union compilation, reducing 
the number of SNe to 307 (i.e., culling ~2.5% of the data); and removed 12 SNe from the 
Union2 compilati o n, red ucing it to 557 SNe (i.e., culling ~2% of the data). (Furthermore, 



Amanullah et al.l (120101 ) applied a 5a 'outlier cut' earlier on in their analysis, which had 
already removed ~6% of their data; but it is not clear if this is the same type of cut as the 
cosmology-dependent outlier rejection done at the end.) Here we specifically consider the 
Union2 compilation with the 12 'outlier' SNe added back in (i.e., 569 SNe in total). Above, 
without these outliers, our JV = 3 fit to d L (z) yielded (w° bs , j° bs ) = (-0.729,0.991); adding 
the outliers back inE], this changes to (u>° bs , jo hs ) = (—0.736, 1.151). Similarly, the N = 4 
fit to d L (z) changes from (w? bs ,j Obs ) = (-0.733,1.050) to (w$ hs ,j$ hs ) = (-0.826,2.633). 
(Also, the N = 5 fit to d L (z) now gives sensible results, (w$ hs ,jg hs ) = (-0.775,1.351).) 
Meanwhile, the N = 3 fit to \n[d v {y)/y\ changes from (w° bs ,j° bs ) = (-0.750,1.338) to 
(w Obs , j° bs ) = (-0.810,2.422). Adding the 12 Union2 outliers back in thus moves the best- 
fit j° bs results much farther away from the ACDM value of jo = 1 than before (though 
under the circumstances, not representing a statistically convincing rejection). What this 
may ultimately prove is not immediately clear, since it is readily apparent from inspection 
of the wild scatter of the outlier SNe points that the quality of those measurements is 



u The twelve outlier points, like the rest of the Union2 SNe data, are publicly available at 
http://supernova.lbl.gov/Union; but the outlier points are only available in a table including all SNe with 
lower precision values and uncertainties, and with a different Hq normalization. Going to lower precision 
for the standard Union2 compilation of 557 SNe only changes the cosmological parameter estimates from 
polynomial fitting by ~0.3 — 5.5%, so for consistency, we do our entire Union2+Outliers analysis (569 SNe) 
using the lower-precision data table. 
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quite mixed, and several of them undoubtedly really are bad data points; but it does show 
how sensitive these cosmological parameter estimates can be to decisions made about data 
handling that are nearly invisible in the final quoted results. 

Now, these jQ hs estimates from polynomial fits to the SNe data alone, while interesting, 
undoubtedly have uncertainties large enough to seriously limit their usefulness. In order to 
get more precise estimates of jQ hs - and to be able to calculate some believable number for 
the statistical uncertainty on that value - we are compelled to combine the SNe data with 
complementary cosmological data sets (e.g., SN+CMB+BAO). It is not trivial to do this in 
a completely cosmographic manner, and so we therefore make the concession of considering 
jo in light of the evolving-EoS Dark Energy analyses done in the Union and Union2 papers. 

For these analyses, the SCP collaboration adopts the form: 

w\z) ee w£ + = w A + wty , (39) 



widely referred to a s the Chevallier-Polarski-Linder (CPL) parameterization ( jChevallier fc Polarski 



200ll ; iLinderi 120031 ) for the Dark Energy EoS. 



Using th e cosmic expansionary history that would result from such a Dark Energy 
(ILinderl 120031 ) to relate these CPL parameters to the cosmographic parameters {wQ hs , jQ hs ) 



obtained from the series expansion of d^(z) (cf. Equation [2HD, a straightforward calculation 
yields: 

w Obs = w A n A , (40) 

which makes obvious sense; and: 

JO ' 8 = {1 + [^0(1 + O] + [j^W]} , (41) 

which gives j^ hs = 1 for the Cosmological Constant case of (wq = — l,w A = 0), for any 
value of as required. Furthermore, note that Equation Hi] is an exact expression for jQ hs 
in terms of (fl A , w^, u> a ), given how they are all defined, even though it has been isolated by 
doing a comparison of various terms between two different series expansions. 



For use in applying this to the SCP Union compilation study in lKowalski et al.l ( 120081 ). 
D. Rubin states (2008, private communication) that their best-fit parameter values for the 
CPL Dark Energy EoS, in a combined analysis of Union SNe, CMB, and BAO data sets, 

; c . „„A _ 1 iq+0.15+0.21 j A _ n 7 o+0.53+0.67 -.i o — 1 O — fl ooo+0.018+0.021 
IS. W Q - -l-lc5_o.13-0.19 ancl w a ~ U - < <3-0.69-0.82> With M M — 1 - \l A - O._8__ 017 _ a020 , 

where each first set of uncertainties includes statistical errors only, and with each latter set 
including statistical plus systematic errors. 



Using these best-fit DE EoS values of (fi A , wt w t) = (0.718, -1.13, 0.73) in EquationiU 
one obtains ~ 2.26, if one simply plugs those values into this formula for the jerk 
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parameter without regard to the probability distributions of those parameters. This result 
for j® hs is fairly similar to the value of 2.49 obtained earlier from our N = 4 polynomial 
fit to this Union SNe data alone, it is demonstrably above the flat ACDM case, and is not 
all that far out of the lower end of the range of j° bs values quoted previously for our best 
simulation runs. 

Using the uncertainties for each of the parameters given above (averaging the ±a values 
in each case), in conjunction with Equation HU one could naively calculate la error bars for 
j° bs of Aj° bs ~ 0.87 (stat errors only), and Aj° bs ~ 1.14 (stat plus syst), actually placing 
ACDM more than la away in both cases. B ut calculating a valid uncertainty on j^ bs is not 



so simple, however: judging by Figure 16 of iKowalski et al.l ( 120081 ). the uncertainties on Wq 
and «; a are clearly not independent; and lacking covariance information, we cannot calculate 
a precise value of Aj^ bs . 

For the Union2 compilation, however, we have both more recent and more complete in- 
formation (D. Rubin, 2010, private communication). For their stat-errors-only analysis, they 
obtain best-fit values of: (Q\,Wq,w£) = (0.727,-1.046,0.160), with a (symmetric) covari- 
ance matrix of (caa, c qo , c aa , c\ Q , CA a , c 0a ) = (1.782 x 10~ 4 , 1.305 x 10~ 2 , 3.053 x 10 -1 , 1.766 x 
10 -4 , 7.254 x 10" 5 , -5.672 x 10~ 2 ); while their stat+syst analysis yields: (Q A ,w£,w£) = 
(0.723,-1.134,0.585), with a covariance matrix of (caa, coo, c aa , cao, c\ a , co a ) = (2.776 x 
10" 4 , 2.794 x 10~ 2 , 3.370 x KT\ -4.747 x 10~ 4 , 5.658 x 10~ 4 , -8.803 x HT 2 ). A straightfor- 
ward calculation (though ignoring the non-Gaussian nature of the parameter distributions, 
particularly for w£) produces final results of: J° bs ± Aj^ bs ~ 1.331 ±0.986 (stat errors only), 
and ~ 2.129 ± 1.292 (stat+syst). Clearly, the uncertainties on the jerk parameter remain 
quite large (with ACDM now back within la), even when using the Union2 SNe compilation 
data (without outliers) in conjunction with the CMB and BAO data sets. Furthermore, 
these two results differ from each other greatly - not only in the size of the error bars for 
j° bs , but also in its best-fit central value - with the only difference between them being 
whether or not one includes systematic uncertainties in the fitting process. 

Summarizing the numerical estimates from this subsection: different analyses of the 
original SCP Union data set yield a variety of values for the jerk parameter, including 
jObs ^ ^ 32 (AT = 3 case) and ~ 2.49 (N = 4 case) for polynomial fits to d L (z) for the 
SNe alone; and j^ bs ~ 2.26 is obtained from a non-cosmographic analysis of an evolving- 
DE EoS with (SN+CMB+BAO) data. For the SCP Union2 data set, alternatively, such 
polynomial fits - specifically, the A^ = 3 and A^ = 4 fits to d^z), and the A^ = 3 fit 
to ln[dp(y)/y], respectively - yield three different best-values for the jerk parameter of: 
j'q bs ~ (0.991, 1.050, 1.338) for the official 557 SNe Union2 compilation with outlier rejection, 
and j$ bs ~ (1.151,2.633,2.422) for the 569 SNe Union2 data set with the ACDM-outliers 
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added back in. Finally, an evolving-DE EoS analysis of Union2 (SN+CMB+BAO) data (557 
SNe) yields j^ bs ± A j° bs ~ 1.331 ±0.986 (stat errors only), and ~ 2.129± 1.292 (stat+syst). 

The conclusions that we may draw from these results can therefore be summed up as 
follows: (1) It seems likely that the jerk parameter lies within a (very approximate) range 
of, let us say, j° bs ~ 0.35 — 3.4; (2) Very large uncertainties remain, with estimates of j° bs 
remaining highly unstable and sensitive to which fitting functions and variables are used to 
obtain them (with many cosmological parameter fits having to be discarded on a purely ad- 
hoc basis), and are also very sensitive to whether or not systematic uncertainties are included 
in the analysis; (3) Various estimates of j° bs are systematically high (somewhat above unity), 
as is favored by our formalism; though the Cosmological Constant requirement of jo = 1 is 
within la of essentially all reasonable fits for the cosmological parameters, and going from 
the Union to the Union2 SNe compilation does appear to move the results somewhat more 
towards ACDM; (4) The range of jerk parameter values obtained from the ~dozen 'best' 
simulation runs with our formalism, j®^ ~ 2.5 — 5.5 (and j° bs ~ 2.6 — 3.8 for the 6 very best 
runs) is somewhat high compared to this most likely range of j° bs as determined from current 
observations; but the large scatter and uncertainties in the observational data, and the major 
theoretical simplifications in our calculations (e.g., the lack of recursive nonlinearities in our 
current models) may plausibly account for this. With further development of our models, 
and with enlarged and improved SNe data sets, it should be possible in the not-too-distant 
future to use observations of the cosmic jerk parameter to definitively support or rule out 
our causally-propagating perturbation formalism as a competitor to Concordance ACDM. 



7. THE COSMIC FATE: EXPANSION VERSUS DISORDER 

One of the signature issues to be addressed by any cosmological paradigm is the question 
of its predictions for the cosmic future: What is the ultimate fate of the universe? While 
the toy model presented in this paper for calculating the effects of causal backreaction is 
too simplified to produce a detailed quantitative answer to this question, it is useful enough 
for discussing the competing physical processes that are involved in determining the final 
cosmic fate. 

For standard FRW cosmologies - without any Dark Energy or Acceleration - the issue 
of the ultimate cosmic fat e is a theoretically sim ple one, depending only upon the total mass- 



energy density, f^Tot (e.g.. lKolb fc Turnerlll990l ). A closed universe (fiiot > 1) will recollapse 



in an 'infinitely' dense "Big Crunch"; whereas a universe with insufficient closure density will 
either expand eternally (open universe, fixot < 1) or asymptotically approach zero expansion 
speed (critical/flat universe, Q^ot — 1), in either case approaching a thermodynamic 'heat 
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death', with an absolute and final cooling of the universe into a "Big Freeze" (or "Big Chill"). 
Cosmic acceleration complicates the issue, making it dependent upon the nature of the Dark 
Energy. If the Dark Energy should be some sort of Quintessence that ultimately loses its 
negative-pressure properties and turns into a form of 'normal' matter, then the cosmic fate 
would revert to the standard case of depending solely upon fiiot- On the other hand, if 
the Dark Energy retains its accelerative potency - such as in the case of the Cosmological 
Constant, A - then the acceleration would never stop, and all matter not tied to t he final large 



clust er containing all bound mass in our region of space ( "Milkomeda" , e.g., ICox fc Loeb 



20081 ) would fly away out of causal contact; the total amount of information (and thus useful 
energy for doing work) available to any observer would therefore be finite even for an eternal 
universe, eventually causing the ability to do computations (an d thus the possibility of 



sustaining life) to ultimately fade away (IKrauss fc Starkmanll2000l ). Or even more severely. 



if the Dark Energy should actually increase in potency over time - such as would be the 
case for Phantom Energy (u> A < —1) - then the universe could end in a "Big Rip", leading 
to a Cosmic Doomsday where all cosmic structures, objects, and even atoms are ultimately 



ripped apart (ICaldwell et al.ll2003l ). Other exotic mo dels naturally lea d to other possibilities 



such as a cyclic universe with a "Big Bounce" (e.g.. IVenezianol 12004 ) . and so on. 



For the causal backreaction formalism that we present here, in which exotic contrib- 
utors such as Dark Energy, Higher Dimensions, etc., are unnecessary, the situation is less 
complex in theory while being more complicated in practice. While there is no more need 
to understand the behavior of any independent substance or modified physics that exerts 
some 'external' control over the rest of the universe, there is now the necessity of calculating 
a fully self-consistent solution for the evolution of an increasingly inhomogeneous cosmos, 
whose own perturbations are causing itself to accelerate (or at least to appear to do so). A 
real solution of this type would require a method of dealing with gravitational nonlinearities, 
along with fully three-dimensional modeling to deal with the ultimate breakdown of our 
smoothly-inhomogeneous approximation. But even without such tools at our disposal here 
we can at least delineate some of the more likely possibilities, to see how they compare to 
the popularly regarded cosmic fates mentioned above, and also to see if any new wrinkle 
appears that modifies these possibilities in some interesting way. 

In terms of our clumping-perturbed metric, Equation [121 the effects of causal backre- 
action completely take over the cosmic evolution as I(t) approaches unity; and the model 
breaks down for I(t) ~ 1 (and is completely meaningless for I(t) > 1), as gravitational 
nonlinearities cause the failure of our approximation of linearly summing together the in- 
dividual Newtonian gravitational potentials from different virialized mass clusters. Two 
questions naturally arise: (1) Will this I(t) function ever actually approach unity? (2) And 
if so, what really happens to the cosmic expansion as a result? 
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Considering the latter question first, to see what the stakes are, it seems likely that 
'real' acceleration due to backreaction, in the sense of qi from Section [2] above, would finally 
become cosmologically important. In other words, the fitting problem of finding an averaged 
cosmology to represent a very inhomogeneous universe would become a dominant concern, 
and accelerated volume expansion due to the continuing causal backreaction, now modified 
by truly nonlinear gravitational effects, would fully control the evolution of the universe. 
In that case, we suppose that it is possible for the result to be a runaway acceleration, as 
potent perhaps as that from a Cosmological Constant (or maybe even as strong as a form 
of Phantom Energy), since the effective value of w° hs would be due to gravitational energy 
rather than due to the exotic nature of any particular physical substance, and hence not be 
subject to the bounds (i.e., w > — 1) normally imposed upon perfect fluids by the dominant 
energy condition. (One could even speculate that such a process may have had relevance 
to the pre-inflation very early universe - given its presumably chaotic and inhomogeneous 
nature - perhaps accounting for some or all of the acceleration usually credited to the Inflaton 
field(s); assuming, of course, that such a process could be shown to be theoretically possible 
given causality restrictions.) 

Alternatively, some authors have argued that backreaction due to matter inhomo- 
geneities cannot be responsible for A- or Phantom-like acceleratio n; or that i f they can 



be, temporarily, then such an acceleration cannot continue eternally. iKasai et al.l (120061 ) use 
the volumetric conservation of pressureless dust to make arguments towards a no-go theorem 
for acceleration from nonlinear backreaction; but their calculations are based upon the same 
cosmologically inapplicable approximations discussed previously in Section [2J Also, they use 
the same function a(t) as both their unperturbed scale factor, and as their measure of cosmic 
volume (i.e., a oc V 1 ^ 3 ) in the perturbed metric - two roles for a(t) that are not simulta- 
neously valid in a model such as ours, for example, since I(t),I(t) ^ alters the evolving 
physical volume of a perturbed region away from a simple proportionality to a 3 (e.g., recall 
Section [5.21) . 



Bose fc Majumdarl (120111 ). on the other hand, argue that whatever specific physical 
process may actually be generating the current cosmic acceleration, that it cannot continue 
eternally, due to backreaction from the future cosmological event horizon which would ex- 
ist in a universe that a ccelerated forever. Despite the use of the Buchert formalism (e.g., 



Buchert fc Ehlersl 119971 ) - along with its unacceptable assumptions for causal backreaction 
(irrotational dust, zero Newtonian backreaction, etc.) - in their calculations of the back- 
reaction from this future horizon, their results are obtained without necessarily presuming 
Buchert-generated acceleration; and they claim that the acceleration must eventually stop 
regardless of the mechanism behind it. Looking closely at their results, however, one sees 
that this claim only really works when one considers a large region made up of subregions 
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(two, in their calculations) that are each individually decelerating; and in that case, all they 
have proven is that the Buchert-generated acceleration of the combined region eventually 
becomes weaker than the total sum of the individual decelerations of the component subre- 
gions. If one of the subregions is actually accelerating on its own through some non-Buchert 
mechanism (such as via our causal backreaction method, nonzero A, etc.), and possesses a 
sufficiently large initial fraction of the total volume - or if all subregions are individually 
accelerating (a > 1, /3 > 1 in their terminology) - then inspection of their Equations 4 and 
9 shows that the overall acceleration do cannot go to zero for positive, finite t, unless the 
'apparent volume fraction' of a component subregion as part of the total were to become 
negative, which is obviously not consistent. 

While the conclusions that one may draw from calculations such as these are debatable, 
and the ideas behind them are still in too early a stage of development to be definitive, they 
still make the important point that there may be 'global' considerations that either prevent 
or ultimately shut down a long-term cosmic acceleration due to effects that go beyond those 
contained within any simplified model like our formalism. It is therefore important to keep 
such caveats in mind when considering the so-called 'ultimate' possible fates of the universe, 
below; but regardless of the eternal validity of any conclusions, at least we can use our model 
to elucidate some reasonably long-term effects of the crucial phase transition which must 
undoubtedly occur when (and if) the causal backreaction becomes nonlinear in strength. 

Now considering our former question - will I(t) really approach unity? And even more 
seriously, will gravitationally nonlinear effects due to causal backreaction come to dominate 
(and in some sense 'permanently' control) the evolution of the Universe? In our simulation 
runs, I(t) increases monotonically, and gets quite close to 1 for the strongest-'accelerating' 
cosmo logical models; it would therefore be useful to estimate how large I(t) can get within 
our formalism. 

In the essential physics of backreaction with causal updating, there are several different 
processes competing against one another - some enhancing the apparent (or ultimately real) 
acceleration, and some damping it - and which way the balance tilts between these factors 
will largely (but as we will see in the end, not completely) determine the ultimate cosmic 
fate. The most obvious factor enhancing the acceleration is the increasing 'dumpiness' of 
the universe, represented in this paper by clumping evolution function ^(t), which increases 
monotonically over time. In our formalism, however, the extent of clumping must saturate 
as \l/(t) — > 1, thus limiting its ability to continue driving the growth of I{t); and this seems 
like a physically realistic constraint on ^(t), representing how virtually all cosmic matter 
should eventually reach a 'maximally clumped' state, in which the universe consists almost 
entirely of steadily expanding voids, punctuated only by a sparse scattering of fixed-size 
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massive clumps that no longer expand with the universe, but which manage to virially 
support themselves against singular collapse. 

But even more important than the increasing dumpiness of any given region of space, 
is how the information about distant dumpiness (i.e., \l/[£ ret (£, a)}) can get to an observer 
from farther and farther distances with increasing t, as the volume of space with significant 
clumping within the observer's causal horizon (a < a max , cf. Equation [TOl) grows rapidly 
with time. This effect is counterbalanced, however, by the continuing expansion of the 
universe, which dilutes such backreaction by taking any given mass 'clump' farther away 
from the observer as t increases (specifically referring to the factor of 1/[£Imd(£)] i n -Rsch(^); cf. 
Equations IBlfTTj) . These opposing factors are represented by a max and (to/t) 2 ^ 3 , respectively, 
in Equation [T2l for I(t). Ignoring the time dependence in [t vet (t, a)} for the moment, one 
approximately gets I(t) oc [(t /t) 2 ^ 3 a^ ax ]; but since a max ~ (t/t ) 1//3 for t £ init (cf. 
Equation [TOl) . we see that the overall time dependence ultimately cancels out of I(t) at 
sufficiently large t. We can interpret this to mean that the increasing backreaction effect 
due to the expanding causal horizon also saturates, due to competition against the cosmic 
expansion itself, with I(t) approaching a constant value at times very late compared to ^ init 
(and very late also compared to the time by which the clumping evolution, ^(t), has mostly 
saturated near unity). What we need to know, of course, is how large this asymptotically 
constant value of I(t) will be; an obviously model-dependent question. 

Treating \I/[i ret (t, a)] as a step function - equal to ~1 for a < a max (with t init ~ 0), 
and zero beyond there - one gets I (to) ~ 6 ^> 1, a result deep into the regime requiring 
a gravitationally nonlinear treatment; but this of course is an unrealistic upper bound on 
the effects of clumping, representing a universe that became fully clumped almost instanta- 
neously, saturating at very early cosmic times (i.e., ^(t) — > 1 at it/to) <C 1). Trying more 
realistic (and continuous) functional forms for ^(t), with a fairly recent saturation of clump- 
ing - in particular, considering the three functions used for our clumping evolution models, 
as defined above in Equations |22r-c - the maximum possible value of I(t) for each clumping 
model can be obtained for t init ~ 0, t — > t , and ^(t Q ) = 1. Doing this, it is easy to show 
that J max (^ S qr) = 3/14, J max (^ L in) = 3/5, and J max (^ M D) = 1- Thus the 'strongest' of our 
clumping evolution models ($md) can definitely produce results strong enough to violate 
the approximation of gravitational linearity; and even stronger models, with more gradual 
clumping growth (that is, an earlier start to significantly large \1/) - given as \l/(t) oc (t/to) N 
with iV < 2/3 - would have J max values that easily exceed unity (e.g., J max = 2 for iV = 1/3, 
/ max = 6 for N = 0, etc.). In short, causal backreaction as prescribed by our formalism 
can clearly be strong enough to not only induce an apparent acceleration, but also to break 
down our linearized gravitational approximation entirely, indicating a possible runaway ac- 
celeration for the real universe, itself. But whether that possibility is realized, or whether 
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the strongest effects of causal backreaction fall short of it, depends extremely sensitively 
upon the precise time-dependence of the evolution of mass clumping and virialization in the 
universe. 

One important moderating factor for backreaction has been completely left out of these 
estimations, however: the recursive nonlinearities discussed previously, representing how 
the metric perturbations due to early inhomogeneities will slow down the causal updating 
that brings in new inhomogeneity information. In addition, some extra volume expansion is 
created (cf. Equation [S] and subsequent discussion) that adds to the ordinary FRW Hubble 
expansion in pulling clumped masses farther apart from one another, thus diluting these 
backreaction effects to a degree greater than that predicted by Equation [12] and numerically 
estimated above. These effects should not only reduce the value of j° bs predicted by our 
simulations (as mentioned in Section I6.2.2p . but may in a larger sense serve as a causal 
backreaction 'regulator', feeding back upon that very backreaction to slow down its effects 
upon the universe whenever the acceleration caused by it begins going into a runaway mode. 
But such a picture of self-regulating cosmic acceleration remains purely speculative until a 
more complete formalism and simulation program exist which can quantitatively model the 
effects of these recursive nonlinearities into the far future. 

Thus the variety of 'ultimate' cosmic fates potentially emerging from our causal back- 
reaction paradigm seem to fall into a few fairly recognizable categories, depending in detail 
upon the evolution of the growth rate of self-stabilized cosmic structures. On the one hand, 
the effects of causal backreaction may weaken or even saturate, leading either to a continuing 
low level of 'acceleration', or to the kind of coasting/mildly slowing behavior of open/flat 
universes - in either case resulting in a similar kind of Big Freeze as would have occurred 
in the normal (f^Tot < 1) FRW case. Alternatively, if the cosmic structure formation evolu- 
tion leads to a stronger, self-sustaining level of backreaction, then it is conceivable that the 
fading-away behavior (the "Big Fade"?) of the rest of the universe as would occur for an 
exponentially-expanding Cosmological Constant universe may not be a bad approximation. 
Or in extreme cases (probably less likely), some form of Big Rip could conceivably occur. 
Qualitatively speaking though, the most reasonable scenario in our paradigm may be some 
kind of stop-and-go accelerating behavior, due to the self-regulating nature of causal back- 
reaction with recursive nonlinearities. This could lead to a situation where the universe is 
always 'riding the edge' of a runaway acceleration, but never quite getting there - perhaps 
in a permanent state of acceleration that oscillates over time but stays relatively close to 
the level of acceleration that we see today, with the same relationship that we currently see 
between effective FLRW observables (i.e., f2^ bs ~ ^M bs t° within an order of magnitude), 
with these conditions remaining generically true for the foreseeable cosmological future. 
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Now, the variety of possible cosmic fates considered here is of course not unique to 
our formalism - each of these possibilities can easily be reproduced by a version of Dark 
Energy with some appropriately evolving equation of state. The difference is that instead 
of being dependent upon the unknown (and perhaps ultimately unknowable, in fine detail) 
theoretical nature of a master Dark Energy field running the rest of the universe, the cosmic 
fate here is due to the normal gravitational forces and motions of matter, which - though very 
complicated in practice - will depend only upon calculations using known laws of physics 
(once the nature of the Dark Matter is known!), and the long-explored mechanisms of general 
relativity. 

And yet, one more wild card still remains, one more consideration regarding the plausi- 
bility of these possible cosmic fates, which will tell us which one really is the most probable 
description of the cosmic future; and the answer is: none of them! In the entire discus- 
sion of our formalism in this paper, there is one 'elephant in the room' that has not been 
analyzed, because it cannot be modeled in this way - and that is the breakdown of our 
smoothly-inhomogeneous approximation itself. As the universe grows older, and ever-larger 
bound structures have time to form, the universe will eventually become inhomogeneous 
and anisotropic on a truly cosmological scale. Thus the Cosmological Principle will break 
down entirely as a usable approximation, and vast cosmic structures will eventually ex- 
ert powerful tidal forces that are strong enough even to threaten bound galactic clusters. 
While this is also not a theoretical phenomenon unique to our specific formalism, once again 
there is a key difference: whereas typically one expects the dominant perturbing influences 
on local structures to be the closer ones (cosmologically speaking), consisting of just one 
or a few definable, huge objects, our causal updating analysis shows that the acceleration 
effects due to backreaction are dominated by the farthest distance out to which one can 
causally see clumped structure. Thus in the distant future, the gravitationally perturbative 
effects acting upon otherwise 'local' structures will be dominated by a multitude of different 
inhomogeneities, located in all directions, that are cosmologically-extreme distances away, 
and that together represent an almost unimaginable amount of mass in clumped structures. 
Such tidal or shearing forces will continue to grow and grow - as long as the causal flow 
of information can continue to come in from ever-more-distant parts of the universe - until 
gravitationally bound local structures are chaotically ripped apart. In far more anisotropic 
fashion than the Big Rip of Phantom Energy - if perhaps less destructively, since some 
form of recursive nonlinearities and dilution will eventually slow the process - this evolu- 
tion will herald the end of the smooth cosmic expansion existing at present, rendering all 
FRW/FLRW models obsolete. The first hints of such behavior in our universe may already 
be making themse lves evident, perhaps accoun ting for the observed phenomenon known as 



the "Dark Flow" ( IKashlinsky et al.ll2008l . I201CH ) - or at least accounting for the causal part 
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of it, if the effect indeed extends to 'superhorizon' scales. One thing is clear, according to 
our causal backreaction formalism: this Dark Flow (assuming it to be observationally real 
and at least partly causally generated) is just the beginning of a greater trend of inhomo- 
geneity and anisotropy that will be acting upon all local cosmic structures, and doing so 
much more rapidly, violently, and chaotically than would be expected in more conventional 
(i.e., non-causal, Poisson-Equation-based) models of structure formation. 

The real fate of the universe will therefore not be a Big Crunch, a Big Freeze, a Big 
Fade, a Big Rip, or even a Big Bounce - instead it will end (if one can call it an 'end') in 
something entirely more appropriate to the entropy-riddled universe that we all have known: 
a Big Mess. 

8. SUMMARY AND CONCLUSIONS 

We may sum up the motivations, methods, and results of this paper as follows. Our 
starting point is the search for an alternative to Dark Energy as an ingredient in the nascent 
Cosmic Concordance, without resorting to the increased theoretical complexity represented 
by modifications to gravity, or the non-Copernican 'specialness' implied by a local cosmic 
void. The problems with Dark Energy are well known, including the coincidence and mag- 
nitude fine-tuning problems for a pure-A Cosmological Constant, and the stability problem 
for a more dynamical form of Dark Energy possessing the ability to cluster spatially. The 
latter (DDE) case requires the ad-hoc addition of some form of new pressure term to support 
it against collapse, given the locally attractive nature of the negative pressure required by 
the DDE to power the cosmic acceleration. Such a new pressure term, if adiabatic (e.g., 
degeneracy pressure), would represent a form of positive pressure contributing to a cosmic 
deceleration that partially or totally nullifies the acceleration meant to come from the DDE; 
and if non-adiabatic (due to some new effects from the DDE Lagrangian), would invalidate 
the Dark Energy as a perfect fluid, thus calling the usual (accelerative) cosmic implications 
of w DE < in the FLRW acceleration equation into question. 

Emphasizing the importance of the fact that 'negative pressure' is locally attractive in 
character - rather than repulsive, as it is often regarded and popularly described - and that 
normal gravitational attraction therefore represents (at least in a non-technical sense) a form 
of negative pressure, we have therefore made a virtue of necessity by recruiting cosmological 
structure formation itself, based upon ordinary gravitational forces, as the driver of the 
observed (possibly apparent) acceleration. This approach, known generally in the literature 
as "backreaction" , removes all of the aforementioned problems by both eliminating the need 
for any Dark Energy species, and by solving the Coincidence Problem by naturally linking 
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the onset of cosmic acceleration with the emergence of cosmic structure, which inevitably 
leads to the creation of observers such as ourselves just in time to see this 'coincidence'. 

Noting that several versions of backreaction have been largely unsuccessful in their 
attempts to account for the observed acceleration, and that powerful arguments have been 
advanced which claim that backreaction as a paradigm cannot be made strong enough to 
succeed at this task, we have described how each of these arguments (and the backreaction 
mechanisms which they address) are functionally invalid. This lack of validity is not generally 
due to mathematical flaws within the arguments themselves, but due to their inapplicability 
to the real universe: all such "no-go" arguments or theorems are based upon one or more 
essential simplifications so restrictive that they eliminate from consideration all of the actual 
physics that is responsible for the crucial 'accelerative' effects. 

Different arguments against backreaction fail for different reasons. No-go conclusions 
drawn from the Raychaudhuri equation are invalid because they neglect vorticity (or more 
precisely, the square of the vorticity, u 2 ), thus ironically dropping the dominant physical force 
in the universe which prevents the collapse of all (non-solid) structures into singularities, and 
which is a positive semi-definite quantity that does not go to zero regardless of how large 
a scale one uses for averaging it. Swiss-Cheese models (typically using Schwarzschild or 
LTB metrics for the holes) underestimate backreaction effects not only because they lack a 
mechanism for modeling vorticity and virialization, but more crucially because they partition 
the cluster-forming universe into a discrete set of non-intersecting, non-interacting volumes 
- a construct completely at odds with how structure formation works in the real universe, 
where overdensities extend their reach far beyond their local domains to cause vast inflows, 
and where any given region of space feels the gravitational effects of many such influences 
combined together, rather than being impacted by one influence while being hermetically 
sealed off from all of the others. Lastly, there are no-go arguments based upon the generally 
Newtonian nature of metric perturbations and mass flows, apparently demonstrating that the 
sum total effects of backreaction must be small (or according to the formalism of Buchert and 
collaborators, identically zero); but such arguments are flawed due to the non-causal basis of 
the calculations (since they drop all signal travel-time delays, gravitomagnetic effects, and 
time derivatives of the perturbation potential), and because they depend upon 'smallness' 
arguments that are irrelevant in the face of the combined perturbative effects of innumerable 
cluster-forming masses acting upon every single region of space within their causal grasp. 

The true utility of such no-go arguments, however, is that gaining an understanding 
of why and where they fail to apply can guide us in the construction of a more successful, 
working form of backreaction. The central challenge is figuring out how best to model 
the crucial cosmic phase transition from smooth matter distributed largely homogeneously 
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throughout space, to most matter being concentrated into a fairly randomly-distributed 
collection of vorticity-stabilized, self-virialized, clumped structures. In this paper, we model 
the situation heuristically and simply by representing the dominant perturbation due to each 
'clump' as being (the Newtonian tail of) a Schwarzschild metric embedded in the expanding 
universe, superposed on top of the background FRW metric - where we utilize the fact 
that the vorticity (or velocity dispersion) stabilizing each structure will add its own positive 
contribution to the volume expansion, as per the Raychaudhuri equation. We then combine 
the perturbative contributions of all such objects by summing them and averaging over 
location and direction, producing a "smoothly-inhomogeneous" universe model that neither 
possesses nor requires a local void or spatial variations. Instead, it possesses only a single 
clumping evolution function that increases over time, representing the ever-growing fraction 
of cosmic matter located in self-stabilized clusters, rather than remaining in the smooth 
background. Lastly, and perhaps most importantly, we introduce the mechanism of "causal 
updating" in order to represent the crucial process of how the information about this growth 
of clustering propagates through the universe to any given observer, coming in from the 
very limits of the causal horizon out to which that observer has had time to see clumped 
structure. 

Employing these physical considerations, we have constructed a formalism of "causal 
backreaction" for calculating how the smoothly-inhomogeneous metric evolves over time, 
given any pre-specified clumping evolution function ^(t). From that metric solution, we can 
calculate a luminosity distance curve as a function of observed redshift z° hs , for comparison 
to residual Hubble curves measured via standard candles (particularly Type la supernovae), 
in order to check if an observed (possibly apparent) acceleration of the right magnitude and 
z° bs -dependence has been observed. Additionally, we can use our formalism to calculate 
several key cosmological observables produced by that function ^(t) - such as t® hs (age of 
the universe), ^m RW (matter density as a fraction of the unperturbed FRW critical density), 
WQ hs (apparent Dark Energy equation of state), j° bs (cosmic jerk parameter), and /° bs 
(characteristic angular scale of the CMB acoustic peaks) - testing whether the essential 
results of an alternative cosmic concordance have been achieved for causal backreaction 
operating upon such a clumping evolution function. Altogether, this adds up to a sterner 
test than the mere production of an apparent acceleration. 

Armed with this formalism, and using a set of 60 clumping evolution functions designed 
for simplicity and using model input parameters determined from straightforward astrophys- 
ical considerations, we have found a ~dozen-plus solutions that fit the SCP Union Compi- 
lation SNe data essentially as well as the best-fit ACDM model. Furthermore, about half of 
these models give good values (within reasonable theoretical and observational uncertainties) 
for all other calculated cosmological parameters. That is, these few 'best' models (obtained 
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without even conducting a full optimization search over the model parameter space) have 
successfully produced an alternative concordance, at least to the extent of solving the cosmic 
Age Problem, achieving spatial flatness (for the unperturbed, pre-clumping FRW universe), 
and being reasonably consistent with the CMB peak positions; and they achieve all of these 
goals - along with explaining the observed acceleration - in a matter-only universe, without 
any form of Dark Energy component. 

Now given that these important objectives - particularly the difficult one of achieving 
flatness by bringing the total cosmic density from Q ~ 0.3 up to Qq ~ 1 - have been major 
pillars in the argument for Dark Energy, our ability to achieve them with matter alone, 
simply by including the ordinary gravitational effects of causal backreaction, would seem 
to turn all forms of Dark Energy (A, Quintessence, etc.) into excess theoretical baggage. 
Nevertheless, since one cannot prove a theory with Occam's Razor alone, we have sought in 
this paper to find ways to distinguish our causal backreaction paradigm from any kind of 
Dark Energy (where possible), or at the very least to distinguish it from the simplest form 
of Dark Energy - i.e., a Cosmological Constant. One generic implication of using backre- 
action to generate the cosmic acceleration is that this places an increased emphasis upon 
the possible existence of larger-scale anisotropics and inhomogeneities than those expected 
from the standard perturbed-FLRW cosmologies. To that end, we have commented above 
about potential observational signs of such deviations from the Cosmological Principle that 
have been seen over the past few years by other researchers, as well as referring to work by 
this author and others using Type la SNe data sets to search for direct signals of possible 
large-scale anisotropics. 

Additionally, there are potential discrepancies for ACDM in particular, including the 
cuspy CDM halo problem and the possible dearth of dwarf satellite galaxies, which might be 
naturally explainable due to clustering-induced feedback occurring in the causal backreaction 
paradigm. Finally, there would of course be alterations to virtually all other cosmological 
observables due to causal backreaction, providing any number of possible signals in the com- 
plementary cosmic data sets that can be used for testing our paradigm; though detailed 
calculations of most of those alterations are beyond the scope of this current paper. One 
altered parameter that we have in considered in detail here, however, is the jerk parameter 
j° bs , which for flat ACDM is required to be exactly unity. For our cosmologically successful 
simulation runs using the causal backreaction model, on the other hand, j° bs > 1 is clearly 
preferred, a general result which we have seen is supported reasonably well (though statisti- 
cally very weakly) by trends in current cosmological data. Quantitatively speaking, the best 
runs in this paper predict an approximate range of j° bs ~ 2.5 — 5.5 (with a more narrow 
range of j'° bs ~ 2.6 — 3.8 for the very best few runs), which is perhaps high compared even 
to the upper range of observational estimates; but as explained previously, such values com- 
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puted here from our simulation runs are almost certainly higher than they should be due to 
current oversimplifications of our model (e.g., the absence of what we have termed "recursive 
nonlinearities" , etc.), which must be remedied with further theoretical development of these 
calculations. 

Having obtained such favorable results for producing an alternative concordance with- 
out Dark Energy, we then considered the 'ultimate' cosmic fates possible for a universe 
dominated by causal backreaction. Though the mathematical toy model introduced in this 
paper is still too simple at this point to yield detailed quantitative predictions about the 
universe in the far future, we expect the possibilities to be quite similar to the usual variety 
of potential fates for a universe dominated by various forms of Dark Energy with its typical 
designer equations of state. But there are two major differences, however: (1) With the 
acceleration (apparent or real) being provided by the detailed processes of structure for- 
mation in the matter itself, rather than from some 'external' material component with its 
own independent equation of state evolution, it seems far more likely that a self-limiting, 
stop-and-go form of acceleration will be the ultimate cosmic fate, rather than any complete 
runaway scenario (such as the Big Rip of Phantom Energy); and: (2) Given the unwonted 
importance here of inhomogeneities at great cosmological distances for determining 'local' 
gravitational behavior, causal backreaction is likely to lead to far more extreme and asym- 
metrical tidal forces on locally bound objects, eventually leading to a far messier and more 
chaotic universe than anything expected from 'normal' models of inhomogeneity evolution, 
given any form of cosmic acceleration this side of a Big Rip. 

Expressed in one sentence, if one asks, "What is the force behind the cosmic accelera- 
tion?" , the answer we would give is that it is not a 'force' at all: rather, the motivating effect 
may be called "the Schwarz" - consisting of the total, summed effects of the Newtonian tails 
of Schwarzschild-like metric perturbations, produced by the virialization of innumerable self- 
stabilized structures filling the universe, with these influences propagating causally towards 
all observers from the extreme edges of their observable cosmic horizons. 

I am grateful to Jacob Bekenstein, Krzysztof Bolejko, Varoujan Gorjian, Wayne Hu, 
Marek Kowalski, Edvard Mortsell, David Rubin, Ran Sivron, David Wiltshire, and Ned 
Wright for brief but helpful communications; and I am especially grateful to Arthur Lue for 
several helpful and clarifying discussions. 
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Fig. 1. — Simplified example of a physically permissible (though inelastic and nonadiabatic) 
merger of three nonrotating objects (a) into two rapidly spinning objects (b). (Not shown: 
ejected material containing mass and energy but zero total angular momentum.) Here, the 
angular momentum averaged over the volume of the combined system, (L), remains ~0, 
though (L 2 ) goes from psO to a large value. 




(r,0 = «>./,) 



Fig. 2. — Geometry for computing the inhomogeneity-perturbed metric at each point along 
the integrated path of a light ray from a supernova to our observation point at Earth. 
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Fig. 3. — Residual Hubble diagrams computed from our apparent acceleration model using 
^Lin(^) clumping functions. Different panels represent initial-clumping epochs of (respec- 
tively): z iQit = (5,10,15,25). Each plot depicts the usual flat ACDM (Q A = 0.73) and 
SCDM cosmologies, shown versus five of our simulation runs using present-day clumping 
parameters of (plotted highest to lowest): ^ = (1.0,0.96,0.92,0.85,0.78). 
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Fig. 4. — Residual Hubble diagrams computed from our apparent acceleration model using 
^Sqr(^) clumping functions. Different panels represent initial-clumping epochs of (respec- 
tively): z iQit = (5,10,15,25). Each plot depicts the usual flat ACDM (Q A = 0.73) and 
SCDM cosmologies, shown versus five of our simulation runs using present-day clumping 
parameters of (plotted highest to lowest): ^ = (1.0,0.96,0.92,0.85,0.78). 
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Fig. 5. — Residual Hubble diagrams computed from our apparent acceleration model using 
^md(^) clumping functions. Different panels represent initial-clumping epochs of (respec- 
tively): z init = (5, 10, 15, 25). Each plot depicts the usual flat ACDM (Sl A = 0.73) and SCDM 
cosmologies, shown versus five of our simulation runs using present-day clumping parameters 
of (plotted highest to lowest): ^ = (1.0,0.96,0.92,0.85,0.78). 
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Fig. 6. — Residual Hubble diagrams for the thirteen 'best' simulation runs, as described in 
the text, of our inhomogeneity-perturbed apparent acceleration model, shown along with 
the best-fit flat SCDM and Concordance ACDM = 0.713) cosmologies. From highest to 
lowest (at 2;° bs ~ 10), the plotted ^mdW curves have the parameters: (\&n> 2init) = (0.96,5), 
(0.92,5), (0.85,5), (0.78,10), (0.78,15), (0.768,14), (0.78,25). From highest to lowest (at 
^obs ^ -[g^ the plotted ^un(t) curves have the parameters: (\l>o ; z- mit ) = (1.0, 10), (1.0, 15), 
(0.96, 15), (1.0,25), (0.96,25), (0.92,25). Shown along with these curves are the SCP Union 
SNe data, here binned and averaged for visual clarity (bin size ALog 10 [l + ^° bs ] = 0.01). 
Each theoretical model is individually optimized in to minimize its Xmt f° r the SCP 
Union SNe data set (see Tabled]). For simplicity, instead of moving the SNe data up or down 
for each different optimized if° bs value, the optimization is depicted here by plotting the 
residual Hubble diagram of the SNe data versus a coasting universe of a single, fixed Hubble 
constant (HQ bs = 72 km s _1 Mpc _1 ), and then displacing each theoretical curve vertically, 
relative to the SNe data, as appropriate for each fit. 
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Table 1: Cosmological Parameters Derived from the 'Best' Runs of our Set of Simulations 





Xlit a 


TJf b 

"Fit 


/o c 
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0.57 


1.12 


69.95 
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1.037 


-0.817 


3.45 


284.2 


(1.0,15) 
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0.55 


1.12 
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41.63 


13.45 
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-0.784 


3.11 


287.6 


(1.0,10) 
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0.314 


0.52 


1.12 


69.27 


42.80 
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-0.746 


2.76 
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(0.96,25) 
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0.55 


1.11 
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41.57 


13.38 


0.968 


-0.759 


2.92 
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(0.96,15) 
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$md Clumping Model Runs 
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0.60 
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(0.96,5) 


315.5 


0.313 


0.63 


1.14 
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1.248 
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Semi- Optimized 
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■= 0.713 


-= i - n 
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311.9 


0.380 
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69.96 
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-0.713 
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Comparison Values from Best-Fit k flat SCDM Model (Q^ — 0, Qm — 1) 




608.2 


3.4E-22 




1.0 


61.35 


61.35 


10.62 
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0.0 
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287.3 



a x|;t computed versus the SCP Union SNe data (jKowalski et al.l 120081 ). As discussed in the text, the $md 
and Vf'Lin runs used here for study are not absolutely optimized (i.e., Xfh minimized) with respect to their 
model parameters, (^0i ^init)i as the ACDM model quoted here has been fully optimized with respect to Q\. 

& Each likelihood probability Ppit is derived from the corresponding Xfh using the Xn d f distribution with 
TVdoF degrees of freedom, where -/VdoF = 304 for our 'J'Lin and $md clumping models, TVdoF = 305 for the 
flat ACDM model, and AW = 306 for flat SCDM. 

c The integrated (Newtonian) gravitational perturbation potential at to, as computed via Equation 1121 
d Each z s quoted here corresponds to z FRW = 1. 

e The H® hs value (given here in km s _1 Mpc _1 ) for each run is found by minimizing its Xfh with respect to 
the SCP Union SNe data. 
•^Each Hq rw is computed relative to the corresponding optimized HQ hs value for that run. 
9A11 t° bs values arc listed here in GYr, and computed assuming no radiation (i.e., VIr = 0). 
h Al\ f2™ w values given here for the ^Lin and $md models are nor malized to n^ s = 0.27. 
I This "Semi-Optimized" $md run is one chosen from a class of low-x| it models, as described in the text. 
J "Best-Fit" for the flat ACDM Model refers here to an optimization over f^A and HQ hs . 
k "Best-Fit" for the flat SCDM Model refers here to an optimization over H® hs . 



